AD-A227 


OTIC  FILE  CWi 


UNCLASSIFIED 
CURlTY  CLASSIFICATION 


3> 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
OMB  H O.  0704-0198 


REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

It.  SECURITY  CLASSIFICATION  AUTHORITY 


b.  DECLASSIFICATION/  DOWNGRADING  SCHEDULE 


lb  RESTRICTIVE  MARKINGS 

3.  DISTRIBUTION  /AVAILABILITY  OF  REPORT 

Approved  lor  public  relaese, 
distribution  is  unlimited 


S.  MONITORING  ORGANIZATION  REPORT  NUM 

A£0SR-TR*  y  0  09  4 


a.  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL 

DEPARTMENT  OF  PHYSICS 

_ [WISE _ _ _  °MIST 

ADDRESS  ( City ,  Stttt.  tnd  ZlPCodt) 

Sackville  Street,  Manchester  M60  1QD 


7b.  ADDRESS  (City,  Stttt.^tnd.ZJP^C^ 

Bolling  AFB ,  DC  ^20332-6448 ^ J 


i.  NAME  OF  FUNOlNG/ SPONSORING 
ORGANIZATION 

AFOSR 


Be  ADDRESS  (City.  Stttt,  tnd  ZlPCodt) 

tx\  L^vC 

Bolling  AFB,  DC  20332-6448 


8b.  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

I'.*"'*  AFOSR-88— 0121 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT 

ELEMENT  NO.  NO 

61102F  '.^>3bV° 


7YT0/A 


1 1  TITLE  (Includt  Security  Cltuificttion) 

POLARIZATION  RADAR  STUDIES  OF  PRECIPITATION:  IMPLEMENTATION  OF  THE 
TECHNIQUE  AND  DATA  INTERPRETATION 


12.  PERSONAL  AUTHORS)  TT  ,  _  _ _  , _ 


ILLINGWORTH,  ANTHONY  JOHN 


13*.  TYPE  OF  REPORT 

FINAL 


13b.  TIME  COVEREO  14.  DATE  OF  REPORT  (Year,  Month,  Dty)  IS.  PAGE  COUNT 

FROM  Feb  88  TOFeb  90  1990  July  25 


COSATI  COOES 


GROUP  |  SUB-GROUP 


M 


IB.  SUBJECT  TERMS  (Continue  on  reverse  if  ntttsstry  tnd  idtntlfy  by  block  number) 

^POLARIZATION  RADAR,  ICE,  HAIL,  RAINDROPS, 
BRIGHT  BAND,  LIGHTNING. 


ABSTRACT  [Continue  on  rtv*rs*  if  ntasssry  snd  identify  by  block  numbtr)  ~x7x  f  ... 

Polarization  radar  observations  provide  information  on  characteristics  of  precipitation 
particles  not  available  with  conventional  weather  radar.  Wo  report  observations  of  four, 
parameters  made  with  the  25m  Chilbolton  dish;  the  largest  steerable  meteorological  rada.:  in 
<I*lese  are  the  first  S-band  measurements  of  the  linear  depolarisation  ratio 
and  the  most  accurate  co-  copolar  correlations  yet  reported.  Ten  publications  describe 
the  work  m  more  detail. 

reP°ft  danonstrates  hew  the  new  parameters  can  be  used  to:  differentiate  ice  fron 
water,  differentiate  the  different  forms  of  ice  (snow,  hail  pellets),  locate  areas  where 
large  hail  is  forming,  and  identify  clouds  posing  a  threat  of  triggered  lightning  before 
natural  lightning  or  breakdown  has  occurred.  lj,f  '  c 


zo  distribution  i  availability  of  abstract 
□  uNCLASSlflEPAjNLlMlTEO  □  SAME  AS  RPT 


lit  NAME  OF  RESPONSIBLE  INOIVIOUAL 


nn  i  »-*•? 


21.  ABSTRACT  SECURITY  CLASSIFICATION 
□  OTIC  USERS  UNCLASSIFIED 


22b.  TELEPHONE  t Includt  Arte  Codt) 

2021  767-4960 _  .. 


rhe  University  of  Manchester  Institute  of  ScienceandTechnology 

PO  Pox  SO.  Manchester  M60  1  QD.  Telephone  061-236  3311.  Teicx  666094 


I. 

j 

i 

i 

l 


afosr-tr- 


0  9  4  1 


9  0 


POLARISATION  RADAR  STUDIES  OF  PRECIPITATION: 
IMPLEMENTATION  OF  THE  TECHNIQUE  AND 
DATA  INTERPRETATION 

Contract  No  AFOSR  88-0121 

FINAL  SCIENTIFIC  REPORT 
15  Feb  1988  -  14  Feb  1990 


Anthony  J  Illingworth 
Dpartment  of  Pure  and  Applied  Physics 


UMIST 

P  0  Box  88 

Manchester  M60  1QD,  UK 

- 

r  < 

T! 

n 

_ _ 

July  1990 

■  .  -  t 

■  •  • '  - . . 

i/ 

:  ■  i'i 

/  Codes 

„  .  ^ iri 

-] 1;:"”  *  color 

'  repy-oduoW  . 

'*•-  -'1  iq  fcLaok  om3  ~l 

,  °  i 

ri  j 

f 

c.r.i/ or 
ial 

i 

Prepared  for: 

EUROPEAN  OFFICE  OF  AEROSPACE  RESEARCH  AND  DEVELOPMENT 

LONDON,  ENGLAND 


AppPot%ifrrru Mi-;  rolea-sj 
dt^trllwtionuulJLsiWd.  y. 


1. 


INTRODUCTION 


This  final  report  describes  work  performed  with  the 
Chilbolton  dual  polarisation  radar  under  contract  number  AFOSR- 
88-0121  during  1988  and  1989.  This  work  constitutes  an  extension 
to  the  radar  work  completed  on  grant  AFOSR-86-0193  covering  the 
years  1986  and  1987  (see  final  report  dated  20  Jan  1988). 

The  Chilbolton  dual  polarization  radar  is  the  largest  steer¬ 
able  meteorological  radar  in  the  world.  It  can  transmit  (and 
receive)  pulses  alternately  polarised  in  the  horizontal  and  ver¬ 
tical  directions.  The  radar  is  fully  operational  and  can 
currently  record  four  independent  polarization  parameters  con¬ 
taining  information  which  enables,  for  the  first  time,  the 
precipitation  particles  to  be  identified.  The  radar  operates  at 
S-band  (10cm)  so  that  propagation  effects,  which  can  be  trouble¬ 
some  at  C-band  and  X-band,  are  virtually  absent. 

We  now  summarise  the  polarization  parameters: 

(a)  Conventional  Radar  Reflectivity,  Z. 

Transmit  and  receive  horizontal  polarisation  and  measure  ZH, 
the  reflectivity  expressed  in  dBZ  relative  to  the  return  power 
for  one  mm  raindrop  per  cubic  meter.  From  Z  alone  we  cannot 
distinguish  ice  from  rain. 

(b)  Differential  Reflectivity,  ZDR 

The  ratio  of  ZH  to  ZV,  the  reflectivities  measured  at  hori¬ 
zontal  and  vertical  polarisations.  ZDR  is  a  measure  of  particle 
shape. 

(c)  The  linear  depolarisation  ratio,  LDR. 

Transmit  vertical  and  receive  the  horizontally  polarised 
component,  that  is,  the  depolarised  signal  ZVH,  and  compare  with 
ZH.  ZVH  is  finite  only  for  oblate  tumbling  particles.  LDR  senses 
fall  mode  and  is  an  excellent  detector  of  wet  ice. 

(d)  The  copolar  correlation  p(H,V). 

The  correlation  between  the  time  series  of  ZH  and  ZV  for 
successive  radar  pulses,  this  parameter  measures,  amongst  other 
things,  the  variety  of  shapes  present. 

The  Chilbolton  radar  has  made  the  first  S-band  observations 
of  LDR  and  the  most  accurate  measurements  of  p(H,V)  yet  reported. 
This  report  shows  how  these  new  parameters  may  used: 

(i)  To  distinguish  between  different  types  of  ice,  and  spec¬ 
ifically,  to  differentiate  between  areas  where  snow  is  forming 
and  those  where  soft  hail  (graupel)  is  present. 

(ii)  To  identify  and  correct  for  errors  in  rainfall  estimates 
caused  by  the  enhanced  return  from  melting  snow  (the  bright 
band ) . 

(iii)  To  identify  where  damaging  hail  may  be  forming. 

(iv)  To  locate  those  clouds  posing  a  threat  of  lightning,  even 
before  any  natural  lightning  has  occurred. 

The  colour  plate  overleaf  shows  how  the  radar  can  be  used  to 
detect  clouds  which  are  likely  to  produce  a  flash  of  lightning, 
even  though  no  natural  lightning  or  breakdown  may  have  occurred. 
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FRONTISPIECE 


RADAR  OBSERVATIONS  OF  A  TRIGGERED  LIGHTNING  EVENT 


The  upper  half  of  the  plot  is  a  vertical  section  through  a 
cloud  of  the  conventional  radar  reflectivity,  Z.  The  track  of  the 
penetrating  aircraft  at  a  height  of  2.3km  is  also  shown. 
Although  values  of  Z  at  this  height  were  only  modest  (40dBZ)  the 
aircraft  triggered  a  lightning  strike  as  it  penetrated  the  cloud. 
On  this  day  no  natural  lightning  had  been  observed  before  this 
time . 


The  lower  half  of  the  plot  displays  one  of  the  new  radar 
parameters,  LDR,  the  linear  depolarisation  ratio,  which  we 
believe  can  identify  those  parts  of  the  cloud  posing  a  risk  of 
triggered  lightning. 

The  different  types  of  ice  in  the  cloud  can  be  identified  by 
the  values  of  LDR  when  they  melt.  Values  of  LDR  near  -15dB  (red 
in  the  figure)  are  indicative  of  oblate  melting  snow  flakes. 
Such  regions  of  the  cloud  should  not  pose  a  lightning  risk.  In 
contrast  to  this,  a  small  region  from  65  to  73km  range  has  much 
lower  values  of  LDR,  near  -26dB  (green  in  the  figure),  where 
there  are  more  spherical  melting  ice  particles  in  the  form  of 
hail  pellets  (graupel).  The  2-D  probe  on  the  penetrating  aircraft 
confirmed  the  existence  of  snow  and  hail  pellets  in  these  two 
regions . 

The  triggered  lighting  event  occurred  when  the  aircraft 
traversed  a  region  where  the  radar  inferred  the  presence  of  small 
hail  pellets.  We  believe  that  the  LDR  parameter  has  the  poten¬ 
tial  to  identify  such  a  region  of  potential  lighting  risk. 
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2.  OBSERVATION  PROGRAM  in  1988  and  1999 

During  1988  and  1989  the  Chilbolton  radar  reported  the  first 
LDR  measurements  at  S-band,  and  also  made  measurements  of  p(H,V) 
of  unrivalled  accuracy. 


(a)  Measurements  have  been  made  on  numerous  days  at  Chilbolton 
in  stratiform  and  convective  precipitation.  The  first  coincident 
sets  of  data  for  the  four  parameters  (Z,  ZDR,  LDR  and  p(H,V)) 
have  been  collected  for  several  different  rainfall  types. 

(b)  On  three  days  coincident  fights  with  the  Met  Office  C-130 
aircraft  have  been  carried  out.  Excellent  data  sets  of  2-D 
probe  observations  of  precipitation  particle  type  with  simultan¬ 
eous  radar  coverage  were  obtained  in  ice  and  in  the  melting 
layer.  Aircraft  co-ordination  worked  well,  with  the  quarter 
degree  beam  radar  hitting  the  aircraft  during  most  cloud  pene¬ 
trations;  this  represents  a  navigation  accuracy  of  about  300m. 

(c)  On  13  July  1988  the  aircraft  triggered  a  lightning  flash 
while  passing  through  a  cloud  which  was  being  scanned  by  the 
radar. 

(d)  From  November  1988  to  January  1989  the  group  from  ETH 
(Zurich,  Switzerland)  located  their  mobile  vertically  pointing 
Doppler  radar  10km  upwind  of  Chilbolton  for  a  co-ordinated  study 
of  the  bright  band.  Vertical  profiles  with  50m  resolution  of  the 
velocity  and  polarization  parameters  were  obtained  through  bright 
bands  in  warm  fronts,  cold  fronts  and  warm  sector  rainbands. 


3.  PUBLICATIONS 

In  our  final  report  on  contract  no  AFOSR  86-0192  we  referred 
to  three  completed  publications: 

(i)  'Detection  of  hail  by  dual  polarisation  radar' 

Nature  320,  431-433,  1986. 

(ii)  'Polarization  radar  studies  of  precipitation  develop¬ 

ment  in  convective  storms' 

Q  J  Roy  Met  Soc,  113,  469-489,  1987. 

(iii)  'Radar  observations  and  modelling  of  warm  rain  initia¬ 

tion  ' 

0  J  Roy  Met  Soc  113,  1171-1191,  1987. 
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A  further  ten  publications  have  been  completed  in  the  two 
years  of  the  current  grant.  Copies  of  these  ten  publications  are 
included  in  this  report  as  they  best  summarise  the  work 
completed.  A  large  amount  of  data  is  still  being  analysed. 


1.  A  J  Illingworth  and  I  J  Cay lor 

'Radar  observations  and  modelling  of  warm  rain  initiation' 

1988  10th  Int  Conf  on  Cloud  Physics,  Bad  Homburg,  Germany. 

2.  A  J  Illingworth 

'The  formation  of  rain  in  convective  clouds' 

1988  Nature,  336,  No  6201,  754-756. 

3.  A  J  Illingworth  and  I  J  Cay lor 

'Polarization  radar  estimates  of  raindrop  size  spectra  and 
rainfall  rates' 

1989  J  Oceanic  and  Atmos  Tech,  6,  939-949. 

4.  A  J  Illingworth  and  I  J  Cay lor 

'Identification  of  precipitation  particles  using  dual  polar¬ 
ization  radar' 

1988  10th  Int  Conf  on  Cloud  Physics,  Bad  Homburg,  Germany. 

5.  I  J  Caylor  and  A  J  Illingworth 

'Dual  linear  polarisation  time  series  as  an  aid  to  hydro¬ 
meteor  identification' 

1988  I GARS S  '88  Symp.  Edinburgh,  Sept.  REF  ESA  SP-284. 

6.  I  J  Caylor  and  A  J  Illingworth 

’Identification  of  the  bright  band  and  hydrometeors  using 
co-polar  dual  polarization  radar’ 

1989  24th  Conf  on  Radar  Meteorology,  Tallahassee,  Florida. 

7.  A  J  Illingworth  and  I  J  Caylor 

'Cross  polar  observations  of  the  bright  band'. 

1989  24th  Conf  on  Radar  Meteorology,  Tallahassee,  Florida. 

8.  S  E  Hopper,  A  J  Illingworth  and  I  J  Caylor 

'Bright  band  errors  in  rainfall  measurement:  identification 
and  correction  using  linearly  polarised  radar  returns' 

Int  Symp  on  Hydrol  Appl  of  Weather  Radar,  Salford,  UK  1989. 

9.  I  J  Caylor,  J  W  F  Goddard,  S  E  Hopper  and  A  J  Illingworth 
'Bright  band  errors  in  radar  estimates  of  rainfall:  identi¬ 
fication  and  correction  using  polarization  diversity' 

1989  COST-73  Int  Semin  on  Weather  Radar  Networking  ,  Brussels. 

10.  I  Frost,  A  J  Illingworth  and  I  J  Caylor 

'Aircraft  and  polarization  radar  measurements  of  a  triggered 
lightning  event' 

1989  Int  Conf  on  Lightning  and  Static  Electricity,  Bath,  UK. 
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The  topics  covered  by  these  papers  can  be  considered  in 
seven  groups.  We  list  these  groups  below,  and  then  in  Section  4 
consider  the  findings  in  more  detail: 

(a)  DEVELOPMENT  OF  WARM  RAIN  (Papers  1  and  2) 

Polarisation  radar  observations  show  that  rain  (in  the 

absence  of  ice)  first  develops  as  a  very  low  concentration  (less 
than  one  per  cubic  meter)  of  large  drops  (above  4mm  diameter). 
Rain  normally  has  several  thousand  drops  per  cubic  meter. 
Modelling  of  USA  and  UK  data  suggests  that  these  first  raindrops 
are  formed  on  giant  nuclei.  This  is  of  relevance  for  theories  of 
rainfall  development  and  implications  for  attempted  weather 
modification.  Such  unexpectedly  large  raindrops  in  weak  echoes 
could  pose  an  erosion  threat  to  returning  space  vehicles. 

(b)  RAINDROP  SIZE  DISTRIBUTIONS  (Papers  3) 

This  paper  presents  new  data  on  the  shape  of  large  rain¬ 
drops,  and  uses  these  to  derive  raindrop  size  and  concentrations 
from  observations  of  Z  and  ZDR.  A  statistical  study  of  the 
variability  of  naturally  occurring  raindrop  spectra  is  also 
given. 

(C)  POLARISATION  PROPERTIES  OF  FLARE  ECHOES  (Papers  4  and  6) 

Flare  echoes  are  artefacts  which  result  from  triple  scatter¬ 
ing  via  the  ground  and  occur  in  low  reflectivity  regions  where  Z 
is  low.  Triple  scattering  only  occurs  in  the  horizontal  polaris¬ 
ation  and  so  flare  echoes  are  accompanied  by  very  high  positive 
values  of  ZDR. 

(d)  INTERPRETATION  OF  CO-POLAR  CORRELATION  (Papers  4,  5,  and  6) 
In  these  papers  we  report  the  most  accurate  measurements  yet 

obtained  of  the  co-polar  correlation.  They  show  that  it  may  be 
used  to  measure  the  breadth  of  the  raindrop  size  spectra,  to 
differentiate  between  snow  and  graupel ,  and  how  it  promises  to 
provide  a  means  of  identifying  where  hail  is  growing  within  a 
cloud . 

(e)  INTERPRETATION  OF  LINEAR  DEPOLARISATION  RATIO,  LDR  (Paper  7) 
This  paper  presents  the  first  measurements  of  LDR  at  S-band, 

and  shows  how  it  can  be  used  to  measure  the  canting  angle  of 
raindrops,  to  identify  the  bright  band,  and  to  distinguish 
between  snow  and  hail  pellets.  A  colour  plot  in  the  paper  shows 
the  very  different  values  of  LDR  for  the  different  types  of  ice. 

(f)  USE  OF  LDR  TO  CORRECT  FOR  BRIGHT  BAND  ERRORS  (Papers  8  and  9) 

The  bright  band  is  region  of  enhanced  radar  reflectivity 
caused  by  melting  snowflakes.  Values  of  Z  are  about  lOdB  higher 
than  in  the  rain,  and  these  high  values  of  Z  can  lead  to  overes¬ 
timates  in  inferred  rainfall  rates  by  up  to  a  factor  of  ten. 
These  two  papers  show  how  the  new  LDR  parameter  can  be  used  to 
identify  these  errors  and  also  suggests  methods  of  correcting 
them. 

(g)  RADAR  OBSERVATIONS  OF  TRIGGERED  LIGHTNING  (Paper  10) 

The  findings  of  this  paper  are  summarised  in  the  frontis¬ 
piece,  the  triggered  lightning  event  confirms  that  the  radar  is 
able  to  identify  those  regions  of  the  cloud  posing  a  lightning 
r  isk . 
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4.  DETAILED  ANALYSIS  OF  THE  SIX  TOPICS 


(A)  WARM  RAIN  DEVELOPMENT  (Papers  1  and  2). 

The  precise  mechanism  by  which  rain  is  formed  in  clouds  in 
which  ice  is  not  present  has  been  the  cause  of  much  discussion. 
The  cloud  drops  formed  by  condensation  are  usually  so  small  that 
they  do  not  possess  enough  momentum  to  collide  with  one  another 
to  grow  to  raindrops.  In  our  earlier  report  (AFOSR  86-0193)  we 
discussed  some  UK  measurements  which  indicated  that  early  echoes 
of  warm  convective  clouds,  although  having  quite  modest  values  of 
Z  of  only  20dBZ,  already  had  large  positive  values  of  ZDR;  these 
findings  indicate  the  presence  of  raindrops  of  3  or  4mm  dia¬ 
meter,  but  because  of  the  low  value  of  the  total  reflectivity, 
these  first  raindrops  must  be  present  in  concentrations  of  much 
less  than  one  for  every  ten  cubic  meters.  In  mature  clouds  rain¬ 
drops  are  present  in  concentrations  of  thousands  per  cubic  meter. 

In  these  papers  we  present  data  on  the  evolution  of  raindrop 
spectra.  We  show  that  such  echoes  are  common  in  the  USA  when  warm 
rain  processes  are  operating.  The  USA  data  come  from  the  CP2 
radar  operating  in  Alabama.  It  appears  that  initially  such  a  low 
concentration  of  large  drops  is  stable  with  very  little 
co 11 is ion- induced  drop  rupture;  when  a  critical  size  is  reached, 
the  drops  do  break  up,  the  concentration  increases,  more 
collisions  occur  and  the  distributions  rapidly  moves  towards  that 
normally  observed  in  mature  clouds. 

Such  warm  rain  processes  should  be  common  in  Florida.  The 
existence  of  such  large  drops  early  in  the  life  of  a  cloud  is 
quite  unexpected.  Aircraft  sampling  these  clouds  with  conven¬ 
tional  probes  would  not  sense  then,  but  they  could  cause  erosion 
to  returning  space  vehicles. 


(B)  RAINDROP  SIZE  DISTRIBUTIONS  (Paper  3) 

The  differential  reflectivity  (ZDR)  measures  the  mean  shape 
f f  hydrometeors  and  provides  an  estimate  of  the  mean  size  of 
raindrops.  Observations  of  ZDR  for  rain  may  be  combined  with  the 
conventional  radar  reflectivity  factor  (Z)  and  fitted  to  any  two- 
parameter  raindrop  size  distribution  and  this  information  used 
to  derive  more  accurate  rainfall  rates.  In  such  work  the  precise 
shape  of  raindrops  is  a  critical  parameter.  Recently  available 
data  suggest  that  large  raindrops  are  more  oblate  than  previously 
believed.  These  new  shapes  support  the  idea  that  ZDR  values  above 
3.5dB  can  be  attributed  to  rain.  Average  values  of  ZDR  as  a  func¬ 
tion  of  Z  obtained  in  heavy  rain  by  the  Chilbolton  radar  agree 
very  closely  with  those  predicted  using  the  new  shapes.  Statis¬ 
tics  are  presented  of  the  natural  variability  of  raindrop  spectra 
in  heavy  rain.  Analytic  expressions  are  proposed  for  computing 
rainfall  rates  from  Z  and  ZDR. 
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(C)  POLARIZATION  PROPERTIES  OF  FLARE  ECHOES  (Papers  4  and  6) 

Flare  echoes  are  caused  by  triple  scattering  from  a  high  Z 
region  at  a  height  h,  down  to  the  ground,  then  back  to  the  radar 
via  the  precipitation,  thus  forming  a  spurious  echo  a  distance  h 
being  the  intense  core.  It  has  already  been  established  that  the 
Doppler  velocities  in  flare  echoes  are  in  error.  We  report  that 
flare  echoes  are  invariably  accompanied  by  very  large  values  of 
positive  ZDR .  These  ZDR  values  do  not  reflect  the  properties  of 
the  hydrometeors,  but  are  an  artefact  of  the  triple  scattering. 
They  arise  because  induced  vertical  dipoles  do  not  radiate  along 
their  axes  down  to  the  ground,  and  ,  as  a  consequence,  all  the 
triple  scattered  power  is  in  the  horizontal,  and  large  positive 
ZDR  results. 

Analysis  of  the  triple  scattering  echoes  at  X  and  S-band 
confirms  that  for  a  given  value  of  Z  the  flare  echo  is  propor¬ 
tional  to  (wavelength) "4 ,  and  is  19dB  more  intense  at  X-band. 
For  some  large  storms  there  is  so  much  attenuation  at  X-band  that 
the  Z  values  on  the  far  side  are  so  reduced  that  the  flare  can  be 
absent.  Observations  of  the  Doppler  velocities  of  flares  using 
ZH  and  ZV  reveal  that  only  the  velocity  derived  from  ZH  are  in 
error.  Doppler  measurements  are  usually  made  using  the  ZH  channel 
but  our  findings  suggest  that  Doppler  problems  would  be  reduced 
if  the  ZV  channel  was  used. 


(D)  INTERPRETATION  OF  THE  CO-POLAR  CORRELATION  (Papers  4,5  and  6) 

This  new  variable  is  defined  as  the  correlation  between  the 
time  series  of  successive  estimates  of  ZH  and  ZV.  For  the 
Chilbolton  radar  64  pulse  pairs  are  transmitted  alternately 
polarised  in  the  horizontal  and  the  vertical i  To  compute  ZDR  the 
mean  values  of  these  ZH  and  ZV  estimates  are  used.  We  have  con¬ 
structed  a  system  for  recording  these  individual  estimates  at  64 
different  gates.  This  system  streams  on  to  magnetic  tape  at 
40kBytes  per  second. 

As  precipitation  particles  reshuffle  in  space,  then  the 
received  signal  fluctuates.  If  the  particles  are  all  the  same 
shape  then  both  ZH  and  ZV  fluctuate  in  the  same  way  and  the 
correlation  is  unity.  If  there  is  a  variety  of  different  shapes 
present,  then  the  correlation  falls. 

It  appears  that  Chilbolton  radar  is  able  to  make  these 
correlation  measurements  with  unrivalled  accuracy.  Values  in 
light  rain  are  over  0.995,  and  fall  to  0.98  in  heavier  rain  when 
larger  oblate  raindrops  coexist  with  the  smaller  spherical  ones. 
There  are  only  a  very  few  reports  in  the  literature  of  values  of 
P(H,V);  observations  with  CP2  and  at  NSSL  (both  in  the  USA)  seem 
to  reach  a  maximum  value  of  0.95.  The  reason  for  the  higher 
values  possible  at  Chilbolton  may  be  as  a  result  of  the  extreme 
purity  of  the  antenna  beam  pattern,  and  the  low  level  of  any 
ground  clutter.  Such  high  correlations  do  agree  with  theoretical 
values  computations  for  rain,  and  because  of  the  accuracy  avail¬ 
able  it  seems  possible  to  estimate  the  breadth  of  the  raindrop 
size  distribution.  If  the  raindrop  size  distribution  is  represen¬ 
ted  as  a  gamma  function  with  an  index  m,  then  the  value  of  m 
reflects  the  breadth  of  the  size  spectra  and  can  be  estimated 
from  p( H, V ) . 
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Very  low  values  of  p(H,V)  for  precipitation  are  confined  to 
the  bright  band  where  snow  is  melting,  and  large  oblate  parti¬ 
cles  flutter  and  rock  from  side  to  side  as  they  fall.  Indeed 

the  low  values  of  p(H,V)  can  be  used  to  identify  the  bright  band. 

Our  computations  show  that  extremely  low  values  of  p(H,V) 
should  occur  when  particles  become  large  enough  that  Mie  scatter¬ 
ing  occurs.  This  could  be  the  basis  of  detecting  hail  large 
enough  to  cause  damage.  Other  methods  of  hail  detection  suffer 
from  ambiguities  due  to  wet  or  dry  growth,  but  this  method  relies 
only  on  Mie  scattering  and  so  triggers  at  a  certain  hail  size  for 
a  given  wavelength. 

Ground  clutter  also  has  low  correlations.  We  expect  returns 
due  to  anomalous  propagation  also  to  have  low  values,  and  suggest 
that  correlation  values  could  be  used  to  detect  and  reject  spur¬ 
ious  echoes  due  to  anaprop. 

Some  advantages  of  the  correlation  technique  are: 

(i)  It  is  unaffected  by  differential  attenuation  and  so,  in  con¬ 
trast  to  ZDR,  should  be  applicable  at  shorter  wavelengths. 

(ii)  p(H, V)  can  be  estimated  using  a  short  time  series,  thus 
overcoming  the  long  dwell  times  sometimes  needed  to  measure  ZDR. 

(iii)  The  correlation  is  unaffected  by  differential  phase  shifts. 


(E)  INTERPRETATION  OF  THE  LINEAR  DEPOLARISATION  RATIO,  LDR 

(Paper  7) 

These  are  the  first  measurements  of  LDR  made  at  S-band.  At 
this  wavelength  the  values  of  LDR  are  unaffected  by  depolarisa¬ 
tion  of  the  incident  beam,  an  effect  which  makes  interpretation 
of  LDR  difficult  at  5cm  and  almost  impossible  at  3cm. 

The  colour  figure  in  this  paper  shows  LDR  measurements  made 
in  two  showers,  one  containing  graupel  and  a  second  with  snow. 
From  this  and  many  other  observations  we  conclude: 

(i)  LDR  values  in  the  rain  are  very  low,  and  are  usually  res¬ 
tricted  by  the  antenna  isolation  to  -32dB.  In  heavier  rain  LDR 
values  are  slightly  higher  and  are  consistent  with  a  mean  canting 
angle  of  raindrops  of  about  5  degrees. 

(ii)  LDR  values  in  dry  ice  are  usually  near  to  the  antenna 
limit,  this  is  because  at  microwave  frequencies  the  ice  parti¬ 
cles  generally  appear  spherical.  However,  in  stratiform  rain  some 
high  ZDR  regions,  usually  accompanied  by  moderate  Z,  can  have 
detectable  values  of  LDR.  Such  regions  are  thought  to  have  single 
crystals  of  ice,  and  because  of  their  higher  dielectric  constant 
and  oblate  shape  they  can  give  a  small  but  finite  cross  polar 
return.  Similarly,  areas  of  graupel  in  convective  clouds  can  have 
a  value  of  LDR  slightly  above  -32dB,  again  this  results  from  the 
more  dense  but  slightly  oblate  ice. 

(iii)  By  far  the  highest  values  of  LDR  are  associated  with  melt¬ 
ing  wet  ice.  Melting  snow  is  associated  with  values  which  are 
higher  than  -20dB,  and  can  be  interpreted  as  being  due  to  wet 
particles  with  an  axial  ratio  of  about  0.55.  Melting  graupel 
gives  LDR  values  of  about  -26dB,  consistent  with  an  axial  ratio 
of  about  0.85.  Aircraft  flights  have  confirmed  these  particle 
inferences.  We  suggest  that  LDR  is  an  excellent  method  of  differ¬ 
entiating  between  snow  and  graupel. 
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(iv)  Ground  clutter  has  values  of  LDR  above  -3dB.  Thus  LDR  has  a 
potential  use  in  identifying  and  rejecting  echoes  due  to  anaprop. 

Technologically  the  advantages  of  LDR  is  that  no  fast  switch 
is  required,  instead  two  receive  channels  are  needed.  One  prob¬ 
lem  is  that  the  low  power  in  the  cross-polar  channel  may  be 
susceptible  to  contamination  by  ground  clutter  if  the  antenna  has 
poor  sidelobe  performance.  Clearly  the  sensitivity  of  LDR  is 
also  limited  by  the  cross-polar  performance  of  the  antenna. 


(F)  USE  OF  LDR  TO  CORRECT  FOR  BRIGHT  BAND  ERRORS  (Papers  8  and  9) 

Papers  8  and  9  analyse  the  increased  values  of  reflectivity 
in  the  bright  band  which  result  from  melting  snow.  Rainfall 
values  derived  from  Z  can  as  a  consequence  lead  to  overestimates 
of  rain  rate  by  up  to  a  factor  of  ten.  These  papers  show  that 
values  of  LDR  above  -20dB  are  an  excellent  detector  of  the  bright 
band.  A  simple  algorithm  can  then  be  used  to  correct  the  value 
of  Z  and  improve  the  rain  rate  estimates.  A  statistical  study  of 
the  Z  and  LDR  data  suggests  that  rainfall  estimates  are  improved 
if  8dB  is  subtracted  from  the  measured  Z  when  the  accompanying 
value  of  LDR  is  above  -20dB. 


(G)  POLARIZATION  RADAR  MEASUREMENTS  OF  A  TRIGGERED  LIGHTING  FLASH. 

(Paper  9) 

The  passage  of  an  aircraft  through  clouds  can  trigger  light¬ 
ing  which  would  not  occur  naturally.  On  13  July  1988  a  triggered 
event  occurred  in  a  cloud  which  was  simultaneously  being  scanned 
by  the  Chilbolton  radar.  The  triggered  lighting  occurred  when 
the  aircraft  was  penetrating  an  area  of  the  cloud  where  LDR 
indicated  the  presence  of  graupel .  Direct  measurements  made  by 
the  aircraft  confirm  the  presence  of  graupel  in  this  region. 
The  frontispiece  to  this  report  displays  the  precise  values  of  Z 
and  LDR  when  the  triggered  event  occurred.  It  is  interesting  to 
note  that  this  cloud  has  only  modest  values  of  Z  and  a  low  echo 
top  and  is  not  an  obvious  candidate  for  electrical  activity. 

This  finding  has  important  practical  implications.  It 
suggests  that,  for  the  first  time,  it  is  possible  to  identify  a 
cloud  which  is  likely  to  trigger  lighting  before  any  natural 
lightning  or  breakdown  as  occurred.  We  are  currently  engaged  in 
further  work  for  EOARD  to  examine  this  hypothesis  further.  Dur¬ 
ing  the  summer  of  1989  we  are  comparing  lightning  location  data 
with  polarisation  radar  observations  to  check  that  lightning  is 
only  observed  in  clouds  where  the  radar  indicates  the  presence  of 
graupel . 
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l  .  INTRODUCTION 

Differential  radar  reflectivity  observations 
of  early  echoes  of  rain  in  warm  convective 
clouds  in  the  UK  and  Alabama,  USA,  indicate 
that  precipitat ioo  initially  forms  as  rain¬ 
drops  greater  than  4mm  in  size  but  present  in 
concentrations  of  less  than  one  per  cubic 
meter.  The  differencial  reflectivity,  ZDR,  (10 
log  (ZH/ZV)),  where  ZH  and  ZV  are  the  radar 
reflectivities  measured  with  horizontal  and 
vertical  polarizations,  respectively.  ZDR  pro¬ 
vides  a  measure  of  the  shape  of  the  hydro¬ 
meteors,  and  because  raindrops  are  oblate  to  a 
degree  which  depends  upon  their  size,  the 
magnitude  of  ZDR  is  a  unique  function  of  rain¬ 
drop  size  (Table  1).  ZDR  is  a  ratio  and  so  is 
independent  of  concentrations,  but  once  the 
mean  size  is  known  then  the  conventional  re¬ 
flectivity,  ZH,  can  be  used  to  derive  an  esti¬ 
mate  of  raindrop  concentration  (SELIGA  and 
BRINGI,  1976). 

2.  RAINDROP  SIZE  AND  CONCENTRATIONS 

If  we  assume  an  exponential  raindrop  size 

dist  r i but  ion : 

N(D)  3  No  exp(  -3.67  D/Do) 
where  N  is  the  concentration  of  drops  of 
diameter  D  and  Do  is  the  eq i vo 1 ume : r ic  dia¬ 
meter,  then  by  sumning  the  contributions  of 
the  various  sizes  of  raindrops  present  we  may 
calculate  the  values  of  ZDR  as  a  function  of 
Do  (Table  1,  final  columns,  spectrum  trun¬ 
cated  at  8mm).  The  value  of  No  is  then  lin¬ 
early  dependent  upon  the  magnitude  of  ZH.  In 
Figure  1  the  solid  lines  are  the  values  of  Z 
and  ZDR  as  Do  varies  but  No  is  kept  constant. 
The  total  drop  concen c r a t i on  is  given  by 
NoDo'3.67.  Many  observations  have  shown  the 
average  raindrop  size  distribution  to  be  that 
pr^pos^d  by  Marshall  and  Palmer  with  No=8000 
ro  mm  *  .  A  long  series  of  radar  observations 
of  Z  and  ZDR  (CAYLOR  and  ILLINCWORTH ,  1967) 

has  confirmed  that  the  average  value  of  ZDR  of 
rain  for  a  given  Z  does  indeed  lie  upon  the 
No=8000  curve  in  Figure  1. 

We  shall  be  discussing  early  echoes  which  have 
values  of  ZDR  of  3  or  4dB,  and  would  be  expec¬ 
ted  to  have  Z  values  of  56  to  67dBZ  for 
Mar sha 1 1 -Palmer  rain.  Instead  the  values  are 
30dBZ  lower,  implying  values  of  No  reduced  by 
a  factor  of  1000.  For  example  from  Figure  1 
we  see  that  a  value  of  Z  of  13d B2  accompanied 
bv  a^ZDR,ot  3dB  implies  a  value  of  No  of  only 
0.8m  ram^  ,  or  a  drop  concentration  (Do*2.5trm) 
of  u .  -m 


3.  CHILBOLTON  RADAR  (UK)  OBSERVATIONS 
A  vertical  section  through  a  young  convective 
cloud  with  anomalously  high  ZDR  is  displayed 
in  Figure  2  (for  other  examples  and  details 
of  the  radar,  see  ILLINGWORTH  et  al  ,  1987  ). 

Values  of  Z  of  30dBZ  are  accompanied  by  a  ZDR 
above  4 . 5dB ,  which  imply  (Figure  1)  large 
drops  (6mm)  in  concen t  rat  ions  of  only  one  per 
cubic  meter.  This  cloud  persisted  for  twenty 
minutes  with  no  great  change  in  character;  the 
2700  data  points  obtained  during  this  period 
are  displayed  in  Figure  3,  where  each  star 
represents  the  average  value  of  ZDR  for  every 
2dBZ  step  in  Z.  The  solid  line  in  Figure  3  is 
the  No=8000  curve  from  Figure  1;  showing  that 
the  low  concentration  of  large  raindrops  per¬ 
sisted  throughout  the  twenty  minute  period. 
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U.  521 

6.70 

4 

4.22 

Values  of  axial  ratio  (from  3EARD  and  CHUANC. 
1987)  and  ZDR  for  various  sizes  of  raindrop. 
The  Mie-Gans  calculations  are  supplied  bv  Dr 
Holt,  Department  of  Mathematics,  University  of 
Essex  ,  and  apply  to  raindrops  at  Q°C  and 
3.0765  GHz. 
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Figure  2.  A  vertical  scan  on  20  June  1980 
displaying  anomalously  high  ZDR. 


Figure  3.  Averaged  2  and  ZDR  over  a  20  minute 
period  for  the  cell  in  Figure  2.  For  compar¬ 
ison  the  model  output  as  a  function  ot  time 
(minutes)  and  the  Z  and  ZDR  for  a  drop  size 
distribution  with  No  =  8000  are  shown. 

4.  GIANT  NUCLEI  MODEL 

To  explain  such  low  concentrations  of  large 
raindrops,  CAYLOR  and  ILLINGWORTH  (1987)  ex¬ 
tended  a  suggestion  (JOHNSON,  1982)  that  each 
large  raindrop  formed  on  an  ultra  giant  nuclei 
( >30um)  ,  present  in  the  background  concentra¬ 
tion  at  this  low  level  (JUNGE,  1972).  Figure  3 
shows  the  prediction  of  the  simple  model  in 
which  Z  and  ZDR  increase  as  the  giant  nuclei 
sweep  out  cloud  water  of  LWC  2  g  m-3 ,  the 
numbers  indicate  time  in  minutes  elapsed  since 
initiation  of  the  mode^  with  nuclei  concentra¬ 
tions  ( m  units. of  m  per  unit  log  interval 
of  radiu9)of  10  for  I0um  radius  to  3X10 
for  lOGum. 


Only  a  few  large  raindrops  grow,  and  the  de¬ 
pletion  of  the  cloud  water  is  negligible.  The 
trodel  is  linear  in  chat  increasing  the  LWC 
merely  changes  the  time  scale  but  not  the 
Z/ZDR  dependence;  in  addition,  increasing  the 
nuclei  concentration  merely  scales  up  the 
values  of  Z  leaving  ZDR  unchanged.  The 
Chilbolton  cloud  persisted  20  minutes  and  it 
is  not  possible  to  give  a  specific  age  to 
various  parts  of  the  cloud.  In  view  of  the 
uncertainty  of  the  absolute  concentration  of 
giant  nuclei,  the  observed  general  dependence 
of  Z  with  ZDR  is  consistent  with  the  model. 

5.  ALABAMA  DATA,  NCAR  CP 2  RADAR 

First  echoes  formed  at  temperatures  above  zero 
on  many  days  during  the  1986  MIST  (Microburst 
Severe  Thunderstorm)  project.  These  early 
echoes  were  also  normally  accompanied  by  an¬ 
omalously  high  values  of  ZDR.  A  particularly 
clear  example  was  on  10  July,  when  the  evolu¬ 
tion  of  a  very  weak  echo  less  than  2km  in 
diameter  was  observed  in  its  entirety.  In  20 
minutes  the  echo  grew  from  zero  to  20dBZ  and 
then  collapsed;  the  top  of  the  radar  echo 
reached  only  5km  (slightly  above  the  freezing 
level),  but  during  this  time  values  of  ZDR  in 
the  cloud  increased  to  4dB,  the  ZDR  core  then 
descended  as  the  rain  fell  to  ground. 

Figures  4  and  5  show  the  specific  values  of 
each  Z/ZDR  data  point  for  two  successive  ver¬ 
tical  sections  through  the  maximum  echo.  The 
numbers  plotted  for  each  point  represent  the 
height  of  the  data  to  the  nearest  kilometer. 
The  evolution  is  consistent  with  the  model  in 
Figure  3,  in  that  the  concentration  of  drops 
is  constant  and  less  than  one  per  cubic 
meter,  5  minutes  after  the  appearance  of  the 
echo,  the  value  of  ZDR  was  3dB,  growing  to  4dB 
in  the  subsequent  three  minutes.  This  is 
consistent  with  the  maximum  drop  size  in¬ 
creasing  from  4.5  to  6mm  (Tab^e  l)  by  sweeping 
out  cloud  water  of  LWC  2  g  m 

6.  CONCLUSIONS. 

We  have  presented  evidence  that  early  echoes 
of  warm  clouds  consist  of  a  very  few  drops 
which  grow  to  a  large  size  by  sweeping  out 
cloud  water.  Because  concentrations  are  so  low 
this  can  happen  without  exhausting  the  supply 
of  liquid  water,  and  also  vich  a  negligible 
number  of  collisions  causing  the  drops  ro 
break  up.  LOW  and  LIST  (1982)  found  that  for 
a  collision  to  cause  shattering,  both  rain¬ 
drops  must  be  larger  than  Iram.  For  a  given  Do 
the  lifetime  of  a  large  raindrop  will  be  in¬ 
versely  proportional  to  No;  in  Figure  1  the 
dotted  lines  join  the  values  of  Z  and  ZDR  for 
lifetimes  of  1000,  10000,  and  100000  seconds 
for  a  5tnm  drop.  The  large  raindrops  in  Figure 
4  and  5  will  survive  for  over  an  hour  before 
collision-induced  rupture.  This  situation 
should  be  contrasted  to  Figure  6  showing  the 
Z/ZDR  scatter  plot  for  a  normal  mature  cloud 
on  6  July  1986,  the  values  are  closer  to  the 
average  No*8000  curve,  and  the  lifetime  of  the 
larger  drops  is  less  than  a  minute. 


Warm  rain  measurements  in  Texas  (CARBONE  and 
NELSON,  1978)  and  in  Hawaii  (BEARD  et  al  , 
1986)  also  indicate  low  values  of  No  and  high 
Do.  Conventional  airborne  instruments  do  not  _ 
generally  have  a  large  enough  sample  volume  to  ^ 
obtain  meaningful  statistics  of  the  very  low  — 
concentrations  inferred  from  the  radar.  In 
Hawaii  one  8mm  raindrop  was  sampled,  suppor¬ 
ting  the  idea  that  such  raindrops  are  indeed 
stable  in  the  atmosphere  if  they  do  not  under¬ 
go  collisons.  Drops  up  to  9mm  are  stable  in 
low  turbulence  wind  tunnels  (PRUPPACHER  and 
BEARD,  19  70). 

We  suggest  that  convective  raindrops  first 
form  on  embryos  which  are  present  in  concen¬ 
trations  of  about  one  per  cubic  meter.  Such 
embryos  could  be  giant  nuclei,  or  alternative¬ 
ly  could  be  ice  crystals.  The  number  of  ice 
crystals  seeding  such  a  cloud  might  be  expec¬ 
ted  to  be  more  variable,  and  to  be  much 
greater  if  the  cloud  top  was  higher,  vet  the 
raindrop  concent  rat  ions  always  seem  to  be  very 
low.  A  raindrop  distribution  of  a  few  large 
drops  would  appear  to  be  very  stable,  but  once 
break-up  occurs,  either  due  to  the  occasional 
collison  or  spontaneously  if  the  drops  become  ZZ 
too  large,  then  many  small  fragments  are  pro-  3 
duced  ,  and  the  distribution  rapidly  and  ir- 
reversibly  changes  to  the  more  normal 
Mar sha 1 l -Palmer  distribution. 
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Figure  4.  Scatter  plot  of  Z  and  ZDR  for  each 
data  point  from  a  young  echo  on  10  July  1986 
at  1301  CDT  indicating  low  concentrations  of 
large  drops . 
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Figure  5.  Values  of  Z  and  ZDR  from  the  same 
cloud  as  Figure  4  observed  3  minutes  later. 
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Figure  6.  Z  and  ZDR  data  points  from  a  mature 
convective  cloud  on  6  July  198b  at  1426  CDT 
indicating  a  Mar sha 1 1 -Palme r  drop  si2e 
distribution. 
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One  of  the  chief  problems  in  cloud  physics'  is  to  explain  bow,  in 
the  absence  of  ice,  small  cloud  droplets  can  coalesce  to  form 
raindrops.  Cloud  droplets  are  produced  in  concentrations  of 
— 10*  m-3  by  condensation  on  sub-micrometre  atmospheric  nuclei, 
but.  because  they  are  so  small,  collisions  between  droplets  are  rare 
events.  In  mature  rainciouds  there  are2  usually  about  1,000  drops 
per  m3;  here,  however,  I  present  differential  radar  reflectivity 
observations  which  indicate  that,  when  ice  is  not  present,  first 
echoes  of  convective  showers  consist  of  large  (>4  mm)  raindrops 
present  in  much  lower  concentrations  (<  1  m').  I  suggest  that 
these  raindrops  form  on  ultra-giant  nuclei  ( >  30  pm  radius)  which 
are  present  in  similarly  low  background  concentrations3.  In  some 
clouds  this  initial  raindrop  spectrum  persists,  whereas  in  other 
cases  it  changes  rapidly  and  irreversibly  to  the  raindrop  distribu¬ 
tion  normally  observed  in  mature  clouds  with  many  more  smaller 
raindrops.  The  latter  may  occur  when  the  break-up  of  the  large 
raindrops  becomes  more  common  and  many  small  satellite  drops 
are  produced. 

Coalescence  of  small  droplets  in  clouds  may  be  triggered  by 
condensation  on  hygroscopic  salt  nuclei  or  by  mixing  of  drier 
air  into  the  cloud;  either  process  could  produce  larger  droplets, 
which  would  then  have  a  reasonable  terminal  velocity  and  be 
able  to  capture  smaller  droplets,  thus  triggering  droplet  growth. 
Paluch  and  Knight4  have  reviewed  the  evidence  for  these  proces¬ 
ses  The  cloud  droplet  concentration  is  at  least  five  orders  of 
magnitude  higher  than  the  raindrop  concentration,  so  that  it  is 
very  difficult  to  make  direct  in-cloud  obsenations  of  the 
occasional  embryonic  raindrop  in  the  radius  range  50-100  pm, 
when  so  many  smaller  droplets  are  also  present.  Measurements 
of  the  conventional  radar  reflectivity,  Z.  are  difficult  to  interpret 
in  terms  of  both  the  size  ( d  i  and  concentration  (  n  i  of  raindrops, 
because  Z  is  given  by  the  product  nd'-  summed  over  all  raindrop 
sizes  Z  is  usually  quoted  in  dBZ,  that  is  in  decibels  over  t  mnT 
mm  A  full  review  of  the  meteorological  applications  of  radar 
is  provided  by  Browning". 

Here  we  report  measurements  of  differential  reflectivity  i  Zn„  t, 
which  is  defined  as 

Z^r  ~  10  l°g  i  Zh  Z\  i  ll* 

where  ZH  and  Z,  are  the  radar  reflectivity  factors  measured 
w  ith  horizontally  and  vertically  polarized  radiation,  respectively 
Z,  )R  is  a  measure  of  the  shape  of  precipitation  particles.  Because 
raindrops  are  oblate  to  a  degree  which  depends  on  their  size. 
Z,  ,K  is  positive  for  large  raindrops  and  its  magnitude  is  a  measure 
of  their  mean  size  With  knowledge  of  the  drop  size.  d.  the  drop 
concentration  n  may  be  calculated'1  from  the  absolute  value  of 
Z.  In  practice  the  raindrop  size  distribution  may  be  well 
approximated  by  a  spectrum  of  the  form 

mdi  =  n,  exp  i -3.6  ~d  d„  i  i2i 

where  d  is  the  equivolumetric  diameter  The  value  of  Z,,R  is 
independent  of  concentration  and  may  be  calculated  as  a  func¬ 
tion  of  d  The  value  of  n,,  may  then  be  found  from  the  observed 
magnitude  of  Z  The  solid  curve  in  Fig.  1  shows  calculated 
v  alues  of  Z  and  ZI)R  for  n„  =  8,000  m  1  mm  this  value  of  n, 
was  originally  proposed  for  stratiform  rain-  but  it  now  appears 
to  describe  the  average  properties  of  most  ram  For  individual 
data  points  and  clouds  not  lying  on  this  average  curve,  the  v  alue 
ol  n  mav  be  derived  by  noting  that  Z  is  proportional  to  n, 
The  total  concentration  is  given  by  n„d„  3.67 

Observations  of  developing  convection  were  made  each  dav 
during  the  MIST  proiect  in  1986  using  the  10-cm  dual  polariz- 
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Fig.  1  Values  of  Z  (radar  reflectivity  i  and  ZDR  I  differential  reflec¬ 
tivity  i  observed  during  the  evolution  of  the  echo  on  10  July  1986 
The  solid  line  represents  the  drop  concentrations  for  average  rain, 
n0  =  8,000 nT3  mm’1  (equation  3).  The  observations  are  average 
v  alues  of  ZDR  for  each  2-dBZ  step  in  Z.  Symbols  represent  observa¬ 
tion  times:  •  12:57,  ~  13:01,  *  13:04.  n„  scales  linearly  with  Z. 
so  inferred  raindrop  concentrations  are  over  30 dB  la  factor  of 
one  thousand  I  below  'average'  The  numbers  indicate  time  in 
minutes  elapsed  since  the  initial  evolution  of  the  radar  responses 
Z  and  ZDR,  during  which  time  ultra-giant  nuclei  sweep  out  cloud 
water  droplets 

ation  CP-2  radar  operated  by  the  National  Center  for  Atmos¬ 
pheric  Research.  The  vaLues  of  Z  and  ZDR  reported  here  are 
accurate  to  better  than  1  dBZ  and  0.2  dB  respectively.  The  first 
echoes  of  5- 10  dBZ  usually  formed  at  temperatures  warmer  than 
0°C  and  were  accompanied  by  anomalously  positive  values  of 
ZDR ,  implying  ( Fig.  1 1  the  presence  of  a  very  low  concentration 
of  unexpectedly  large  raindrops  and  confirming  more  limited 
UK  observations'’'*  made  in  1983  and  1984  with  the  Chilbolton 
polarization  radar 

A  particularly  well  observed  case  occurred  on  10  July  1986, 
when  a  series  of  90  vertical  (RHIl  and  horizontal  'PPIl  radar 
sections  was  made  through  an  isolated  echo  from  its  initiation 
with  Z  =  5  dBZ  at  12:54  c ST  i  =  GMT -6  hours i  to  a  maximum 
of  Z  =  20  dBZ  ( accompanied  by  values  of  ZDR  near  4  dB).  and 
terminating  at  13:  10  as  the  echo  fell  to  the  ground  and  disap¬ 
peared.  During  this  time  the  echo  was  unsheared,  its  top  never 
rose  above  5  km  altitude  and  the  horizontal  cross-section  never 
exceeded  2  km  in  diameter.  It  is  not  possible  to  plot  the  hundred 
or  so  data  points  obtained  for  each  scan,  but  a  summary  of  the 
evolution  is  provided  in  Fig.  1.  in  which  the  average  value  of 
ZOH  for  each  2-dBZ  step  in  Z  is  plotted  for  three  scans  at  12 :  5". 
13:01  and  13  :04  These  scans  are  all  at  the  same  azimuth  and 
pass  through  the  maximum  echo  We  note  that  the  maximum 
v  alues  of  Z  and  ZI>R  increase  w  ith  the  passage  of  time,  but  that 
the  average  values  all  lie  on  a  similar  curve  P'ots  for  the 
intermediate  scans  show  a  smooth  progression  along  this 
evolutionary  curve  Increases  in  Z  ceased  at  13:04,  as  the  main 
echo  fell  below  cloud  base  (inferred  from  a  nearby  sonde  ascent  i 
and  the  plots  for  the  scan  at  13  :0’  are  within  0.5  dB  of  the 
13  :04  data. 

The  evolution  shown  in  Fig.  I  is  compatible  with  a  situation 
in  which  a  very  few  raindrops  grow  by  sweeping  out  cloud  liquid 
water.  In  3  min  the  maximum  ZnR  increased  from  2.4  dB  to 
4dB.  corresponding  to  an  increase  of  the  maximum  diameter 
for  monodispersed  drops  from  -4  mm  to  6  mm.  this  occurring 
by  collection  of  cloud  droplets  with  a  liquid  water  content  ol 
2gm  '  The  concentrations  are  very  low ;  for  example,  a  Z  value 
of  15  dBZ  accompanied  by  a  ZnR  of  3dB  is  equivalent,  for 
monodispersed  drops,  to  the  presence  of  4  5-mm  drops  at  a 
concentration  of  0.03  m  .  or.  for  an  exponential  distribution, 
to  a  value  of  n„  of  0  5  m  mm  1  The  difference  in  concentra¬ 
tions  predicted  by  the  exponential  and  the  monodispersed 
spectra  is  negligible'  compared  with  the  deviation  from  the 
jverage  curve  in  Fig  1  Analysis  of  the  sonde  ascent  predicts 
a  maximum  adiuhatis  liquid  water  content  of  "'tm  for  a 


Fig.  2  The  evolution  of  the  echo  observed  on  28  June  1986.  Solid 
line  and  averages  follow  the  same  key  as  in  Fig  1.  Times  are: 

•  52.21,  ”  12:2'.  -  12:40.  The  dashed  line  joins  the  values  of  Z 
and  ZDR  that  correspond  to  a  5  mm  raindrop  which  has  a  lifetime 
of  1. 000s  before  shattering  as  a  result  oi  collisions  For  a  given 
ZDR  the  lifetime  is  inversely  proportional  to  Z 

temperature  of  8‘“C.  Above  this  level,  much  drier  air,  with  a 
mixing  ratio  of  less  than  2  g  m'J  at  5  °C,  restricted  cloud  growth. 

The  data  are  in  very  good  agreement  with  the  time  evolution 
predicted  by  a  model"  in  which  ultra-giant  nuclei  collect  cloud 
water  and  grow  to  raindrops  This  model  extends  earlier  work10 
which  showed  that  such  a  mechanism  could  produce  a  high 
radar  reflectivity  within  15-20  min  Our  new  model  predicts  the 
accompanying  values  of  ZDR.  but  does  not  consider  the  vertical 
structure  that  may  result  from  sedimentation.  The  nuclei  need 
not  be  hygroscopic  but,  because  they  are  so  large,  they  collect 
the  cloud  droplets  with  reasonable  efficiency.  The  numbers  in 
Fig  1  indicate  the  time  in  minutes  elapsed  since  initiation  of 
the  mechanism  proposed  in  our  model,  for  a  cloud  liquid  water 
density  of  2  g  m  "  with  concentrations  tin  units  of  m~!  per  unit 
log  interval  of  radius  l  of  1,000  for  nuclei  of  10-q.m  radius  to 
0.03  for  100-p.m  radius  These  nuclei  concentrations  are  those 
quoted  by  Junge'  for  the  background  level  found  above  the 
trade-wind  inversion  and  are  similar  to  the  maritime'  nucleus 
spectrum  used  by  Johnson111  The  drop  concentrations  inferred 
on  10  July  are  exceptionally  low,  values  a  factor  of  ten  higher 
were  inferred  on  other  days  in  the  MIST  programme.  It  may  be 
that  or  vuch  days  nucleus  concentrations  closer  to  the  ground 
were  higher  than  the  background  levels  used  in  the  model 
When  the  convection  was  vigorous,  the  intensifying  echo  top 
usually  rose  above  the  0  °C  isotherm  and  the  spectra  in  the  rain 
tended  towards  that  usually  observed  in  mature  clouds,  but  on 
one  occasion  1 28  June  19861  it  was  possible  to  observe  the 
evolution  of  an  echo  to  over  50dBZ  without  ice  forming  in  the 
cloud  The  average  values  of  ZOR  for  each  2-dBZ  step  in  Z  for 
three  successive  scans  are  shown  in  Fig.  2.  The  time  resolution 
is  rather  poor  because  the  radar  was  executing  360s  surveillance 
>vans.  At  12:21  the  maximum  Z  value  was  17dBZ,  associated 
with  a  ZnR  value  of  ~1.5  dB,  and  inferred  concentrations  are 
more  than  two  orders  of  magnitude  (equivalent  to  20  dB  I  below 
the  curve  for  average  rain.  Six  minutes  later,  at  12:27,  Z  had 
reached  nearly  40dBZ  and  ZnR  was  —3  dB.  indicating  that  the 
biggest  drops  had  grown  in  diameter  by  2  mm  The  concentration 
of  these  large  drops  is  still  two  orders  of  magnitude  below 
average  By  12  40  the  highest  Z  was  over  55  dBZ,  but  the  picture 
is  then  dramatically  different,  the  inferred  raindrop  size  distribu¬ 
tion  being  that  found  in  mature  clouds. 

To  support  our  contention  that  drop  collisions  are  responsible 
tor  this  abrupt  change  in  the  spectra,  the  dashed  line  in  Fig  2 
represents  values  of  Z  and  Zt)R  for  which  a  5-mm  droplet  would 
hav e  an  av  erage  lifetime  of  1 .000  s  before  it  collided  with  another 
droplet  greater  than  I  mm  in  diameter  C  ollisions  between  two 
raindrops  which  are  both  larger  than  I  mm  are  generallv  fol¬ 
lowed  bv  disruption  and  shattering  to  produce  mans  small 


fragments"  For  a  given  ZDR,  the  value  of  du  is  constant  and 
the  number  of  drops  is  proportional  to  Z;  accordingly,  the 
lifetime  of  each  large  raindrop  is  inversely  proportional  to  Z. 
We  note  that  in  Fig.  1  the  data  lie  tn  a  region  where  the  large 
drops  have  a  lifetime  of  several  hours  before  they  suffer  collision- 
induced  rupture.  In  Fig.  2  the  initial  concentrations  are  some¬ 
what  higher,  and  by  12:27  the  drop  lifetimes  are  slightly  shorter 
than  20  min;  finally,  at  12 :40h  the  spectrum  has  changed  to  the 
mature  mode,  with  a  high  concentration  of  smaller  drops,  so 
that  the  large  drops  contributing  to  a  ZDR  of  3-4  dB  would  have 
a  lifetime  of  only  a  few  seconds  and  would  be  continuously 
regenerated. 

Some  measurements 1:1 3  made  from  aircraft  within  developing 
clouds  have  also  indicated  raindrop  spectra  different  from  those 
found  in  mature  clouds,  although  Z  values  were  higher  and  the 
spectra  were  not  as  extreme  as  those  reported  here.  Measure¬ 
ments  in  Texas'3  revealed  low  values  of  n0  in  early  echoes  and 
an  implied  absence  of  collision-induced  break-up.  The  authors 
attributed  the  low  values  of  n0  to  drop  sorting  in  the  updraft, 
but  this  does  not  seem  compatible  with  the  time  evolution  of 
the  ZDR  RHI  scans  from  the  MIST  data.  In  Hawaii13  raindrops 
8  mm  in  size  were  observed.  The  existence  of  such  large  rain¬ 
drops  had  been  questioned  because  in  Marshall- Palmer  rain 
they  would  rapidly  suffer  collision-induced  break-up.  but  the 
authors  showed  that  when  total  drop  concentrations  are  low, 
the  lifetimes  are  sufficiently  long.  The  very  low  concentrations 
inferred  from  the  radar  data  of  the  MIST  project  would  be 
difficult  to  confirm  by  aircraft  penetrations  because  of  the  limited 
sample  volume  of  most  precipitation-sizing  probes. 

Other  mechanisms  for  the  production  of  the  raindrops  must 
be  considered  There  is  always  the  possibility  that  a  few  ice 
crystals  ( which  would  have  a  Z  value  below  -20  dBZ  and  would 
be  undetectable  by  radar)  could  fall  from  higher  levels,  melt, 
and  then  grow  into  large  raindrops.  This  seems  unlikely  :  the  air 
was  very  dry  at  these  colder  levels  and  no  bright  band  in  ZDR 
was  observed  (It  is  our  experience  that  as  small  ice  crystals 
start  to  melt  there  is  a  sudden  jump  in  ZDR  as  they  become  wet 
and  change  their  dielectric  constant,  followed  by  an  equally 
sudden  drop  to  zero  ZDR  as  they  melt  completely  to  become 
small  spherical  raindrops  ! 

In  his  initial  report  on  the  existence  of  ultra-giant  nuclei, 
Junge'  remarked  that  they  were  so  rare  that  they  could  not  be 
important  in  the  w  arm-rain  process.  Such  a  view  seems  eminently 
sensible,  but  this  new  radar  evidence  shows  that  warm  rain  does 
indeed  first  form  in  such  low  concentrations  It  may  be  that  a 
fortunate  combination  of  collection  efficiencies  and  droplet  sizes 
results  in  just  one  cloud  droplet  in  10"  being  able  to  capture  its 
neighbours  and  grow  into  a  raindrop.  It  is  evident  that  further 
modelling  and  observations  are  needed,  but  1  propose  mean¬ 
while  that  cloud  droplet  capture  by  ultra-giant  nuclei  is  a  simple 
and  equally  plausible  mechanism. 

I  thank  Professors  Bring!,  Forbes,  Fujita  and  Wakimoto  for 
allowing  access  to  the  MIST  data,  and  staff  at  NCAR  for  their 
assistance  during  my  visit  This  research  was  also  supported  by 
NERC  and  EOARD  NCAR  is  supported  by  the  National 
Science  Foundation 
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ABSTRACT 


The  differential  reflectivity  ( Zpp )  measures  the  mean  shape 
of  hydrometeors  and  provides  an  estimate  of  the  mean  size  of 
raindrops.  Observations  of  Zpp  for  rain  may  be  combined  with  the 
conventional  radar  reflectivity  factor  (Z)  and  fitted  to  any  two 
parameter  raindrop  size  distribution  and  this  information  used  to 
derive  more  accurate  rainfall  rates.  In  such  work  the  precise 
shape  of  raindrops  is  a  critical  parameter.  Recently  available 
data  suggest  that  large  raindrops  are  more  oblate  than  previously 
believed.  These  new  shapes  support  the  idea  that  Zqr  values 
above  3 . 5dB  can  be  attributed  to  rain.  Average  values  of  ZDp 
as  a  function  of  Z  obtained  in  heavy  rain  by  the  Chilbolton  radar 
agree  very  closely  with  those  predicted  using  the  new  shapes. 
Statistics  are  also  presented  of  the  natural  variability  of 
raindrop  spectra  in  heavy  rain.  Analytic  expressions  are 
proposed  for  computing  rainfall  rate  from  Z  and  Zqr. 


INTRODUCTION 


1  . 

The  radar  reflectivity  factor,  Z.  of  rain  measured  by  a 
conventional  radar  is  not  a  unique  function  of  the  rainfall  rate 
or  the  mean  size  of  the  raindrops,  but  is  given  by  the  product 
ND^,  where  N  is  the  concentrat ion  of  drops  of  diameter  D,  summed 
over  all  the  drop  sizes  present.  Many  different  empirical  rela¬ 
tionships  have  been  proposed  relating, the  rainfall  rate  (R)  to 
the  radar  reflectivity  (Z).  The  use  of  one  of  these  formulae  is 
essentially  equivalent  to  assuming  a  constant  raindrop  size 
distribution.  Such  relationships  reflect  an  average  long  term 
dependence,  but  the  actual  rainfall  predicted  by  an  individual 
value  of  Z  is  typically  in  error  by  a  factor  of  two  because  of 
natural  fluctuations  in  the  raindrop  size  spectra  (Wilson  and 
Brandes ,  1979  )  . 

Polarization  diversity  radar  provides  a  second  observable 
parameter,  the  differential  reflectivity,  which  is  related  to  the 
mean  drop  size,  and  when  this  information  is  used  in  conjunction 
with  Z,  more  accurate  rainfall  rates  may  be  derived.  The 
differential  reflectivity  (  ZqR  )  is  defined  as: 

ZDR  =  10  log  (  Zpj/  Z\j )  (1) 

where  Z R  and  Zy  are  the  radar  reflectivity  factors  measured  at 
horizontal  and  vertical  polarizations,  respectively.  For 
spherical  particles,  such  as  small  raindrops,  and  Zy  are 

equal,  and  ZDR  is  zero.  Larger  raindrops  are  oblate  to  a  degree 
which  depends  upon  their  size,  so  ZR  exceeds  Zy;  ZDR  is  positive, 
and  its  magnitude  is  related  to  the  mean  raindrop  size  present. 

When  observations  of  both  Z  and  ZqR  are  available,  then  the 
data  can  be  fitted  to  any  two  parameter  raindrop  site  distribu¬ 
tion.  Raindrop  size  distributions  are  commonly  of  the 
exponential  form: 


N  ( D ) 


NQ  exp  (-3.67  D/D0) 


(2  j 


where  D  is  the  equivolumetric  raindrop  diameter  of  the  oblate 

drops  and  D0  is  the  volume-median  drop  diameter.  Marshall  and 

Palmer  (1948)  reported  that  D0  is  a  function  of  rainfall  rate  and 
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that,  on  the  average,  N0  is  8000m  mm  .  Seliga  and  Bringi 
(1976)  proposed  that  the  value  of  D0  may  be  derived  from  the 
magnitude  of  Zqr,  and  N0  may  then  be  found  if  Z  is  known. 

The  differential  reflectivity  technique  relies  on  a  precise 
knowledge  of  the  shapes  of  different  sizes  of  raindrops.  The 
Chilbolton  radar  (Cherry  and  Goddard,  1982)  can  measure  Zgp  to  an 
accuracy  of  about  O.ldB  which  is  equivalent  to  a  change  in  rain¬ 
drop  axial  ratio  of  about  0.01.  Most  Zqr  measurements  have  been 
interpreted  using  the  laboratory  based  equilibrium  raindrop 
shapes  of  Pruppacher  and  Pitter  (1971)  who  quote  errors  of  about 
+  0.03  for  the  axial  ratios  of  the  drops  up  to  4mm  diameter. 

In  a  series  of  comparisons  of  Chilbolton  radar  data  with  a 
ground  based  disdrometer  (Goddard  and  Cherry,  1984;  and  Goddard, 
Cherry  and  Bringi,  1983),  significantly  improved  agreement  be¬ 
tween  the  radar  and  disdrometer  derived  rainfall  rates  was  ob¬ 
tained  if  the  drops  below  3mm  in  size  were  made  slightly  more 
spherical  than  the  shapes  of  Pruppacher  and  Pitter.  The  proposed 
changes  in  drop  shapes  were  within  the  experimental  error  of  the 
Pruppacher  and  Fitter  laboratory  measurements;  for  example  the 
axial  ratio  of  a  1mm  drop  was  changed  from  0.98  to  1.00.  Recent¬ 
ly,  more  accurate  measurements  of  the  shape  of  naturally  occurr¬ 
ing  raindrops  using  airborne  shadowgraph  instruments 
1 Chandrasekar  et  al ,  1988)  have  confirmed  that  the  shapes  in¬ 

ferred  from  the  radar  are  correct  and  that  the  smaller  drops  are 
indeed  slightly  more  spherical  than  suggested  by  Pruppacher  and 


Pitter.  New  laboratory  experiments  (Beard  and  Ochs,  1968)  on  1.2 
and  1.3mm  drops  also  support  this  conclusion. 

In  this  paper  we  direct  attention  to  heavier  rainfall  with 
values  of  Z  above  40dBZ  and  rainfall  rates  generally  above  10mm 
hr  ^ .  The  precise  shapes  of  raindrops  larger  than  4mm  are  uncer¬ 
tain,  but  are  important  for  Z^p  measurements  in  heavier  rain. 
Doviak  and  Zrnic  (1984)  state  that  the  maximum  value  of  Zpp  due 
to  rain  should  be  3.5dB.  This  assertion  is  based  upon  Green's 
(1975)  formula  for  the  axial  ratio  of  raindrops,  which  is  very 
close  to  the  measurements  of  Pruppacher  and  Pitter.  Using  these 
shapes  a  6mm  drop  would  be  associated  with  a  Z^p  of  3 . 8dB  and  an 
8mm  drop  4.7dB;  summing  the  weighted  contribution  of  the  various 
drops  in  a  Marshal  1 -Palmer  raindrop  size  distribution  leads  to 
the  conclusion  that  the  ZDp  of  rain  cannot  exceed  3.5dB.  Doviak 
and  Zrnic  (1984)  suggest  that  higher  values  of  Z^r  result  from  a 
measurement  related  problem  or  are  due  to  targets  other  than 
raindrops  and  suggest  they  arise  from  ice  particles  which  have 
very  oblate  shapes  when  they  melt. 

Values  of  Zpp  above  3 . 5dB  and  as  high  as  6dB  have  been  re¬ 
ported  in  convective  showers  (Illingworth  et  al,  1987)  and  in 
view  of  the  above  argument  it  seems  logical  to  attribute  such 
observations  to  melting  ice  particles.  On  closer  examination 
other  possibilites  arise.  Based  on  a  statistical  analysis  of  the 
high  values  of  Zgp  accompanying  intense  echoes  at  altitudes  where 
rain  is  to  be  expected,  Caylor  and  Illingworth  i  1986  ;  hypothe¬ 
sised  that  raindrops  larger  than  4mm  are  more  oblate  than  predic¬ 
ted  by  Pruppacher  and  Pitter.  This  suggestion  was  prompted  by 
two  findings.  Laboratory  measurements  'Rasmussen  et  al,  1984  ) 
indicate  that,  unless  the  ice  particles  are  initially  extremely 
oblate  then  the  shapes  during  melting  are  n.  :no  i  e  distorted  than 


those  of  the  same  sized  water  drop.  Secondly,  a  limited  series  of 
aircraft  measurements  of  large  drops  (Cooper  et  al,1983)  implied 
more  oblate  forms.  Beard  and  Chuang  (1987)  have  recently  pub¬ 
lished  more  accurate  computations  of  the  shape  of  large  drops  and 
in  Section  2  of  this  paper  we  use  these  new  shapes  to  calculate 
the  values  of  to  be  expected  for  various  raindrop  size  dis¬ 
tributions.  The  predicted  variation  of  the  average  values  of  Zqr 
as  a  function  of  Z  using  these  new  shapes  is  compared  with  the 
actual  average  values  observed  "’ith  the  Chilbolton  radar.  The 
variation  in  the  naturally  occurring  spectra  around  these  mean 
values  is  discussed  in  Section  3.  The  implications  of  these 
results  for  rainfall  measurement  are  explored  in  Section  4. 
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THE  DIFFERENTIAL  REFLECTIVITY  OF  HEAVY  RAIN 


(a)  Drop  Shapes  and  Differential  Reflectivity 

The  axial  ratios  for  drops  above  4mm  proposed  by  Pruppacher 
and  Pitter  (1971)  are  compared  with  the  more  recent  values  ob¬ 
tained  by  Beard  and  Chuang  (1987)  in  Table  1.  The  values  of 
Pruppacher  and  Pitter  have  considerable  scatter  when  experimental 
-rrors  are  considered,  with,  for  example,  drops  between  7  and 
8mm  having  ratios  in  the  range  0.51  to  0.69.  The  uncertainties 
in  axial  ratios  quoted  for  Beard  and  Chuang ' s  model  are  less 
than  0.015.  For  raindrops  larger  than  2mm  Beard  and  Chuang  find 
that  the  axial  ratio  is  a  nearly  linear  function  of  size  with  a 
fall  in  axial  ratio  of  about  0.06  for  every  millimeter  increase 
in  size.  The  9mm  drop  shape  entry  in  the  table  is  taken  from  the 
graph  in  their  paper  and  is  consistent  with  the  shape  measured  by 
Pruppacher  and  Beard  (1970)  in  the  low  turbulence  UCLA  wind 
tunnel.  The  axial  ratio  for  a  10mm  diamn' drop  is  derived  by 
extrapolation  of  the  linear  relationship  and  is  not  supported  by 
direct  evidence. 

Table  1  show-,  the  different  predicted  magnitudes  of  Zrjp  for 
the  two  sets  of  drop  shapes  using  the  simple  Ray 1 e i gh-Gans 
tr.ecry.  This  is  strictly  only  applicable  to  the  small  drops,  and 
the  Mie-Gans  theory  should  be  used  when  the  drop  diameter  becomes 
a  reasonable  fraction  of  the  radar  wavelength.  We  are  indebted  to 
Dr  Hoit  of  the  Department  of  Mathematics,  Essex  University  for 
the  Mie-Gans  calculations  in  Table  1,  which  were  carried  out  for 
water  drops  at  -iOJC,  O^C  and  +10°C  and  the  Chi  1  bo  1  ton  frequency 
or  .  J  765GHz .  Beard  i  1  9  76)  found  that  changes  in  pressure  and 
temperature  have  ^  negligible  effect  on  raindrop  snape  over  ail 
reasonable  values  to  be  found  in  the  atmosphere .  The  values  of 
•  Mp  using  the  exact  theory  are  only  si  ight  ly  uiqt.e:  :  nan  for 


Rayleigh  scattering  apart  from  drops  above  8mm  where  the  diver¬ 
gence  is  considerable.  Similarly,  we  note  that  temperature 
effects  on  scattering  for  a  given  drop  shape  are  negligible  for 
drops  smaller  than  8mm. 

Rasmussen  et  al  (1985)  report  laboratory  experiments  of  drop 
shapes  in  the  presence  of  a  vertical  electric  field.  Appreciable 
changes  in  axial  ratio  occurred  only  for  a  field  above  200  kv/m 
and  drops  4mm  larger  than  diameter;  the  measurements  did  not 
encompass  drops  larger  than  5mm  or  field  orientations  other  than 
the  vertical.  It  is  generally  believed  that  the  field  required 
to  trigger  lightning  in  the  atmosphere  is  about  300  kV/m; 
regions  where  the  field  exceeds  200  kV/'m  are  likely  to  be 
localised  in  both  time  and  space,  and  are  also  more  likely  to 
occur  where  the  hydrometeors  are  in  the  ice  phase.  The  effects 
of  electric  fields  on  Zqr  are  likely  to  be  small.  However,  if 
the  lightning  does  occur  in  a  region  with  large  raindrops,  a 
localised  sudden  change  in  Zpp  might  occur.  Such  a  step  change 
in  ZpR  might  be  detected  if  the  radar  was  not  scanning  but  was 
dwelling  on  the  lower  regions  of  an  active  storm.  In  view  of  the 
changing  sensitivity  of  differently  sized  large  raindrops  to 
field  changes  apparent  in  the  results  of  Rasmussen  et  a  1 . ( 1 9 8 5  )  . 
a  quantitative  interpretation  of  any  such  Zpp  jumps  will  probably 
not  be  possible. 

(b)  Statistical  Analysis  of  Radar  Data 

A  simple  prediction  of  the  values  of  Zqr  to  be  expected  in 
rain  is  displayed  in  Figure  1  for  a  Marshal  1 -Palmer  raindrop  size 
distribution  with  a  value  of  N'c  of  8000m  'mm  1  .  Tiie  two  curves 
are  obtained  using  the  two  drop  shape  models  in  Table  1  for 
different  values  of  D0  in  Equation  2  at  C°C  and  terminating  the 
raindrop  spectrum  at  8mm.  The  curves  diverge  for  values  of  Z 


above  40dBZ;  the  Pruppacher  and  Pitter  shapes  lead  to  a  maximum 
predicted  value  of  3.5dB,  whereas  the  more  oblate  shapes  suggest 
values  of  up  to  4 . 5dB  are  possible. 

The  crosses  in  Figure  1  represent  the  average  values  of  Z^r 
for  each  2dBZ  step  in  Z,  obtained  from  measurements  in  rain  made 
with  the  Chilbolton  radar  during  summer  convective  storms  on  many 
different  days.  Isolated  echoes  smaller  than  2km  in  diameter  were 
not  included  in  the  analysis.  During  the  observation  period,  the 
freezing  level  was  at  3km  or  slightly  above,  and  so,  to  minimise 
ambiguities  due  to  melting  ice  particles,  data  were  only  an¬ 
alysed  for  altitudes  below  1km. 

It  is  well  established  that  in  severe  storms  with  intense 
reflectivities,  values  of  ZqR  close  to  zero  can  extend  down  to 
the  ground;  these  signatures  are  interpreted  as  hail  reaching  the 
ground  (e.g.  Illingworth  et  al ,  1986).  These  features  are  usually 
quite  localised  and  easy  to  identify,  and  so  any  cell  having 
values  at  1km  altitude  where  ZqR  was  below  ldB  with  Z  above  40dBZ 
was  excluded  from  the  analysis.  An  example  of  the  average  values 
of  Zqp  for  every  2dBZ  step  in  Z  for  such  a  cell  is  shown  in 
Figure  2,  where  the  effect  of  the  hail  in  reducing  the  mean 
values  of  ZqR  for  values  of  Z  above  55dBZ  is  quite  dramatic. 

The  ZDR  values  for  the  averaged  data  in  Figure  1  are  higher 
than  expected  and  it  has  been  argued  that  values  of  Z^R  above 
3 . 5dB  should  be  attributed  to  raindrops  containing  melting  ice 
cores  (e.g.  Doviak  and  Zrnic,  1984;  Beard  and  Chuang , 1987  )  .  The 
evidence  from  the  Chilbolton  radar  does  not  support  this  view. 
Values  of  ZDR  of  up  to  9dB  have  been  observed  in  the  bright  band 
and  are  thought  to  result  from  melting  snowflakes  and  ice  crys¬ 
tals,  but  these  only  occur  in  stratiform  clouds  with  a  clearly 
defined  melting  layer.  In  vigorous  convective  storms  our  observa- 


tions  of  many  vertical  profiles  of  Z  and  Zpp  indicate  that,  as 
ice  particles  melt,  the  value  of  Zj)p  gradually  increases  from 
the  zero  value  associated  with  graupel  and  hail,  and  reaches  a 
maximum  value  reflecting  the  equilibrium  shape  of  the  raindrop 
when  melting  is  complete.  A  typical  profile  is  shown  in  Figure 
3.  The  Zqp  profiles  in  convective  clouds  do  not  show  a  maximum 
associated  with  the  final  stages  of  melting  followed  by  a 
collapse  to  a  less  oblate  raindrop  containing  no  ice.  Labora¬ 
tory  studies  of  the  melting  of  ice  support  this  suggestion; 
Rasmussen  et  al  (1984)  show  that  ice  spheres  smaller  than  9mm  do 
not  shed  water,  but  that  as  melting  proceeds  they  become  pro¬ 
gressively  more  oblate,  finally  assuming  the  shape  of  the  equiva¬ 
lent  diameter  water  drop.  For  spheres  larger  than  9mm  diameter, 
Rasmussen  et  al  observed  that  an  unstable  torus  of  liquid  water 
built  up  around  the  equator.  However,  these  ice  spheres  were 
tethered  within  a  wind  tunnel  but  a  particle  falling  in  the  free 
atmosphere  would  be  able  to  tumble  and  would  undoubtedly  shed 
such  a  torus. 

Reflectivity  data  below  250m  altitude  were  also  excluded 
from  the  analysis  in  Figure  1,  as  these  were  likely  to  be  affec¬ 
ted  by  ground  clutter.  Even  with  this  height  restriction  the 
quarter  degree  beamwidth  Chilbolton  radar  still  detected  the 
occasional  ground  return,  but  the  Zqr  of  such  returns  is  easy  to 
recognise,  as  the  values  of  Zqr  lie  in  the  range  +5dB  to  -5dB  and 
change  by  several  dB  between  one  300m  long  gate  and  the  next 
i Hall  et  al,  1984).  To  suppress  this  clutter,  a  ZDp  reading  was 
rejected  if  the  change  in  Zjjp  (  dB )  at  a  neighbouring  gate  was 
mm  p  than  60  . 

The  average  values  of  Z^p  obtained,  after  the  removal  of 
ground  clutter  and  ambiguities  due  to  melting  ice,  are  those 


displayed  in  Figure  1.  These  data,  together  with  the  sample  size 
and  standard  deviations  are  summarised  in  Table  2.  If  we  assume  a 
normal  distribution  then  the  standard  errors  in  the  mean  are  only 
significant  for  the  values  of  Z  in  the  two  highest  ranges;  these 
error  bars  are  plotted  in  the  figure. 

The  use  of  the  Beard  and  Chuang  shapes  improves  the  agree¬ 
ment  between  the  model  and  the  measurements,  and  we  shall  use 
these  shapes  from  now  onwards.  The  errors  in  axial  ratio  of 
0.015  are  equivalent  to  a  change  in  Zqr  for  the  raindrops  of 
0.2dB,  which  is  still  larger  than  the  estimated  error  of  O.ldB  in 
the  Zqp  measurement  made  by  the  Chilbolton  radar. 

(c)  Analytic  Forms  of  the  Raindrop  Spectrum 

We  have  few  direct  observations  of  the  concentration  of  the 
largest  raindrops.  The  ground-based  disdrometer  described  by  Joss 
and  Waldvogel  (1967)  is  widely  used  for  radar  validation  experi¬ 
ments,  but  the  largest  drop  size  that  can  be  recorded  is  5.3mm; 
even  if  the  size  range  was  increased,  a  negligible  number  of 
these  larger  drops  would  be  detected  because  of  the  small  sample 
area  of  the  instrument.  The  data  displayed  in  Figure  1  were 
obtained  for  storms  within  100km  range  of  the  Chilbolton  radar, 
but  the  number  of  such  storms  which  would  be  above  a  single 
immobile  ground-based  disdrometer  would  be  extremely  small. 
This  restriction  does  not  apply  to  airborne  instruments,  and 
Wiliis  (1984)  reports  raindrop  spectra  measured  inside  a 
hurricane  where  the  rain  rate  was  computed  to  be  ib9mm  hr.  but 
unfortunately  no  data  are  presented  for  drops  greater  than  4mm. 
This  device  sampled  2  nr  of  cloud  every  lOsecs  t  1  .  3km  of  flight 
tracki;  in  the  heaviest  rain  it  sampled  about  20  raindrops  in  the 
4  to  5 mm  size  range;  extrapolating  the  exponential  distribution 


indicates  that  such  a  device  would  not  obtain  a  significant 
sample  of  8mm  drops.  Recently,  Beard  et  al  (1986)  have 
reported  finding  large  raindrops  within  warm  shallow  clouds  over 
Hawaii.  On  one  occasion  the  drop  diameter  reached  8mm  and  the 
average  concentration  of  raindrops  in  the  5  to  6mm  range  for  the 
observations  was  about  one  per  cubic  meter.  It  ‘seems  that 
earlier  suggestions  that  collision  induced  break-up  would  limit 
the  maximum  raindrop  size  to  a  few  millimeters  have  been 
overemphasised. 

Although  it  appears  that  large  raindrops  may  well  be  found 
in  heavy  rain,  and  that  such  large  raindrops  would  affect  the 
values  of  Z  and  Zqr,  it  seems  very  difficult  to  obtain  an  in¬ 
dependent  means  of  sampling  such  drops  to  check  any  radar  infer¬ 
ences.  We  shall  adopt  a  different  approach,  and  test  various 
forms  of  the  raindrop  size  distribution  to  see  which  agrees  most 
closely  with  the  radar  data.  To  compute  the  curves  in  Figure  1, 
we  assumed  a  Dmax  of  8mm,  a  constant  N0  of  8000  m-^mm  and 
allowed  D0  to  vary.  We  shall  also  explore  the  effect  of  fitting 
the  data  to  a  gamma  distribution  (Ulbrich  and  Atlas,  1984): 

N ( D )  =  Nq  Dm  exp( - ( 3 . 67+m )D • D0 )  (3) 

The  third  parameter,  m,  reflects  the  breadth  of  the  size  distri¬ 
bution.  The  use  of  simplified  size  distributions  is  partially 
justified  by  the  averaging  of  microscale  features  by  the  radar 
compared  to  disdrometer  distributions. 

Let  us  first  consider  the  effect  of  truncating  the  exponen¬ 
tial  distribution.  Both  Joss  and  Waldwogei  (1969;  and 
Sachidananda  and  Zrnic  (1988)  used  a  maximum  value  of  6mm  but,  in 
view  of  the.  above  discussion,  larger  values  such  as  the  9mm  used 
in  the  Z_r;R  calculations  of  Goddard  and  Cherry  11984)  should  be 


more  realistic. 


Drops  of  9mm  are  stable  in  a  low  turbulence  wind 


tunnel  (Pruppacher  and  Beard,  1970)  and  so,  in  the  absence  of 
collisions  with  other  raindrops,  it  is  reasonable  to  suppose  they 
would  exist  in  natural  clouds.  In  Figure  4  values  of  Z  and  Zqr 
are  plotted  for  maximum  drop  sizes  of  8  and  10mm,  with  the  value 
of  N0  kept  at  8000  m  ^mm  ^  throughout .  The  agreement  of  the 
predictions  using  the  10mm  diameter  cut-off  with  the  average 
radar  data  is  remarkable. 

Figure  4  also  explores  the  effect  of  the  value  of  m  in  the 
gamma  function.  When  m  =  0  the  gamma  function  reduces  to  the 
simple  exponential,  bur  some  measurements  indicate  that  natural 
raindrop  spectra  contain  rather  fewer  of  both  the  very  large 
drops  and  the  very  small  drops  than  would  be  the  case  for  an 
exponential,  and  that  a  value  of  m  =  2  is  more  appropriate 
(Ulbrich,  1983).  Again  N0  is  kept  constant  at  8000  m  ^mm  ^ ,  and 
the  values  for  Dmax  =8  and  10mm  with  m=2  are  also  plotted  in 
Figure  4.  The  agreement  with  observation  is  not  as  good;  the 
lower  predicted  values  of  Zqr  at  high  rain  rates  reflect  the 
lower  concentrations  of  the  large  raindrops  in  the  narrower 
spectrum . 

If  four  variables  are  used  in  the  raindrop  spectrum  then, 
for  only  two  observables,  an  infinite  combination  of  values  -is 
possible.  There  is  no  reason,  per  se,  to  believe  that  the 
analytic  form  of  the  raindrop  distribution  should  be  the  same  for 
ail  rain  rates,  for  example,  N!0  can  be  adjusted  to  give  the 
required  value  cf  Z  for  any  value  of  Zqr.  In  the  absence  of  any 
other  radar  observable  parameters,  the  most  sensible  choice  for 
summertime  convective  clouds  in  the  UK  seems  to  be  a  pure  ex- 

_  :  _  i 

ponentrai  with  a  variable  D0,  but  keeping  \'Q  at  8000  m  "mm  1  and 
DmdX  -  6  or  10mm.  A  polynomial  which  fits  these  theoretical 


computations  and  gives  Z  correct  to  0.2dBZ  is  given  in  Table  3a. 

Figure  5  compares  the  predicted  values  of  Z^r  as  a  function 

of  0o  for  a  pure  exponential  raindrop  spectrum.  The  polynomials 

which  fit  the  curves  for  the  Beard  and  Chuang  shapes  to  an 
0  C  S 

accuracy  of  O^jS^B  in  are  given  in  Table  3b.  For  comparison 
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purposes  we  also  plot  the  relationship  [Zqr  =  0.45  D0  ] 
suggested  by  Seliga  et  al.(1986)  for  Green's  (1975)  shapes. 

Because  this  eguation  is  constrained  to  pass  through  the  origin 
and  can  have  no  points  of  inflexion,  it  cannot  represent  the 
dependence  very  accurately. 

So  far  we  have  used  the  values  of  computed  at  0°C.  Some 

values  at  +10°C  and  -10°C  are  presented  in  Table  1.  Although 

these  computations  are  not  complete,  and  the  change  in  Z^p  for 
the  full  Mie-Gans  theory  is  not  a  simple  function  of  temperature, 
the  values  in  Table  1  indicate  that  there  would  only  be  a  small 
change  in  the  curves  in  Figure  1  over  this  temperature  range. 
The  values  at  -10°C  are  of  interest  when  interpreting  transient 
narrow  columns  of  positive  Zqr  which  extend  up  to  2km  above  the 
freezing  level  in  vigorous  convection.  We  have  argued 

i  Illingworth  et  al ,  1987  )  that  these  are  due  to  supercooled  rain¬ 
drops  . 

These  effects  of  temperature  at  S  band  (10cm)  should  be 

small,  but  seem  to  be  more  serious  at  C  band  (5.5cm).  Computa¬ 

tions  at  5.5cm  wavelength,  but  using  the  Pruppacher  and  Pitter 
drop  shapes  (Meischner.  Bringi  and  Jank,  lQ88i,  show  a  pronounced 
resonance  for  6mm  drops  with  a  maximum  Z^p  of  near) y  8dB  at  20  C, 
failing  to  4dB  for  7  and  8mm.  The  maximum  is  much  less  marked  at 

0  C.  The  use  of  the  Beard  and  Chuang  shapes  will  modify  these 

values,  but  they  indicate  that  severe  problems  could  arise  in  any 
quantitative  analysis  of  high  ZDp  regions  observed  at  C-band. 


3. 


VARIABILITY  OF  NATURAL  RAINDROP  SPECTRA 


In  the  previous  section  we  concentrated  upon  obtaining 

average  values  of  ZqR  for  a  given  magnitude  of  Z.  Particular 

observations  of  Z  and  ZqR  do  not,  of  course,  lie  on  this  average 

curve;  it  is  this  fact  which  enables  ZqR  to  be  used  as  an 

independent  observable.  The  curve  in  Figure  1  for  the  Beard  and 

Chuang  drop  shapes  represents  the  best  available  estimate  of  the 
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values  of  Z  and  ZqR  for  a  value  of  N0  =  8000  m  mm  in  the 
exponential  raindrop  size  spectrum  of  Equation  (2).  For  any 

particular  observation  of  Z  and  ZqR,  the  value  of  D0  can  be 

derived  from  ZqR,  and  N0  found  from  Figure  1,  using  the  fact  that 
Z  scales  linearly  with  N0 . 

Histograms  of  the  frequency  of  occurrence  of  different  values 
of  ZqR  for  each  interval  of  2dBZ  in  Z  are  plotted  in  Figure  6. 
The  mean  values  of  ZqR  of  each  histogram  are  those  plotted  in 
Figure  1.  The  width  of  the  histograms  can  be  used  to  estimate 
tne  variability  of  naturally  occurring  raindrop  spectra.  Histo¬ 
grams  for  values  of  Z  less  than  30dBZ  are  not  presented  because 
the  resolution  of  obtainable  by  the  radar  is  insufficient; 

histograms  for  Z  above  60dBZ  are  not  shown  because  the  sample 
sizes  are  too  small.  An  asymmetry  of  the  distributions  in  Figure 
6  is  evident,  with  a  long  tail  showing  that  values  of  ZpR  many 
times  the  mean  occasionally  occur.  This  implies  that  values  of 
N’o  are  generally  close  to  8000  m  ^mm  ^ ,  but  that  on  occasion  much 
lower  values  occur.  Figure  4  shows  that  non-exponential  drop 
size  spectra  such  as  gamma  functions  will  also  contribute  to  the 
observed  variability.  but  that  the  long  tail  in  ZDR  cannot  be 
explained  without  invoking  low  drop  concentrations.  We  have 
drawn  attention  (Illingworth  et  al.  1^87)  to  the  anomalously  high 
values  of  7. dp  observed  in  some  isolated  first  echoes,  but. 


because  of  their  limited  spatial  extent,  such  data  are  not 
included  in  the  'average'  data.  The  unusual  characteristics  of 
these  first  echoes  may  be  re-emphasised  with  reference  to  Figure 
6.  In  these  isolated  clouds  Z  values  of  30  and  40dbZ  were  accom¬ 
panied  by  Zpp  values  of  3  and  4dB,  respectively;  the  statistical 
analysis  presented  in  Figure  6  shows  that  for  these  two  values  of 
Z,  98%  of  the  data  in  the  sample  were  accompanied  by  Zpp  values 
of  a  lower  magnitude. 

In  view  of  the  skewness  of  the  Z^p  histograms  in  Figure  6, 
the  histograms  were  re-computed  for  log(Zppj.  The  results  are 
displayed  in  Figure  7  and  a  greater  symmetry  of  the  distributions 
is  apparent.  This  is  confirmed  by  Table  2  which  provides  the 
values  of  standard  deviations,  skewness  and  kurtosis  for  the 
histograms  in  Figures  6  and  7.  The  skewness  (defined  as  the 
ratio  of  the  third  moment  to  the  cube  of  the  standard  deviation) 
of  the  log(ZQp)  distributions  is  about  one  tenth  that  for  the 
linear  ZDp  histograms  for  the  range  of  Z  from  30  to  50dBZ;  above 
50dBZ  the  difference  is  less  clear,  but  this  may  be  affected  by 
the  smallness  of  the  sample  size.  The  kurtosis  (the  ratio  of  the 
fourth  moment  to  the  fourth  power  of  the  standard  deviation) 
should  be  3  for  a  normal  distribution;  higher  values  arise  if 
the  wings  of  the  distribution  are  wider  than  for  a  normal  distri¬ 
bution.  Again,  the  values  of  kurtosis  for  the  logfZpp) 
histograms  in  the  Z  range  30-50dbZ  are  close  to  3,  whereas  the 
linear  Z^p  histograms  have  much  greater  values.  Kurtosis  values 
for  Z  above  60dBZ  are  not  meaningful  because  of  the  small  sample 
size.  The  data  in  Table  2  are  for  clouds  on  many  different  days. 
No  systematic  differences  'were  apparent  when  the  observations  of 
individual  clouds  were  considered. 

The  statistics  in  Table  2  of  the  variability  of  spectra  may 


be  of  use  in 
tical  studies 
in  the  lower 


interpreting  conventional  radar  data  and 
of  problems  related  to  electromagnetic 
troposphere . 


for  statis- 
propagation 


THE  VARIABILITY  OF  RADAR  DERIVED  RAINFALL  RATES 


4  . 

(a)  Derivation  of  Rain  Rate 

In  this  section  we  compare  various  proposed  relationships 
for  deriving  rainfall  rates  when  both  Z  and  Zqr  are  available. 
Figure  8  compares  the  values  of  rain  rate  (R)  as  a  function  of  Z 
and  ZpR  using  the  original  shapes  of  Pruppacher  and  Pitter  and 
the  new  shapes  of  Beard  and  Chuang  for  an  exponential  size  dis¬ 
tribution.  For  a  given  ZpR  the  values  of  R  scale  linearly  with 
Z.  On  the  Figure  curves  for  10  and  100  mm./ hr  are  plotted.  The 
new  drop  shapes  lead  to  slightly  lower  rainfall  rates  for  a  given 
Z  ZqR  combination;  for  example,  for  a  ZqR  of  3  and  3.7  dB  the 
reductions  are  26%  and  444  respectively.  The  effect  of  using  a 
Dmax  of  10mm  is  marked  with  the  dashed  line;  this  uncertainty 
over  choice  of  truncation  limit  only  becomes  significant  for  ZqR 
above  3dB,  at  4dB  the  reduction  in  deduced  rain  rate  is  about 
25  . 

An  analytic  approximation  to  the  theoretical  curves  in 
Figure  8  for  a  rain  rate  of  1  mm  hr  1  using  the  Beard  and  Chuang 
shapes  modified  for  small  drops,  is  given  in  Table  3c.  Measure¬ 
ment  errors  in  ZpR  affect  the  precision  with  which  rain  rates  may 
be  estimated.  Because  of  the  changing  gradient  of  the  curves  in 
Figure  8,  the  errors  in  the  rainfall  estimate  will  increase  -for 
lower  values  of  ZqR.  When  ZqR  is  zero  it  provides  only  an  upper 
bound  on  the  mean  drop  size,  and  so,  for  a  given  Z,  the  ZpR  in¬ 
formation  merely  sets  a  lower  limit  on  the  rainfall  rate.  The 
data  presented  from  the  Chilbolton  radar  in  this  paper  have  a 
standard  error  of  about  O.ldB  in  ZqR. 

Several  other  theoretical  relationships  have  been  published 
recently  which  relate  the  rainfall  rate  to  the  observed  values  of 
Z  and  Zpp  and  these  are  compared  in  Figure  3.  The  curves  are  for 


a  rain  rate  of  1  mm  hr  ^ ,  and,  as  in  Figure  8,  the  rain  rates 
scale  linearly  with  Z.  Seliga  et  al .  (1986,  Equation  35)  use 
Green's  (1975)  drop  shapes  and  suggest  for  ZqR  in  dB : 

Z(dBZ)  =  27.10  +  10.40  log(ZDR)  ( 0 . 2<ZDR<0 . 7 )  (4a) 

Z(dBZ)  =  27.98  +  16.70  log(ZDR)  ( 0 . 7<ZDR<2 . 6 )  (4b) 

Because  we  now  believe  that  small  drops  are  slightly  more 
spherical  than  Green's  formula.  Equation  4  overestimates  the 
rainfall  by  about  1 . 5dB  (40%)  for  low  values  of  ZDR,  and  con¬ 
versely  for  large  values  of  ZqR  the  Green  shapes  are  insuffic¬ 
iently  oblate  and  the  rain  rate  is  underestimated  by  about  1 . 5dB 
(30%).  Seliga  et  al .  analyse  the  drop  size  distributions 
measured  during  a  three  hour  heavy  rainfall  event  and  compute  the 
theoretical  values  of  Z  and  ZqR  wuich  would  have  been  observed, 
and  show  that  using  the  two  observables  Z  and  ZqR  should  give  a 
better  estimate  of  rainfall  rate  than  an  empirical  Z-R  relation¬ 
ship.  This  argument  demonstrates  the  validity  of  the  technique 
but  does  not  prove  that  the  drop  shapes  used  for  Equation  4  are 
correct . 

Z  and  ZqR  radar  data  for  part  of  a  single  15km  x  15km  PPI 
scan  are  analysed  by  Sachidanada  and  Zrnic  (1987).  The  data  are 
compared  -with  the  theoretical  relationship  for  an  exponential 
drop  size  distribution  with  a  maximum  size  of  6mm  using  Green's 
(1975)  drop  shapes.  For  a  rain  rate  of  1mm  hr  *  their  Equation  S 
with  ZqR  in  dB  may  be  expressed  as 

Z( dBZ  )  =  21.65  +  4.86  ZDR  ! 5  ! 

They  claim  that  this  linear  relationship  fitted  the 
theoretical  curve  to  within  3  ,  or  0.05dB;  in  view  of  the  other 
theoretical  curves  in  Figure  9,  this  seems  difficult  to 


unde r st and . 


They  comment  that  the  observations  or  the  scatter  of 


with 


Z  and  Zqr  for  the  single  PPI  scan  seem  large  when  compared 
the  scatter  simulated  by  generating  random  values  for  the  drop 
size  distribution  parameters  in  the  ranges: 

30  <  Nq  <  30,000  mm  ^  m  ^ 

-2  <  m  <  +  3 

4  <  Dmax  <  6  ( mm ) 

0.27  <  D0  <  1.23  (mm) 

Accordingly,  the  observations  from  one  part  of  the  scan 
which  seemed  to  have  a  particularly  high  scatter  were  discounted, 
although  we  note  that  the  highest  values  of  D0  used  in  the 
simulation  is  equivalent  to  a  rain  rate  of  only  6mm  hr-1  for 
Marshal  1 -Palmer  rain.  Noticing  that  the  mean  value  of  Z  for  a 
given  ZqR  for  the  remaining  data  seemed  rather  higher  than 
expected,  they  suggested  that  the  raindrops  were  very  oblate  and 
proposed  a  relationship  between  the  axial  ratio  and  drop  diameter 
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a  b  =  1  -  0.64  0*'“  (6) 

A  slight  dip  in  the  mean  values  of  Z  for  higher  values  of  ZDR  was 

interpreted  in  terms  of  drop  break  up.  For  the  more  oblate 

shapes  they  suggested  a  revised  relationship  for  R  =  1mm  hr  ^ : 

Z  =  22.39  +  3.49  ZDR  (7) 

Our  experience  with  the  Chilbolton  data  is  that  the  averages  from 
a  single  PPI  are  highly  variable  and  that  such  generalisations 
embodied  in  Equations  6  and  7  are  not  justified. 

(b)  Variability  of  Rain  Rates  for  a  Given  Z 

In  order  to  obtain  an  estimate  of  the  variability  of  natural 
rainfall  for  a  given  Z,  the  histograms  of  Chilbolton  data  for  Z^R 
i n  Figure  6  have  been  replotted  in  Figure.  10.  The . po 1 ynomi a  1  in 

Table  3c  has  been  used  to  convert  the  values  of  ZqR  to  a  rainfall 
rate.  The  data  are  summarised  in  Table  4,  .vhere.  for  a  given  Z, 


the  mean  ram  rate  and  its  standard  deviation  are  presented. 


We 


note  that  the  standard  deviation  of  the  inferred  rain  rate  is 
somewhat  greater  than  half  the  mean  rate,  but  that  the 
distributions  are  skewed  so  that,  on  occasion,  the  inferred  rain 
rate  is  many  times  that  expected  from  a  simple  Z-R  relationship. 
The  data  is  taken  with  a  sample  volume  which  is  typically  a  300m 
cube;  if  this  volume  was  increased  then  averaging  would  be 
expected  to  narrow  the  histograms  in  Figure  10. 

Before  applying  the  relationships  in  Table  2  the  accuracy  of 
the  observations  must  be  considered.  The  Chilbolton  radar 
estimates  Zpp  to  an  accuracy  of  about  O.ldB,  introducing  an  error 
in  rain  rates  of  about  12s  (that  is,  0.5dB)  for  values  of  Zpp 
above  IdB  (Figure  9).  When  ZqR  values  approaching  zero  are 
common,  as  is  the  case  in  Table  4  when  Z  is  below  -iOdBZ,  then  the 
errors  in  the  derived  rainfall  rate  are  larger. 

The  fundamental  accuracy  limiting  the  Zpp  measurement  should 
be  analysed  before  calculating  rainfall  rates.  Even  if  we  assume 
there  are  no  instrumental  errors  due  to  mismatches  of  the 
channels  for  horizontal  and  vertical  po 1 ar l zat ion ,  and  that  the 
beam  is  filled  with  rainfall  with  no  severe  reflectivity 

gradients  across  the  sample  volume,  the  accuracy  of  the  Zpp 
estimate  is  limited  by  the  finite  dwell  time  of  the  radar. 
Bringi  et  ai .  (1983)  show  that  aoout  40-60  independent  sample 
pairs  of  Zjj  and  Zy  are  required  if  ZDR  is  to  be  measured  with  a 
standard  error  as  low  as  O.ldB.  The  time  to  independence  is 
related  to  the  width  of  the  Dopple-  spectrum  of  the  target  and  is 
usually  in  the  range  20  to  40  msecs  at  10cm  wavelength..  The 
ehiibolton  ZDp  data  are  obtained  with  a  time  integration  of  210 
msecs  and  spatial  integration  over  four  "’5m  range  gates.  If 

raster  scanning  is  employed  with  no  spatial  avei aging  then  the 


number  of  independent  samples  will  fall.  If  the  number  of  sam¬ 
ples  is  only  6,  for  example,  the  standard  error  in  Zqr  will  be 
about  0.5dB.  Even  for  large  values  of  Zqr  this  would  lead  to  an 
error  of  nearly  lOOi  in  the  derived  rain  rate. 


DISCUSSION  AND  CONCLUSIONS 


5  . 

The  precise  shape  of  raindrops  which  should  be  used  for  the 
interpretation  of  Zpp  measurements  has  been  a  matter  of  con¬ 
siderable  debate.  The  radar  results  reported  in  this  paper  are 
in  good  agreement  with  the  with  the  new,  more  accurate,  shapes 
proposed  for  the  large  drops  by  Beard  and  Chuang  (1987)  and  the 
Goddard,  Cherry  and  Bringi  (1983)  shapes  for  drops  below  3mm.  In 
the  absence  of  any  third  observable  reflecting  the  value  of  m  in 
the  gamma  distribution,  the  most  reasonable  spectrum  for  inter¬ 
preting  the  Z  and  Zpp  observations  is  the  exponential  of  Equation 
(2)  with  a  Dmax  =  8mm.  In  heavy  rain  with  Z  values  above  50dBZ 
there  is  an  indication  that  a  Dmax  =  10mm  could  be  used. 

Based  upon  these  new  drop  shapes  and  for  an  exponential 
distribution  of  raindrop  sizes,  the  analytic  expressions  in  Table 
3  appear  most  appropriate  for  calculating  D0  from  the  observed 
value  of  Zqp  and  for  deriving  N0  and  the  rainfall  rate  from  Z  and 
Zr)p  data.  Particular  attention  should  be  paid  to  errors  intro¬ 
duced  into  Zgp  by  sampling  limitations. 

In  Illingworth  et  al .  ( 1987)  it  was  argued  that  values  of 
Zpp>  of  6dB  were  due  to  a  monodispersed  distribution  of  large 
raindrops.  and  the  new  drop  shapes  discussed  in  this  paper 
support  their  interpretation.  From  a  radar  analysis,  Cayicr  and 
Illingworth  (1986)  hypothesised  that  large  raindrops  must  be  more 
oblate  than  was  commonly  accepted;  the  more  oblate  shapes  have 
now  been  confirmed  independently,  and  give  more  consistent  agree¬ 
ment  with  the  radar  data.  However,  the  extreme  oblateness  of  the 
very  large  drops  proposed  by  Cay lor  and  Illingworth  i  1936  )  are 
•:  v.  exaggerated,  and  those  or  Beard  and  1'huang  :  1  ?>>7  >  should  be 
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The  average  values  of  differential  reflectivity  for  eacn  2d8Z 
step  in  Z  observed  m  heavy  rain  by  the  Chilboiton  radar.  The 
lert-hand  half  of  the  table  displays  the  mean,  standard  devia¬ 
tions  'Oi,  skewness  and  kurtosis  of  the  distributions  of  Znp 
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TABLE  4 
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The  average  values  of  radar- infer  red  rainfall  rate  ' R  ;  for  each 
2dBZ  step  in  Z,  observed  in  heavy  rain  by  the  Chiiboiton  radar. 
Histograms  of  the  data  are  shown  in  Figure  10.  Op.  is  the  standard 
deviation  and  n  the  sample  size. 
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Figure  1 

The  variation  of  Z  with  Z^p  for  an  exponential  drop  size  distri¬ 
bution  using  two  drop  shape  models.  The  solid  line  is  for  the 
Beard  and  Chuang  (1987)  drop  shapes,  and  the  dashed-dot  line  is 

for  the  Pruppacher  and  Pitter  (1971)  shapes.  In  both  cases  Dmax 

-  3  - 1 

=  8mm  and  No  =  8000m  mm  .  The  crosses  are  averaged  Chilbolton 
radar  data  collected  from  summer  convective  clouds. 

Figure  2 

An  example  of  the  effect  due  to  hail  showing  the  significant 
decrease  in  Z^p  at  high  Z  values.  X  represents  the  averaged 
Chilbolton  data  from  two  PPIs  through  a  storm  on  6  July  1983, 
where  a  hail  shaft  was  identified.  The  two  lines  are  theoretical 
Z-Zpp  relationships  and  are  identical  to  those  in  Figure  1. 

Figure  3 

A  typical  vertical  profile  through  a  vigorous  convective  cloud 
showing  the  monotonic  increase  in  Zpp  as  the  ice  particles  fall 
through  the  freezing  level  at  3  km  and  subsequently  melt. 

Figure  4 

The  relationship  of  Z  to  Zqr  for  variations  in  the  parameters  of 
the  gamma  drop  size  distribution.  The  solid  lines  are  m  =  0, 

Dmax  =  8mm  and  Dmax  =  10mm.  The  dotted  line  is  m  =  2,  Dmax  = 
8mm  and  the  dashed  line  is  m  =  2,  Dmax  =  10mm.  The  Eeard  and 
Chuang  drop  shapes  were  used  and  No  =  8000m  ^mm  * .  The  crosses 
are  averaged  Chilbolton  data  as  in  Figure  1. 


Figure  5 


A  comparison  of  the  various  proposed  relationships  between  Zpp 
and  D0. 

Solid  line:  Beard  and  Chuang  drop  shapes  for  Dmax  =  8  and  10mm. 

Dash-dot:  Pruppacher  and  Fitter  shapes  for  Dmax  =  8mm. 

Dashed  line:  Zqr  =  0.45  Do*"'**’  (Seliga  et  al.l986j. 

Figure  6 

Histograms  of  the  number  of  Zqr  values  associated  with  2dB 
intervals  in  2.  The  data  were  observed  for  summer  convective 
clouds  in  1983  and  1984.  The  number  in  the  corner  of  each  histo¬ 
gram  indicates  the  centre  of  the  particular  2  interval  in  units 
of  dBZ . 

Figure  7 

The  histograms  of  Figure  4  replotted  with  logarithmic  spaced  Zqr 
bins. 

Figure  8 

The  variation  of  Z  with  expected  for  constant  rainfall  rates. 

The  lower  set  of  curves  is  for  a  rainfall  rate  of  10mm  hr  and  the 
upper  set  is  for  a  rate  of  100mm  hr.  The  solid  line  is  the  Beard 
and  Chuang  drop  shapes  and  an  exponential  drop  size  distribution 
with  Dnax  =  8mm,  while  the  dashed  line  is  for  the  Dmax  =  10mm. 

The  dotted  line  shows  the  rainfall  expected  with  the  Pruppacher 
and  Fitter  shapes  for  an  exponential  distribution  with  Dmax 
-  8  mm. 


Figure  9 


A  comparison 
rainfall  rates 

i  ,  - 1 

1mm  hr 

Solid  lines: 

Dotted  line: 
Dashed  lines: 


of  various  suggested  relationships  for  deriving 
from  Z  and  Zqr.  The  curves  are  for  rain  rates  of 

This  study  using  the  modified  Beard  and  Chuang 
shapes . 

Seliga  et  al.(1986),  Green's  shapes,  Equation  4. 
Sachidananda  and  Zrnic  {1987),  Green's  shapes. 
Short  dashes  -  Equation  5;  Long  dashes  -  Equation 
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Figure  10 

Histograms  of  the  rainfall  rates  associated  with  2dB  intervals  in 
Z.  The  data  are  from  Figure  4,  using  Table  3  to  compute  the 


rainfall  rate  from  Z  and  Zpp. 
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I.  INTRODUCTION 

The  S  band  Chi  1  bo  It  on  radar  in  the  UK  has  a 
quarter  degree  beamwidth  and  can  transmit 
pulses  (separation  1.6msec)  which  are  alter¬ 
nately  polarised  in  the  horizontal  and  verti¬ 
cal  direction.  For  co-polar  reception  we  can 
obtain  the  radar  reflectivity  factors  for  the 
two  polarizat ions  ,  ZH  and  ZV,  and  then  derive 
ZDR,  the  differential  reflectivity  (.10  log 
’ZH  ZV ) ) .  ZDR  is  a  wasure  of  mean  shape, 
which  for  raindrops  gives  drop  size.  We  shall 
also  discuss  results  of  the  time  series  data 
obtained  by  recording  the  power  received  from 
each  transmitted  pulse.  The  cross-polar  re¬ 
turn,  ZHV  ,  and  the  linear  depolar izat ion 
LDR= 10 log( ZHV ,ZV) )  can  also  be  obtained, 
which  provides  information  on  particle  fail 
mode.  The  Chilbolton  antenna  can  detect  values 
oi  LDR  as  low  as  -32dB.  Results  1  reported 
elsewhere'  show  how  the  canting  angle  of  rain¬ 
drops  and  the  mean  axial  ratio  of  rumbling 
graupel  may  be  derived  from  LDR . 

J.  DIFFERENTIAL  REFLECTIVITY,  ZDR 

The  application  of  2  and  ZDR  to  estimate  both 
*ne  size  and  c oncent r at l on  of  raindrops  is 
well  <nown  e.g.  ILLINGWORTH  et  a: ,  198'  '  . 

Until  recent  Iv  PRUPPACHER  and  FITTER  19'Z' 
prcvided  the  best  estimate  of  drop  shapes, 
email  adjustmnets  to  these  shapes  now  appear 
appropriate.  F  r  m  ZDR  radar  data  -TODOARD  et 
a:  argueu  tnat  millimetre  drops  are 

slightly  more  spherical;  this  suggestion  has 
been  confirmed  ov  CHANDRASEKAR  et  al  ;  1986  . 
T?  •>xplam  ZDR  va  I  les  measured  in  neavy  rain, 
IA R  .iiiu  ILLINGWORTH  i  1 9 15 7  postulated  that 
i:  .-ps  larger  man  *tnm  must  he  more  .'Plate  than 
t n e  *pp a : ner  and  fitter  shapes.  BEARD  and 

n  *NG  .  4  h '  nav«*  recent  iv  cinfirme:  that 
larger  ir 'ps  ire  more  'plate.  Direct  me  as  .  re  - 
~+“ '  -  s  o:  -*d  .i 1 1  i  br  ;  am  ir  >p  -napes  ir^  now  ;n 
wr.'-n-n*  with  ra-tar  inferences,  and  anv  nat 
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3.  TIME  SERIES  CORRELATION  OF  ZH  AND  ZV 
Figure  1  shows  time  series  data  obtained  at  a 
single  75m  gate  in  210  milliseconds  from  64 
pulse  pairs  transmitted  with  alternate  hori¬ 
zontal  and  vertical  po lar i sat  ions .  Successive 
estimates  fluctuate  as  the  particles  reshuffle 
in  space.  ZDR  is  normally  derived  by  com¬ 
paring  the  linear  averages  of  these  64  esti¬ 
mates,  as  shown  by  the  straight  lines  in  the 
Figures.  Figure  1  shows  a  zero  ZDR  for  small 
spherical  raindrops.  Additional  information  is 
available  from  p(H,V),  the  correlation  of  ZH 
and  ZV,  which  here  is  equal  to  0.99;  for 
spheres  the  theoretical  value  is  unity.  Some 
statistical  properties  of  p(H,V)  are  discussed 
by  3RINGI  et  ai  (  1983) . 
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ptH,V)  will  be  less  than  one  it  the  scattering 
amplitudes  for  the  two  polarizations  from  each 
particle  arrive  at  the  antenna  with  a  non 
constant  amplitude  ratio  or  with  a  phase  dif¬ 
ference.  A  distribution  of  differently  sized 
particles  will  be  the  most  important  effect  in 
reducing  p(H,V),  A  p(H,V)  of  0.432  in  the 
melting  layer  (Figure  2)  confirms  this.  The 
bright  band  in  Z  is  quite  pronounced  in  the 
vertical  profiles  plotted  in  Figure  3;  the 
maximum  value  of  Z  is  lodBZ  higher  than  in  the 
rain  below,  indicating  chat  Che  particles  must 
be  of  low  density  ice  and  are  probably  snow¬ 
flakes.  The  ZDR  profile  shows  a  bright  band  in 
ZDR  200m  below  that  in  Z.  This  depression  of 
the  ZDR  bright  band  is  a  common  occurrence  in 
stratiform  clouds,  and  increases  as  Z  becomes 
larger;  the  melting  snowflake  apparently 
reaching  its  maximum  degree  of  oblateness 
after  it  has  started  to  collapse.  The  low 
value  of  ?(H,V)  in  the  bright  band  is  probably 
caused  by  the  coexistence  of  half-melted  snow¬ 
flakes  and  smaller  raindrops. 

A  spectral  analysis  (Figure  4)  of  the  time 
series  in  the  melting  laver  (Figure  2)  shows 
an  interesting  component  from  40-60Hz  which  is 
greater  in  ZH  Chan  ZV,  and  is  not  observed  in 
tne  time  series  for  the  snow  above  or  the  rain 
below  the  bright  band. 

We  are  currently  developing  Monte  Carlo  mo¬ 
dels  t>  predict  the  reduction  in  p(H,V)  for 
the  following  hydrometeors: 

a''  Monod  i  s  persed  ablates  -  small  phase  dif¬ 
ferences  between  the  amplitude  pair  from 
-ach  particle  due  to  losses  or  particle 
asymme t  ry . 

a  Monod i spe rsed  tumbling  aolates 
' e . g.  graupe 1 ) . 

;  ? aly oi sper sed  f  ambling  aolates 

e  .  g.  graupel  . 

d ?  recess  mg  and  nutating  polydispered 
ablates  'e.g.snow  and  melting  snow1. 

e  Differently  shaped  aligned  oblates 
(  e  .  g .  ra  i  n  '  . 

In  ram  with  a  ZDR  of  ).5d3  we  have  found  a 
significant  tall  of  p ( H , V  J  to  about  ) . 9  5 ; 
Out  nave  not  yet  *st all  shed  if  p'H,V)  :an  ae 
ised  as  an  independent  e  timate  at  tne 
r^ait’i  or  the  size  ( i  r  i  out  .  on  . 

The  '..'west  val  les  at  p '  Vi  .  ■  in  the  i:*  m 

'uvective  clouds,  shoal  i  ?ccur  when  tne  -;f— 
terence  m  the  horizontal  and  vertical  ;:ar- 
termg  imp  1 1 :  ud.- s  or  the  particle  19  largest; 
f n ; s  na -  we  11  ac  :u r  when  slightly  as pne r 1 ; a  1 
car:  »•:»-*  enter  the  Mie  region,  and  :  oul  i  -.«►  a 
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Figure  2.  Time  series  pair  observed  in  the 
melting  layer  of  a  stratiform  cloud  showing 
a  lower  correlation. 


v«r*tlCll  A'OM;*  0*  Sfl  3«CHMr  1M'  *t  06430*  V 
’•tf*  ’  *  Ml*  «  J  u  •  1 34  XMMg  Oat  ■  M  it 


9  94* 


Figure  3.  Profile  of  the  vertical  structure 
of  a  stratiform  ram  cloud  from  which  the  rime 
series  in  Figures  i  and  2  are  taken. 
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4.  ZDR  FLARE  ARTEFACTS  AND  IMPLICATIONS  FOR 
DOPPLER  ANALYSIS 

WILSON  and  RE ‘CM  (  1938)  analyse  Che  Doppler 
properties  of  a  radar  artefact  called  a 
"flare  echo",  which  extends  downrange  of  some 
incense  radar  storm  echoes.  Ic  is  caused  by 
triple  scattering  from  a  high  Z  region  at 
height  h,  down  to  the  ground,  then  back  to  the 
radar  via  the  precipitation,  thus  forming  a 
spurious  echo  a  distance  h  behind  the  intense 
core.  This  flare  has  a  Doppler  velocity  of 
equal  magnitude  to  chat  of  the  core  but  with  a 
reversed  sign. 

Figure  5,  from  the  Chi  1  bo  It  on  radar,  has  a 
flare  echo  where  ZDR  reaches  *9dB.  At  a  range 
of  37km  Z  reaches  70d3Z,  and  the  flare  may  be 
identified  by  the  column  of  positive  ZDR  in¬ 
clined  at  -*5deg  to  the  vertical,  extending 
from  the  grouno  at  38km  range  to  an  altitude 
ot  ukm  at  a  range  of  91km.  We  believe  that 
this  high  value  of  ZDR  arises  because  che 
triple  scattering  of  the  flare  echo  is  con¬ 
fined  to  che  ZH  channel;  for  the  vertical 
polarization  the  induced  dipoles  within  the 
cloud  should  not  radiate  along  their  3xes  down 
to  the  ground. 

Clearly  such  values  of  *9dB  cannot  oe  identi¬ 
fied  m  terns  of  hvdr one teors ,  but  m  less 
obvious  cases  positive  values  of  ZDR  towards 
the  back  of  high  Z  regions  cou’d  be  erroneous¬ 
ly  interpreted  as  hvdrometeors.  ZDR  artefacts 
caused  by  flare  echoes  should  be  much  more 
wi  iespread  ana  intense  at  C  and  X  band  than  at 
5  band. 

Doppier  measurement  s  are  usually  made*  using 
the  ZH  chjnnel  .  but  these  arguments  suggest 
mat  Doppler  proolens  associated  with  flare 
echoes  would  be  reduced  if  tne  ZV  cnannel  was 
used.  An  analysis  of  the  MAYPOLE  iata,  ob¬ 
tained  Dv  the  SCAR  J P 2  S-banu  radar,  reveals 
mat  the  Dippler  derived  velocities  from  ZH 
and  ZV  are  normally  identical  apart  tram  flare 
regions.  Fortner  studies  are  und^rwav  to  see 
if  such  argurw»nts  are  valid  at  C  and  X  band. 
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ABSTRACT 

We  present  radar  observations  of  precipitation  par¬ 
ticles  using  a  new  polarization  parameter,  p(H,V) , 
which  is  the  correlation  of  the  time  series  of  ZH 
and  ZV  {the  radar  reflectivity  factors  measured  with 
horizontal  and  vertical  polarization).  p(H,V)  is 
very  close  to  unity  in  most  precipitation.  In 
stratiform  precipitation  low  values  of  p(H,V)  are 
round  only  in  the  melting  layer,  and  so  may  be  used 
as  a  simple  means  of  identifying  the  bright  band. 

We  suggest  that  in  convective  clouds  p(H,V)  will  be 
low  only  where  hailstones  are  large  enough  to  enter 
the  Mie  scattering  region.  Our  observations  are 
made  at  S-band  ,  but,  because  p(H,V)  should  be  little 
affected  by  propagation  problems,  the  technique 
could  oe  implemented  at  shorter  wavelengths. 

Keywords:  Polarization,  T ime -Ser ies ,  Differential- 

Reflectivity,  Bright-Band,  Hail,  Correlation 


The  5  sand  Chilbolton  radar  in  the  UK  has  a  quarter 
aeiiiree  beamwiath  and  can  transmit  pulses  (separa¬ 
tion  l.otnsec!  which  are  alternately  polarized  m 
tne  horizontal  and  vertical  direction.  From  t'ne  co- 
p:iar  reception  we  can  obtain  the  radar  reflec¬ 
tivity  factors  for  the  two  polarizations,  ZH  and  ZV, 
ana  tn er.  derive  ZDR,  the  differential  reflectivity 
‘■10  log  (ZH  ZV)).  ZDR  is  a  measure  of  mean  shape, 
wmcn  for  raindrops  gives  drop  size.  The  applica¬ 
tion  of  Z  and  ZDR  to  estimate  both  the  size  and 
concentration  at  raindrops  is  well  known;  interpre¬ 
tation  of  the  ZDR  of  ice  vs  more  difficult  (Hall  et 
a .  ,  19  8*-.,  Illingworth  et  al  ,  1987',  but  measurements 
m  nature  convective  clouds  with  the  Chilbolton 
radar  consistently  show  ZDR  values  within  O.ldB  of 
zero.  We  shall  discuss  results  of  the  time  series 
data  obtained  by  recording  the  power  received  from 
each  transmitted  pulse. 

Figure  1  shows  time  series  data  obtained  at  a  single 
on;  gate  in  210  milliseconds  from  pulse  pairs 
transmitted  with  alternate  horizontal  and  vertical 
polarizations.  Successive  estimates  fluctuate  as  the 
particles  reshuffle  in  space.  ZDR  is  normally 
derived  by  comparing  the  linear  averages  of  these  6* 
estimates,  as  shown  by  the  straight  lines  in  Figure 
i,  in  tru>  case  ZDR  is  zero  for  the  small  spherical 
raindrops  in  the  light  rain.  Additional  information 


is  available  from  p(H,V),  the  correlation  of  ZH  and 
ZV,  which  here  is  equal  to  0.99;  for  spheres  the 
theoretical  value  is  unity.  Some  statistical  proper 
ties  of  p(H,V)  are  discussed  by  BR1NGI  et  al  (1983/ 
Sachidananda  and  Zrnic  (1985)  show  theoretically 
that  values  of  p(H,V)  in  rain  should  be  close  to 
unity . 


Figura  1.  Tiaa  Sanaa  In  light  rain  on  16/12/87 
Z^iS.MdBZ  Zon— O.OldB  rHV-0.991  tK-28.8aaac 


FlQura  2.  Tiaa  Sanaa  In  bright  band  on  16/12/07 
ZH"2l.23d8Z  Zoq-3.34dB  r^jy-0.432  tH-22.4aaac 
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One  factor  reducing  p(H,V)  will  be  if  the  various 
particles  in  the  sample  volume  have  different  ratios 
of  scattering  cross-sections  for  the  two  polariza¬ 
tions.  A  distribution  of  particles  with  different 
shapes  will  be  the  most  important  effect  in  reducing 
p(H»V).  A  p(H,V)  of  0.432  in  the  melting  layer 
(Figure  2)  confirms  this.  The  bright  band  in  Z  is 
quite  pronounced  in  the  vertical  profiles  plotted  in 
Figure  3;  the  maximum  value  of  Z  is  l5dBZ  higher 
than  in  the  rain  below,  indicating  that  the  parti¬ 
cles  must  be  of  low  density  ice  and  are  probably 
snowflakes.  The  ZDR  profile  shows  a  bright  band  in 
ZDR  200m  below  that  in  Z.  This  depression  of  the 
ZDR  bright  band  is  a  common  occurrence  in  stratiform 
clouds,  and  increases  as  Z  becomes  larger;  the 
melting  snowflake  apparently  reaching  its  maximum 
degree  of  (electromagnetic)  oblateness  after  it  has 
started  to  collapse.  The  low  value  of  p(H,V)  in  the 
bright  baud  is  probably  caused  by  the  coexistence  of 
half-melted  snowflakes  and  smaller  raindrops. 

The  power  spectra  for  the  time  series  for  ZH  and  2V 
in  rain  are  identical  (e.g.  Figure  4),  but  in  the 
melting  layer  (Figure  5)  there  are  usually  some 
higher  frequency  components  which  are  greater  in  ZH 
than  ZV.  These  may  result  from  the  greater  contri¬ 
bution  of  the  oblate  snowflakes  with  low  terminal 
velocities  leading  to  a  broader  Doppler  spectrum  in 
ZH ,  or,  alternatively,  fluttering  of  the  oblate 
semi-aligned  snowflakes  modulating  ZH  more  chan  ZV, 

Our  measurements  in  the  ice  in  convective  clouds 
generally  show  high  values  of  p(H,V) .  When  values 
of  decorrelation  times  (t  )  fall  below  20  msec  it  is 
necessary  to  interpolate  the  data  so  that  the  ZH  and 
ZV  sampling  vs  effectively  simultaneous.  We  expect 
low  values  of  p(H,V)  to  be  confined  to  regions 
containing  large  hail.  Ground  observations  of  hail 
show  that  the  axial  ratio  rarely  exceeds  0.9.  Lab¬ 
oratory  scattering  experiments  on  such  particles 
which  are  large  enough  to  enter  the  Mie  region, 
reveal  that  their  scattering  cross-sections  should 
change  by  up  to  lOdB  as  they  tumble;  this  should 
give  rise  to  low  values  of  p(H,V). 

The  melting  layer  observations  lend  support  to  this 
proposed  hail  detection  algorithm.  We  observe  high 
p(H,V)  in  heavy  rain  with  a  ZDR  of  2 . 5dB  although  a 
distribution  of  drop  shapes  must  be  present,  but  an 
equally  high  ZDR  in  the  melting  layer  is  accompanied 
by  low  p(H,V).  We  are  investigating  theoretically 
if  this  implies  such  a  large  range  of  scattering 
amplitudes  chat  we  need  to  invoke  Mie  scattering 
from  the  larger  half  melted  snowflakes  in  the  bright 
band  . 
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1 .  INTRODUCTION 

Polarization  radar  measurements  can  provide 
information  on  the  characteristics  of  precipita¬ 
tion  particles  which  is  unavailable  from  the 
simple  reflectivity  measurement.  We  discuss 
means  of  measuring  raindrop  sizes  and  concentra¬ 
tions,  raindrop  canting  angles,  and  we  outline 
techniques  for  differentiating  between  the 
various  forms  of  frozen  hydrometeors  (snow, 
graupel  and  hail)  and  for  rapid  identification 
of  the  bright  band. 

The  S-band  (10cm)  Chilbolton  radar  (Goddard 
and  Cherry  ,  1987)  situated  in  Hampshire  in  the 
UK  has  a  25m  dish  (quarter  degree  beamwidth)  and 
is  the  largest  steerable  meteorological  radar  in 
the  world.  The  radar  can  transmit  (and  receive) 
pulses  every  1.6msec  which  are  alternately  pol¬ 
arized  in  the  horizontal  and  vertical  direc¬ 
tions.  The  pulses  are  0.5usec  long  with  a  peak 
power  of  500kW,  and  the  return  is  digitised 
every  500nsec  to  give  a  75m  range  resolution.  In 
this  paper  we  present  observations  (introduced 
by  Illingworth  and  Cay  lor,  1988a)  of  (a)  The 
differential  reflectivity,  ZDR,  which  measures 
mean  hydrome^eor  shape,  and  (b)  the  co-polar 
correlation  p*(H,V) ,  which  reflects  the  variety 
of  particle  shapes  present.  A  companion  paper 
(Illingworth  and  Caylor,  1989)  discusses  a  third 
parameter,  (c)  the  Linear  Depolarization  Ratio, 
LDR,  which  is  affected  by  the  fall  mode  of  the 
particles.  We  shall  also  report  results  of  the 
cross-polar  correlation  and  the  frequency  spec¬ 
tra  of  the  co-  and  cross-polar  time  series. 

2.  DEFINITION  OF  TERMS  AND  EXAMPLES  OF  DATA 

Figure  V  shows  time  series  data  for  210 
msecs  obtained  at  a  single  75m  gate  from  64 
pulse  pairs  transmitted  with  alternate  horizon¬ 
tal  and  vertical  polarizations.  Reception  is  co- 
polar  (or  parallel)  to  the  transmitted  polariza¬ 
tion  so  the  64  samples  are  of  ZH  and  ZV,  the 
reflectivities  for  horizontal  and  vertical  pol¬ 
arization  respectively.  Successive  estimates  of 
ZH  and  ZV  fluctuate  as  the  particles  re-arrange 
in  space,  the  decorrelation  time  of  each  time 
series  being  related  to  the  width  of  the  Doppler 
spectrum  of  the  scattering  particles. 

The  differential  reflectivity,  ZDR,  is 
de  '  ned  as 

ZDR  *  10  log  (ZH/ZV)  (1) 

where  ZH  and  ZV  are  the  averages  (linear  power) 
of  the  64  samples,  as  shown  by  the  straight 
lines  in  Figure  l;  this  data  is  for  heavy  rain 


and  has  a  ZDR  of  2.76dB.  ZDR  is  particularly 
useful  for  distinguishing  ice  from  rain  and  for 
measuring  the  mean  shape  of  the  raindrops  (see 
Section  3). 


Time  Series  on  13  July  1988  at  155425  UT 

Ray-148  Gate-30  Az-253.0deg  El-2.4deg  R-l3.5km 

Zj^-46 . 05d8Z  ZDR-2.76dB  rHV-0.992  tH- 19.2msec 


Figure  1 

ZH  and  ZV  time  series  taken  in  heavy  rain. 

Additional  information  can  be  derived  from 
the  correlation  of  the  two  time  series  in  ZH  and 
ZVj  which  we  call  the  co-polar  correlation 
(p  (H,V)).  Previously  the  average  values  of  ZH 
and  ZV  have  been  computed  in  real  time  and  the 
time  series  data  has  not^been  recorded.  For  the 
heavy  rain  in  Figure  1  p  is  0.99  (  see  Sections 
4  and  5). 

3.  ZDR,  THE  DIFFERENTIAL  REFLECTIVITY 

In  rain  ZDR  provides  a  measure  of  the  mean 
raindrop  9hape.  Raindrops  become  more  oblate 
with  increasing  size,  so  ZDR  is  positive  and 
the  magnitude  of  ZDR  is  related  to  the  mean 
raindrop  size.  The  application  of  Z  and  ZDR  for 
estimating  both  the  size  and  concentration  of 
raindrops  is  well  known  (e.g.  Seliga  and  Sringi , 
1976).  Until  recently  Pruppacher  and  Pitter 
(1971)  provided  the  best  estimate  of  drop 
shapes.  Small  adjustments  to  these  shapes  now 
appear  appropriate.  From  ZDR  radar  data  Goddard 
and  Cherry  (1983)  argued  that  millimeter  sized 
drop9  are  slightly  more  spherical  than  previous¬ 
ly  believed;  this  suggestion  has  been  confirmed 


bv  the  aircraft  measurements  of  drop  shapes  by 
Chandrasekhar  et  al  (1988)  and  the  laboratory 
measurements  of  Beard  and  Ochs  (1988).  To  ex¬ 
plain  high  ZDR  values  measured  in  heavy  rain 
Cay  lor  and  Illingworth  (1987)  postulated  that 
drops  larger  than  4mm  must  be  more  oblate  than 
the  Pruppacher  and  Pitter  shapes;  Beard  and 
Chuang  (1987)  have  recently  confirmed  that  lar¬ 
ger  drops  are  indeed  more  oblate.  Direct  meas¬ 
urements  of  equilibrium  drop  shapes  are  now  in 
agreement  with  radar  inferences,  and  any  natur¬ 
ally  occuring  drop  oscillations  do  not  appear  to 
bias  the  radar  measurements.  An  analysis  of 
these  drop  shapes  and  their  effect  on  estimates 
of  drop  size  distributions  and  rainfall  rates 
is  given  by  Illingworth  and  Caylor  (1988b).  Z 
and  ZDR  measurements  indicating  that  first 
echoes  in  warm  convective  clouds  consist  of  a 
low  concentrat ion  of  large  raindrops  are  dis¬ 
cussed  by  Illingworth  (1988). 

The  theoretical  precision  of  the  ZDR  meas¬ 
urement  is  analysed  by  Bringi  et  al  (1983),  and 
Figure  2  (from  their  equation  7)  predicts  how 
the  number  of  independent  samples  and  the  co- 
polar  correlation,  p,  limit  the  accuracy  of  ZDR. 
For  the  tin*  aeries  in  figure  1,  there  are  10 
independent  samples,  p  (H,V)  i9  0.99,  so  the 
predicted  standard  error  in  the  ZDR  estimate  is 
0.25dB.  To  increase  the  accuracy  of  rainfall 
estimates,  ZDR  must  be  measured  to  an  accuracy 
of  O.ldB.  Figure  2  implies  that  40  independent 
samples  are  needed  if  p*0.99;  the  Chilboiton 
radar  achieves  this  by  spatial  averaging  over 
four  adjacent  75m  range  gates. 

Interpretation  of  the  ZDR  of  ice  is  more 
difficult  (Hall  et  al,  1984,  Illingworth  et  al 
1987)  because  ZDR  varies  with  axial  ratio, 
shape,  fall  mode  and  dielectric  constant.  Figure 
3  shows  the  values  of  ZDR  as  a  function  of  axial 
ratio,  for  Rayleigh  scattering  from  oblate 
spheroids  of  water  and  ice  with  their  major  axes 
aligned  in  the  horizontal.  Our  experience  with 
the  Chilboiton  radar  is  that  ZDR  values  in  con¬ 
vective  clouds  are  within  O.ldB  of  zero,  in¬ 
dicating  tumbling  ice,  or,  at  the  other  extreme, 
aligned  graupel  with  a  mean  axial  ratio  of  less 
than  0.95.  In  stratiform  clouds  the  values  of 
ZDR  in  ice  are  variable,  and  although  zero  is 
the  moat  common,  values  of  up  to  3d8  occur, 
particularly  where  Z  tends  to  be  low.  High  ZDR 
values,  locally  up  to  6dB,  can  be  found  in  the 
melting  layer. 

4.  MEASUREMENT  PROBLEMS  FOR  THE  CO-POLAR 

CORRELATION 

Two  practical  problems  arise  in  estimating 
the  correlation:  the  finite  signal  to  noise 
ratio  (SNR)  and  the  non-simultaneous  sampling  of 
ZH  and  ZV.  Bringi  et  al  (1983)  show  that  for  a 
finite  SNR  the  measured  value  p*  is  related  to 
the  true  value  of  p  by: 

p*  *  p  (1  ♦  1/SNR)°  <1  ♦  ZDR/SNR)0-5  (2) 

For  values  of  ZDR  close  to  OdB,  a  true  value  of 
p-I,  is  reduced  to  0.999  foe  an  SNR  of  30dB, 
and  0.99  for  20dB.  The  reduction  of  the  correla¬ 
tion  due  to  non-simultaneous  sampling  is  given 
by  (Zrnic  et  al ,  1988): 

p*  *  p.  p(l>  (3) 


where  p(l)  is  the  decorrelation  of  the  H  or  the 
V  tin*  series  after  one  pulse  delay.  In  prac¬ 
tice,  because  of  measurement  errors,  Equation  3 
can  occasionally  lead  to  corrected  values  of 
p(H, V)  which  are  slightly  larger  than  unity.  We 
prefer  to  use  a  third  order  polynomial  fit  and 
interpolation  to  estimate  the  coincident  values, 
although  it  is  computationally  more  intensive. 
Sorting  of  many  hundreds  of  time  series  in  rain 
at  altitudes  between  250m  and  l  km  and  with  a  SNR 
better  than  20dB,  reveals  that  for  a  decorrela¬ 
tion  time  (tau)  of  13msec  the  average  value  of 
p*2  is  0.95,  but  rises  to  0.97  after  interpola¬ 
tion;  the  average  improvement  using  equation  3 
is  similar..  For  a  tau  of  25maec  the  average 
value  of  p*  is  0.97  which  rises  to  0.98  after 
interpolation;  for  longer  values  of  tau  correc¬ 
tions  for  non-simultaneous  sampling  are  not 
significant . 

The  drop  in  correlation  due  to  Doppler 
effects  will  be  more  serious  at  shorter  wave¬ 
lengths,  although  this  can  be  mitigated  by  the 
use  of  a  higher  PRF,  The  extreme  narrowness  of 
the  Chilboiton  beam  minimises  Doppler  broadening 


Independent  Samples 

Figure  2 

Standard  deviation  of  ZDR  a a  a  function  of  the 
number  of  independent  samples  and  correlation. 


Figure  3 

ZDR  as  a  function  of  axial  ratio 
for  oblate  spheroids. 


due  to  wind  shear,  and  data  are  always  taken 
with  antenna  elevations  below  15  .  If  elevations 
any  higher  are  used  then  the  Doppler  spread 
introduced  because  of  the  raindrops*  terminal 
velocities  will  lead  to  very  short  decorrelation 
times  and  compromise  any  measurments  of  the  co- 
polar  correlation. 

In  their  theoretical  study  Sachidananda  and 
Zrnic  (1985)  predict  that  the  correlation  in 
rain  should  be  at  least  0.995.  Jameson  and  Dave 
(1988)  and  Bebbington  et  al  (1987)  have  demon¬ 
strated  that  the  linear  correlation  may  be  de¬ 
rived  from  circular  polarization  measurements. 
In  this  case  simultaneous  reception  of  the  two 
circularly  polarized  is  possible,  but  it  is  not 
obvious  how  the  accuracy  of  the  derived  value 
for  linear  correlation  is  limited  by  the 
transformation  from  circular  to  linear  and  by 
the  algorithms  which  correct  for  propagation 
effects . 

5.  CO-POLAR  CORRELATION  RESULTS 

The  measured  values  of  linear  co-polar 
correlation  in  most  precipitation  are  not  signi¬ 
ficantly  different  from  unity.  Theoretically  in 
rain  with  a  large  ZDR  there  should  be  a  reduc¬ 
tion  in  correlation  due  to  the  variety  of  drop 
shapes.  We  fiod  tha^  when  ZDR  is  2.5dB  the 
average  value  of  the  p  is  about  0.97,  but  these 
high  ZDRs  tend  t«'  occur  in  more  turbulent  condi¬ 
tions,  and  it  is  not  clear  how  much  of  the 
reduction  is  due  to  the  inadequacy  of  our  cor¬ 
rection  algorithm.  Ideally,  the  drop  in  correla¬ 
tion  should  br  related  to  the  width  of  the 
raindrop  spectra,  and  could  provide  a  value  of  m 
if  a  gamma  function  is  used  to  describe  the 
raindrop  spectrum,  but  any  effect  seems  to  be 
masked  by  measurement  errors. 

An  examp 1*  of  a  time  series  measured  in 
snow  i9  jhow-  in  Figure  4.  In  this  case  the 
value  of  p*  reaches  the  exceptionally  high  value 
of  0.999,  aide'!  no  doubt  by  the  SNR  of  50dB  for 
these  data.  The  only  region  where  appreciable 
reductions  ocir  is  in  the  fright  band  as 
demonstrated  in  Figure  5  where  p  falls  to  0.31. 
The  low  correlation  io  the  bright  band  is  proba¬ 
bly  caused  by  the  coexistence  of  large  oblate 
wet  snowflakes  and  spherical  raindrops. 

Conventional  radars  scanning  in  P?I  mode 
are  used  to  estimate  rainfall,  but  when  the  beam 
intersects  the  jright  band  the  rainfall  is  over¬ 
estimated.  Smi.i  (1986)  has  suggested  algorithms 
for  bright  baid  identification  which  involve 
comparing  PPI  -  at  two  elevations.  Our  ^ata  in 
stratiform  clo.d  show  that  a  value  of  p  below 
0.8  is  an  excellent  detector  of  the  bright  band; 
the  efficienc  of  the  algorithm  is  easily  veri¬ 
fied  by  examining  RHI  scans  where  the  bright 
band  can  be  clearly  seen  from  the  Z  structure. 
The  use  of  the  correlation  method  has  several 
advantages : 

(a)  It  is  unaffected  by  differential  attenua¬ 
tion  and  so,  in  contrast  to  ZDR,  should  be 

applicable  at  shorter  wavelengths. 

(b)  Only  a  s^ort  time  series  is  required  to 

estimate  p  ,  thus  overcoming  the  long  dwell 

times  needed  to  measure  ZDR. 

(c)  The  correlation  is  unaffected  by  differen¬ 
tial  phase  shifts. 


We  suggest  that  in  convective  clouds  the 
value  of  p  should  fall  wheo  the  hail  particles 
are  large  enough  for  Mie  scattering  to  occur. 
This  idea  is  based  upoa  the  observations  that 
all  hail  particles  are  siightly  aspherical  and 
tumble  a a  they  fall.  Consequently  the  vertical 
and  horizontal  cross  sections  of  the  individual 
tumbling  hailstones  will  vary  by  up  to  15dB 
(Atlas  and  Vexler ,  1963).  At  10cm  Mie  scattering 
should  be  rare,  but  for  a  3cm  radar  the  drop  in 
correlation  should  be  an  indication  that  hail 
has  reached  a  diameter  of  about  1cm.  In  addition 
to  the  advantages  itemized  for  the  melting  band, 
this  algorithm  should  work  equally  well  for  wet 
and  dry  hailstones.  Thus  far  the  maximum  Z  we 
have  observed  in  ice  is  about  40dBZ  (but  see, 
Zrnic  et  al ,  1988). 


Time  Series  on  29  May  1988  at  160418  UT 

Say-51  Gate-62  Az-290.0deg  El-11. 3deg  9-9. 3km 

ZH-25.47dBZ  Z0n-0.0ldB  rHV-0.999  tH-22. 4msec 


Figure  4 

Time  senes  for  ice  above  the  melting  layer. 


Time  Series  on  29  May  1.988  at  160418  UT 

Ray-24  Gate-62  Az-290.0deg  El-5.4deg  R-9.3km 

ZH-25. 43dBZ  ZDR-2.52d8  rHV-0.312  tH-12.9maec 


Figure  5 

ZH  and  ZV  time  series  in  the  bright  band. 


6. 


2DR  FLARE  ARTEFACTS  AND  IMPLICATIONS  FOR 
DOPPLER  ANALYSIS 


Wilson  and  Reum  (1988)  Analyse  the  Doppler 
properties  of  a  radar  artefact  called  a  "flare 
echo"  which  extends  dovnrange  of  some  intense 
radar  storm  echoes.  It  is  caused  by  triple  scat¬ 
tering  from  a  high  2  region  at  a  height  h,  down 
to  the  ground,  then  back  to  the  radar  via  the 
precipitation,  thus  forming  a  spurious  echo  a 
distance  h  behind  the  intense  core.  The  Doppler 
velocity  of  this  flare  is  related  to  the  term¬ 
inal  velocity  of  the  particles  in  the  core. 

Figure  6,  from  the  Chilbolton  radar,  has  a 
flare  echo  where  2DR  reches  *9dB.  At  a  range  of 
87km  2  is  70dB2,  and  the  flare  may  be  identified 
by  Che  column  of  positive  2DR  inclined  at  45  deg 
to  the  vertical,  extending  from  the  ground  at 
88km  range  to  an  altitude  of  4km  at  a  range  of 
91km.  We  believe  that  this  high  value  of  ZDR 
arises  because  the  triple  scattering  of  the 
flare  echo  is  confined  to  the  2H  channel;  for 
the  vertical  polarization  the  induced  dipoles 
within  the  cloud  should  not  radiate  along  their 
axes  down  to  Che  ground. 

Clearly,  such  values  of  +9dB  cannot  be 
interpreted  in  terms  of  hydromeceors ,  but  in 
less  obvious  cases  positive  value  of  2DR  towards 
the  back  of  high  2  regions  could  be  erroneously 
ascribed  to  hydrometeors.  2DR  artefacts  caused 
by  flare  echoes  should  be  much  more  widespread 
and  incense  at  C  and  X  band  chan  at  S  band. 

Doppler  measurements  are  usually  made  using 
Che  2H  channel,  but  these  arguments  suggest  that 
Doppler  problems  would  be  reduced  if  the  2V 
channel  was  used.  An  analysis  of  the  MAYPOLE 
data,  obtained  by  the  NCAR  CP2  S-band  radar, 
reveals  that  the  Doppler  derived  velocities  from 
2H  and  2V  are  normally  identical  apart  from  in 
the  flare  regions  where  there  are  large 
differences . 
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1  .  INTRODUCTION 

In  this  paper  we  present  cross-polar  radar 
measurements  in  precipitation  and  discuss  the 
vertical  profiLes  of  three  polarization  para¬ 
meters  through  the  bright  band.  The  specifica¬ 
tion  of  Che  S-band  Chilbolton  radar  and  an  in¬ 
troduction  to  the  various  polarization  parameters 
is  given  in  a  companion  paper  (Caylor  and 
Illingworth,  1989,  referred  to  as  C&I),  which 
concentrates  on  the  differential  reflectivity 
(ZDR),  and  the  correlation  between  the  co-polar 
received  signals  for  vertically  and  horizontally 
polarized  t ransrai ssion  (p“(H,V)).  We  now  con¬ 
sider  the  additional  information  available  in  the 
cross-polar  (or  orthogonal)  linearly  polarized 
returns  . 

Figure  1  shows  time  series  data  for  210 
msecs  obtained  at  a  single  75m  gate  from  64  pulse 
pairs,  for  each  pulse  the  reception  alternates 
between  the  co-polar  and  the  cross-polar.  The 
linear  depolarization  ratio  LDR  is  defined  as: 

LDR  -  10  log  (ZVH/ZH)  (1) 

where  ZVH  is  the  average  (linear  power)  of  the  64 
estimates  of  the  cross  polar  return  (transmit 
vertical,  receive  horizontal),  and  ZH  is  the  same 
average  for  the  co-polar  (horizontal)  return. 
For  the  time  series  in  Figure  1  the  value  of  LDR 
is  about  -3ldB.  LDR  senses  particle  fall  mode 
and,  as  we  shall  see,  is  an  excellent  detector  of 
wet  ice.  It  is  also  possible  to  estimate  the 
cross-polar  correlation  (pZ(VH,H))  between  the 
two  time  series,  which  for  Figure  1  is  only 
-0.08. 

2.  LINEAR  DEPOLARIZATION  RATIO 

The  earliest  LDR  data  were  taken  by  Browne 
and  Robinson  (1952)  at  wavelengths  of  3.2cm  and 
3mm  who  detected  an  enhanced  cro99-polar  return 
from  the  bright  band.  More  recent  observations  at 
3cm  are  reported  by  (Herzegh  and  Conway,  1986). 
Herzegh  and  Jameson  (1988)  propose  a  correction 
algorithm  for  the  propagation  problems  that  occur 
at  these  wavelengths,  which  are  caused  by  the 
progressive  depolarization  of  the  incident  radia¬ 
tion  as  it  passes  through  preci pi  tat  ion ,  and 
result  in  a  rise  in  the  apparent  values  of  LDR 
with  increasing  range. 

The  first  measurements  of  LDR  at  10cm  wave¬ 
length  were  made  with  the  Chilbolton  radar  in  the 
suaner  of  1988  (Goddard  et  al ,  1988).  At  thi9 
wavelength  propagation  problems  appear  to  be 
minimal,  and  theoretical  interpretation  is  sim¬ 
plified  because  Rayleigh  scattering  may  be 
assumed  for  most  precipitation  particles. 

Precipitation  particles  will  depolarize  in¬ 
cident  radiation  only  when  they  are  aspherical 
and  have  no  axi9  of  symmetry  parallel  to  the 
incident  polarization.  It  we  consider  oblate 


particles,  then  the  strength  of  the  depolarized 
signal  is  dependent  upor.  the  ratio  of  the  polar¬ 
izabilities  along  their  major  and  minor  axes, 
that  is  to  say,  it  is  a  function  of  the  value  of 
ZDR  the  particles  would  have  if  they  were 
aligned.  For  randomly  tumbling  particles  the 
value  of  LDR  is  given  by  (Atlas  et  al ,  1933): 

LDR  *  (1  -  yZDR  ♦  Z DR ) / ( 3  ♦  yZDR  +  8ZDR)  (2) 

where  ZDR  is  the  value  the  particles  would  have 
if  they  were  aligned.  Equation  (2)  is  plotted 
(solid  line)  in  Figure  2.  The  cross  polar  isola¬ 
tion  of  the  Chilbolton  antenna  is  so  good  that 
LDR  levels  down  to  at  least  -32dB  can  be  measured 
(dotted  line  in  the  Figure).  The  values  of  ZDR  as 
a  function  of  axial  ratio  for  oblate  spheroids  of 
water,  solid  ice,  graupel  and  snow  are  given  in 
Figure  3  of  C&I.  If  randomly  tumbling  particles 
are  to  give  a  measurable  LDR,  their  value  of  ZDR 
(if  aligned)  must  be  at  least  0.8dB.  This 
suggests  that  because  of  its  low  dielectric  con¬ 
stant  the  depolarization  of  sn^w  should  be  very 
low;  for  a  density  of  O.lg  cm  the  axial  ratio 
of  tumbling  snow  would  have  to  be  less  than  0.15 
to  be  detectable.  For  tumbling  graupel  (density 
0.5)  to  give  a  measurable  LDR,  an  axial  ratio  of 
less  than  0.73  is  required;  this  may  occur 
occasionally .  As  soon  as  these  particles  become 
wet,  the  dielectric  constant  will  rise  and  the 
cross-polar  return  will  be  much  larger.  The 
values  of  LDR  for  spheroids  with  a  constant 
canting  angle  are  also  plotted  in  Figure  2,  the 
two  dashed  curves  are  for  angles  of  5°  and  10° 
randomly  distributed  about  the  vertical  axis; 
these  values  are  lower  than  for  the  same  parti¬ 
cles  randomly  tumbling,  and  should  be  applicable 
to  raindrops. 

Time  Series  on  13  July  1988  at  155512  UT 

Ray-20  Sate-30  Az-253.0deg  El-2.4deg  R-13.5k» 

ZH-48.37dBZ  LDR— 30.67d8  r^vH*-0-005  tH-9.6awec 


Figure  1.  The  co-  and  cross-polar  time 
series  in  heavy  rain 


Values  of  LDR  in  light  rain  are  limited  by 
the  isolation  of  the  antenna  to  less  than  -32dB, 
with  some  measurements  as  low  as  -34dB.  When  the 
rain  has  a  ZDR  of  above  l,5dB  it  is  possible  to 
detect  a  cross-polar  return  which  tends  to  in¬ 
crease  as  ZDR  rises.  First  results  indicate  that 
values  are  consistent  (Figure  2)  with  a  constant 
canting  angle  of  10  randomly  distributed  about 
Che  vertical  axis,  or,  more  realistically,  for  a 
Gaussian  distribution  of  such  canting  angles  with 
a  standard  deviation  of  about  5°.  Such  low  can¬ 
ting  angles  reduce  the  values  of  ZDR  of  the  rain 
by  less  than  O.IdB  (Goddard  et  al ,  1988).  These 
LDR  values  are  significantly  lower  than  the  -26 
to  -29dB  reported  at  3cm  by  Herzegh  and  Conway 
(1986);  the  3cm  values  could  be  affected  by  pro¬ 
pagation  effects,  lower  antenna  isolation,  or 
non-Ravleigh  scattering  by  the  large  raindrops. 

Values  in  ice  again  tend  to  be  at  the  anten¬ 
na  limit.  During  1988  no  regions  of  ice  with  Z 
exceeding  40dBZ  were  sampled,  but  on  some  occa¬ 
sions  LDR  values  of  -30dB  were  observed,  consis¬ 
tent  with  tumbling  graupel  (density  0.5)  with  an 
axial  ratio  of  about  0.66,  or,  more  likely,  solid 
ice  with  an  axial  ratio  of  0.8.  Much  higher 
values  are  to  be  expected  for  hail  in  wet  growth 
(see  measurements  by  Herzegh  and  Conway,  1986). 
When  hailstones  become  so  large  chat  Mie  scat¬ 
tering  theory  applies,  the  values  of  LDR  should 
also  rise  because,  even  for  axial  ratios  of  0.9, 
the  differences  in  the  polarizability  along  the 
major  and  minor  axes  of  the  hailstone  can  be  a 
factor  of  ten  (Atlas  and  Wexler ,  1963).  By  far 
the  most  interesting  LDR  data  were  obtained  for 
melting  ice  particles  as  described  below. 

3.  VERTICAL  SECTIONS  AND  PROFILES 

The  two  colour  plates  show  RHI  scans  of  Z, 
ZDR  and  LDR.  In  Plate  l  we  believe  the  ice 
particles  to  be  snow,  whereas  in  Plate  2  it  is 
suggested  that  graupel  is  present.  The  data  for 
these  two  plates  are  taken  at  very  short  range  to 
maximise  the  power  in  the  cross-polar  channel.  In 
Plate  l  the  system  noise  at  a  7km  range  is  equi¬ 
valent  to  a  Z  of  about  -20dBZ  and  at  14km  range 
to  -!4dBZ;  if  there  is  to  be  significant  power 
in  the  cross  polar  channel  for  all  types  of 
precipitation,  then  the  Z  value  in  the  main 
channel  should  be  at  lease  15dBZ  and  20dBZ  at 
these  two  ranges,  respectively.  An  additional 
advantage  of  operating  at  short  ranges  is  that  at 
10km  the  beamwidth  is  only  50m.  The  ZDR  and  LDR 
data  are  not  recorded  simultaneously,  but  for 
sequential  scans  are  separated  by  about  30secs; 
comparisons  were  only  made  when  the  Z  values  had 
not  changed  appreciably  in  this  period. 

LDR  measurements  are  very  sensitive  to 
ground  clutter  as  shown  in  Plate  l,  where  the  LDR 
values  below  300m  are  high  and  very  variable  from 
gate  to  gate.  The  ground  scatters  isotropically 
(Skolnik,  1980)  and  only  a  small  ground  component 
is  required  to  dominate  the  cross-polar  signal 
from  the  precipitation.  The  ZDR  data  in  Plate  l 
are  unaffected  because  the  co-polar  signals  are 
so  much  higher;  ground  clutter  is  very  easv  to 
identify  from  ZDR  (Hall  et  al,  1984). 

In  Plate  1  the  LDR  values  in  the  dry  ice  and 
che  rain  (ZDR  up  to  IdB)  are  both  at  the  antenna 
limit.  The  most  striking  feature  is  the  uniformly 
high  values  of  LDR  of  about  -16dB  in  the  melting 


Figure  2.  Values  of  LDR  for  canted  and 


tumbling  oblate  spheroids 


LDR  (dB) 


Figure  3.  Vertical  profiles  at  9.3km  for  the 
RHI  in  Plate  1. 

Vertical  Profile  on  29  May  1988  at  160418  UT 
AZ  -  290 . 00deg  Gate  16  at  range  9.3km 
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layer,  which  are  equivalent  to  wet  tumbling 
spheroids  with  an  axial  ratio  of  0.55.  Figure  3 
depicts  a  vertical  profile  of  Z,  ZDR,  LDR  and  the 
co-polar  correlation  p  (H,V)  at  a  range  of  9.3km 
for  the  data  in  Plate  l.  From  this  and  other 
profiles  we  note  that  Che  bright  band  in  LDR  and 
Z  cover  the  same  constant  height  range,  whereas 
Che  bright  band  in  ZDR  occurs  up  to  5‘  ho  lower, 
with  the  larger  depressions  occuring  where  Z  is 
greater.  The  values  of  correlation  are  very  high 
in  the  ice  (Figure  5  of  C&I  is  the  time  series  at 
a  height  of  1.85km),  the  slightly  lower  values  in 
the  rain  below  are  probably  measurement  related 
and  caused  by  the  shorter  decorrelation  time.  The 
low  levels  of  correlation  are  found  only  in  the 
bright  band  and  coincide  with  the  Z^R  bright 
band,  at  1.3km  LDR  and  Z  are  high  but  p  is  still 
unity.  All  these  data  point  to  the  melting  of  low 
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We  conclude  that  the  particles  are  snow¬ 
flakes,  but  the  evidence  on  the  change  of  fall 
mode  of  the  flakes  as  they  melt  is  harder  to 
disentangle.  When  melting  first  starts  at  1.5km 
height  the  flakes  are  undergoing  a  predominantly 
tumbling  or  spinning  motion  (see  the  lab  studies 
of  Mitra  et  al  (1988)),  although  the  slight  rise 
in  ZDR  to  0.3dB  implie|  that  there  is  a  small 
degree  of  alignment,  p  is  still  unity  because 
the  particles  are  the  same  shape.  At  about  1.2km, 
the  degree  of  alignment  becomes  much  greater ,  as 
the  central  pores  of  the  oblate  snowflakes  start 
to  fill  in  with  water,  and  the  fall  mode  is 
stabilised.  Over  the  next  200m  some  flakes 
collapse  to  spherical  raindrops,  causing  a  fall 
in  Z,  and  p  falls  because  we  have  a  variety 
of  /jlace  and  spherical  shapes.  Subsequently, 
below  1  km  all  the  half-melted  flakes  collapse  to 
more  spherical  raindrops,  ZDR  falls  and  p  rises 
as  the  particles  once  again  all  have  the  same 
shape  . 

2 

The  low  values  of  p  in  the  bright  band 
limit  the  accuracy  of  the  ZDR  measurement.  Even 
chough  the  shorter  decorrelation  time  in  the 
bright  band  allows  us  to  take  20  independent 
samples  in  the  210msec  time  series,  Figure  2  of 
C&I  (from  Bringi  et  al ,  1983)  shows  that  if  p  is 

0.6,  then  the  standard  deviation  of  ZDR  will  be 
about  IdB.  This  should  be  contrasted  with  the 
situation  in  the  rain  or  snow,  where,  although 
only  ten  samples  are  gathered  in  a  single  time 
series,  for  p*0.99  the  ZDR  error  is  only  0.2dB. 
For  the  profiles  in  Figure  3  the  values  of  Z,  ZDR 
and  LDR  are  obtained  by  spatial  averaging  over 
four  adjacent  ?5m  range  gates;  reducing  the  error 
in  ZDR  to  C.ldB  and  Z  to  0.7dBZ. 

The  second  colour  Place  shows  a  very  diff¬ 
erent  situation.  There  is  a  monotonic  increase  in 
ZDR  when  melting  occurs  and  no  maximum  is  evi¬ 
dent,  and  the  LDR  maximum  on  melting  is  about  - 
26dB  wit1;  the  altitude  of  this  maximum  varying  by 
up  to  500m.  The  correlation  is  high  everywhere, 
with  the  slight  reduction  from  unity  in  Che  rain 
again  being  a  measurement  problem.  From  the 
colour  plate  we  note  that  the  LDR  is  finite  in 
regions  of  the  rain  where  ZDR  is  high;  this  has 
been  interpreted  in  terms  of  finite  canting  an¬ 
gles  in  the  previous  section.  Some  areas  of  the 
ice  have  a  finite  return  also,  values  of  -30dB 
indicating  higher  density  more  oblate  particles. 
Figure  4  is  a  vertical  profile  at  a  range  of 
13. 8km  through  the  rain  with  the  highest  ZDR  in 
the  Plate,  where  the  rain  has  a  ZDR  between  2  and 
2 . 5dB  and  LDR  is  about  -29dB.  There  is  strong 
evidence  Chat  the  ice  particles  are  high  density 
slightly  oblate  graupel  particles: 

(a)  An  LDR  of  -26JB  implies  tumbling  wet  parti¬ 
cles  with  axial  ratio  0.83.  Dry  ice  parti¬ 
cles  of  this  shape  would  not  give  a  de¬ 
tectable  LDR  even  for  a  density  of  0.9. 

(b)  ZDR  increases  monotonies  1 ly  -  the  oblate 
graupel  particles  melt  and  progressively 
assume  the  9hape  of  the  equilibrium  rain¬ 
drop,  the  collapse  stage  to  less  oblate 
raindrops  noted  in  Figure  3  is  absent. 

(c)  The  correlation  ;s  high  everywhere  indica¬ 
ting  that  there  is  not  a  mixture  of  parti¬ 
cles  having  very  different  shapes. 

After  correcting  3cm  data  for  propagation 
effects,  Bringi  et  al .  (1986)  derived  an  LDR  of 

-24  to  -26dB  for  melting  graupel. 


Vertical  Profile  on  13  July  1988  at  155425  UT 
AZ  -  253. OOdeg  Gate  31  at  range  13. 8km 


Figure  4.  Vertical  profiles  at  13.8km  for  the 
RHl  in  Plate  2. 


The  values  of  LDR  and  Z  start  to  rise  at 
2.4km,  and  LDR  reaches  it  maximum  value  after 
about  200m,  where  ZDR  is  still  only  slightly 
positive.  The  profiles  are  consistent  with  tum¬ 
bling  graupel  particles  which  stabilise  as 
melting  progresses.  Another  feature,  absent  in 
Plate  1,  is  that  the  altitude  of  the  LDR  bright 
band  varies  by  up  Co  500m,  with  the  maximum 
heights  occurring  where  2  is  greatest.  The  peak 
values  of  LDR  appear  to  be  independent  of  alti¬ 
tude.  The  height  at  which  ZDR  starts  to  increase 
follows  the  same  pattern.  This  variation  of  the 
altitude  of  the  LDR  bright  band  is  difficult  to 
explain.  It  may  be  that  where  Z  is  greater  the 
temperature  of  the  surface  of  the  graupel  is  a 
few  degrees  higher  and  so  wet  growth  starts 
sooner;  this  could  arise  if  the  high  Z  regions 
have  higher  LWC  and  larger  graupel  particles. 

The  different  magnitudes  of  LDR  in  the 
bright  band  should  be  a  reliable  indicator  of  the 
type  of  precipitation  in  the  ice  above.  Although 
the  LDR  measurement  requires  good  antenna  perfor¬ 
mance,  the  very  high  values  of  LDR  in  the  bright 
band  should  be  reasonably  easy  to  detect.  The  LDR 
measurement  itself  can  be  simpler  than  ZDR,  in 
that  no  rapid  switching  of  the  polarization  of 
the  transmitted  power  is  required,  and  the  two 
receive  channels  do  not  need  to  be  so  precisely 
matched . 

We  are  currently  analysing  2D  probe  data 
taken  from  the  Met  Office  C-130  when  flying  near 
or  just  above  the  melting  layer  on  several  days 
in  1988.  On  13  July  1988  the  aircraft  was 
struck  by  lightning  while  penetrating  a  region 
with  a  Z  below  40dBZ.  The  strike  occured  where 
the  LDR  bright  band  indicated  a  6km  region  of 
graupel  embedded  in  a  widespread  area  having  an 
LDR  31*1  ow  bright  band.  This  is  consistent  with 
current  ideas  on  electrification;  cloud  electri¬ 
fication  occurs  where  graupel  and  ice  crystals 
exist,  and  t r iboelect r ic  charging  of  the  aircraft 
by  collisions  with  crystals  should  further  en¬ 
hance  the  electric  field  so  that  a  discharge  may 
be  triggered  by  the  aircraft.  In  regions  of  snow 
bright  band  any  small  ice  crystals  would  tend  to 
be  captured  by  the  snowflakes. 


4.  CROSS  CORRELATION  AND  THE  SPECTRA  OF  THE 
TIME  SERIES. 

A  typical  power  spectrum  of  the  two  co-polar 
time  series  (Fig  5  in  C&I)  in  the  snow  bright 
band  for  the  profile  in  Figure  3  is  displayed  in 
Figure  5.  Higher  frequency  components  are  present 
in  the  ZH  spectrum  than  in  that  for  ZV.  The  value 
of  ZDR  was  2.76dB,  so  ZH  is  weighted  more  by  the 
oblate  half  melted  snowflakes  while  ZV  is  reflec¬ 
ting  more  the  signal  from  the  raindrops .  We 
suggest  that  the  high  frequency  fluctuations  in 
ZH  are  caused  by  the  rapid  oscillatory  and  spin¬ 
ning  motion  of  the  melting  snowflakes  (Mitra  et 
al  ,  1988)  modulating  the  values  of  the  scattering 
cross  sections;  the  slower  variation  in  ZV  are 
caused  by  the  changes  in  phase  of  the  scattered 
signals  as  the  particles  move  in  space.  The  9ame 
features  are  present  when  several  time  series  are 
averaged  together.  Ten  series  and  a  dwell  time 
of  2.5-seconds  are  usually  needed  to  obtain  relia¬ 
ble  spectra;  this  leads  to  very  slow  scan  rates 
for  a  quarter  degree  beamwidth  antenna. 

The  values  of  the  correlation  between  the 
co-polar  and  the  cross-polar  returns  (p  vVH,H) 
are  generally  not  significantly  different  from 
zero  and  do  not  appear  to  contain  useful  infor¬ 
mation.  However  when  the  antenna  limit  is  reached 
the  correlation  can  rise  to  0.3  to  0.5,  presuma¬ 
bly  caused  by  the  signal  from  the  main  channel 
leaking  through  the  antenna  into  the  cross-polar 
channel.  The  spectra  for  the  two  tin*.  series  in 
Figure  l  in  heavy  rain  are  plotted  in  Figure  6  , 
This  illustrates  the  widespread  finding  that  the 
ZVH  spectrum  has  much  higher  components  than  the 
ZH  spectrum,  and  that  the  low  frequency  maximum 
in  the  ZH  spectrum  is  absent  in  the  cross-polar 
spectrum.  The  ZVH  spectrum  reflects  the  rare  of 
oscillations,  spinning  and  tumbling  of  the  parti¬ 
cles.  The  high  frequency  tail  in  the  ZVH  spectrum 
above  80Hz  is  probably  caused  by  system  noise; 
tor  this  spectrum  the  ZH  signal  to  noise  ratio 
was  ;6dB  and  chat  for  the  cross-polar  channel  was 
25dB. 

5.  CONCLUSIONS 

3ased  on  the  wo- <  presented  in  this  paper 
and  CAI  we  suggest: 

(a)  Low  values  of  the  co-polar  correlation  can 
be  used  for  rapid  identification  of  the 
aright  band. 

(b)  Low  values  of  the  co-polar  correlation  in 
vigorous  convective  clouds  could  indicate 
the  presence  of  Large  hail. 

'  c )  LDR  measurements  in  rain  are  consistent 
with  a  Gaussian  distribution  of  canting 
angles  with  a  standard  deviation  of  5°. 
fd)  LDR  values  in  the  bright  band  car  be  used 
for  distinguishing  between  snow  and  graupel. 
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Figure  5.  Spectra  of  the  co-polar  time  series  in 
the  melting  layer. 

Power  Spectrum  on  29  Hay  1980  at  160418  UT 
Ray  -  24  Gate  -  62  Range  -  9.3km 
Az  -  290.0deg  £1  -  6.7deg 


Figure  6.  Cross  and  co-polar  spectra  for  the 
time  series  in  Figure  1. 
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BRIGHT  BAND  ERRORS  IN  RAINFALL  MEASUREMENT:  IDENTIFICATION 
AND  CORRECTION  USING  LINEARLY  POLARIZED  RADAR  RETURNS 
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ABSTRACT 

A  principle  source  of  error  in  rainfall  rates  derived  from  the  radar 
reflectivity  (Z)  is  caused  by  the  enhanced  radar  return  due  to  melting 
snowflakes  in  the  bright  band.  We  present  observations  made  by  the  narrow 
beamwidth  S-band  Chilbolton  radar  which  involve  transmission  of 
horizontally  and  vertically  polarized  radiation  and  reception  of  the 
co-polar  and  cross-polar  return  signals.  The  linear  depolarization  ratio 
(LDR)  is  defined  as  the  ratio  of  the  cross-polar  to  the  co-polar  return. 
The  high  values  of  Z  in  the  bright  band  are  accompanied  by  values  of  LDR 
above  -18dB;  in  echoes  where  no  bright  band  is  present  the  values  of  LDR 
are  everywhere  below  -20dB. 


1.  INTRODUCTION 

The  conventional  radar  reflectivity,  Z,  is  proportional  to  NI?  ,  where  N  is 
the  concentration  of  particles  of  diameter  D,  summed  over  all  sizes.  Z  is 

usually  expressed  in  units  relative  to  the  signal  from  a  1mm  raindrop  per 
cubic  meter,  even  though  from  Z  alone  it  is  not  possible  to  distinguish 
rain  from  ice  and  the  reflectivity  of  a  raindrop  is  about  7dB  higher  than 
the  equivalent  mass  of  ice.  Neither  can  Z  be  used  to  differentiate  between 
the  various  forms  of  frozen  hydrometeors  (snow,  hail,  hailstones  etc),  or 
to  measure  the  sizes  and  concentrations  of  raindrops.  Z  is  usually  con¬ 
verted  into  a  rain  rate,  R,  using  an  empirical  relation  of  the  form: 

Z  =  284  R  l'47  (1) 

which,  for  a  given  R,  is  equivalent  to  assuming  a  constant  size 

distribution  of  raindrops.  Some  of  the  problems  in  estimating  rainfall 

using  Equation  1  are  demonstrated  in  Figure  1,  which  is  a  vertical 
section  of  Z  through  stratiform  rainfall  observed  by  the  quarter  degree 
beamwidth  Chilbolton  radar. 

The  most  notable  feature  in  Figure  1  is  the  layer  of  enhanced  reflectivity 
or  'bright  band’  at  about  1.1km  altitude.  At  a  range  of  45km  the  value  of 
Z  reaches  43dBZ  leading  (via  Equation  1)  to  an  estimated  rain  rate  of  about 
18mm/hr.  Values  below  150m  are  affected  by  obscuration  of  the  radar  beam, 
but  from  200m  to  800m,  in  the  rain,  the  value  of  Z  is  only  28dBZ, 
equivalent  to  a  rainfall  rate  of  only  1.7  mm/hr.  The  bright  band  leads  to 
an  overestimate  of  R  by  a  factor  of  ten. 

The  enhanced  return  in  the  bright  band  is  caused  by  large,  low  density  wet 
snowflakes  which  reflect  microwaves  as  if  they  were  giant  raindrops.  Dry 
snowflakes  have  a  lower  return  because  of  their  low  dielectric  constant; 
below  the  bright  band  Z  falls  due  to  two  factors;  when  the  snowflakes  melt 
completely  they  collapse  to  smaller  raindrops  and,  secondly,  the  raindrop 
concentration  falls  as  the  terminal  velocity  increases. 

Other  problems  are  evident  from  Figure  1.  The  Z  values  at  altitudes 
greater  than  2.5km  are  lower  than  in  the  rain,  much  of  the  growth  and 
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aggregation  of  the  ice  occurring  at  lower  levels.  Consequently,  the  higher 
elevation  scans  of  a  one  degree  beamwidth  radar  will  underestimate  rainfall 
at  larger  distances  (1),  while  sampling  of  the  rain  alone  with  the  lower 
elevation  scans  without  ground  clutter  and  /or  obscuration  is  difficult, 
because,  in  the  UK,  the  melting  layer  is  usually  below  2km. 

Ground  based  rain  gauges  can  provide  localised  real-time  information  for 
correcting  bright  band  errors  (2).  Most  radar  networks  scan  in  PPI  mode  to 
obtain  a  complete  spatial  coverage,  and  in  truly  stratiform  rain  the  bright 
band  should  be  recognisable  as  a  concentric  ring  of  enhanced  reflectivity 
centred  on  the  radar.  Smith  (3)  has  suggested  an  automatic  means  of 
identifying  the  bright  band  by  comparing  the  range  of  the  maximum  values  of 
Z  at  a  particular  azimuth  for  two  different  beam  elevations.  However,  in 
practice  the  height  of  the  bright  band  can  change,  the  precipitation  is 
never  stratiform  and,  in  the  UK  at  least,  quite  vigorous  showers  often  have 
bright  bands.  In  this  paper  we  demonstrate  a  means  of  uniquely  identifying 
the  bright  band  by  analyzing  the  cross-polar  radar  return. 

2.  THE  CHILBOLTON  POLARIZATION  RADAR 

The  Chilbolton  radar  operates  at  S-band  (10cm)  and,  with  a  25m  dish,  is  the 
largest  steerable  meteorological  radar  in  the  world,  having  a  beamwidth  of 
only  a  quarter  of  a  degree.  Earlier  reports  (4,5,6)  have  considered  the 
implementation  and  interpretation  of  differential  reflectivity  (ZDR),  we 
now  analyze  observations  made  in  1988  of  a  new  parameter,  the  linear 
depolarization  ratio  ( LDR) . 

The  differential  reflectivity  (ZDR)  provides  an  estimate  of  mean 
hydrometeor  shape.  It  is  defined  as 

ZDR  =  10  log( ZH/ZV)  (2) 

where  ZH  and  ZV  are  the  radar  reflectivities  measured  at  horizontal  and 
vertical  polarizations  respectively.  For  small  raindrops  or  tumbling  ice 
particles,  ZH  and  ZV  are  equal  and  ZDR  is  zero.  The  theoretical  values  of 
ZDR  for  oblate  particles  with  their  minor  axes  aligned  in  the  vertical  are 
plotted  in  Figure  2;  ZDR  increases  with  greater  oblateness  and  higher 
dielectric  constant.  In  heavier  rain  ZDR  is  positive  and  reflects  the  mean 
shape  (and  hence  the  size)  of  the  raindrops.  The  ZDR  of  ice  is  more 
complex  (6).  Because  of  the  low  dielectric  constant,  dry  snowflakes  have  a 
ZDR  close  to  zero,  but  wet  snowflakes  can  have  high  positive  values. 
Graupel  tends  to  tumble  and  so  be  associated  with  a  zero  ZDR  value. 

The  linear  depolarization  ratio,  LDR,  is  a  measure  of  the  hydrometeor  fall 
mode  and  appears  to  be  an  excellent  indicator  of  wet  ice.  It  is  defined  as 

LDR  =  10  log(ZVH/ZH)  (3) 

where  ZVH  is  the  (horizontal)  cross-polar  return  from  a  vertically 

polarised  transmitted  pulse,  and  ZH  (as  in  Equation  2)  is  the  co-polar 
(horizontal)  return  for  horizontally  polarized  transmission.  A  cross 
polar  return  occurs  only  when  oblate  hydrometeors  fall  with  their  major  or 
minor  axis  at  an  angle  to  the  vertical.  Computations  of  LDR  for  tumbling 
oblate  spheroids  (Figure  3)  and  are  consistent  with  the  Chilbolton 
observations  (7).  Snowflakes  have  such  a  low  dielectric  constant  that  even 
if  they  are  very  oblate  their  LDR  is  below  the  antenna  limit  of  -32dB, 
oblate  dry  hail  or  graupel  could  have  a  value  up  to  -20dB  if  the  axial 
ratio  were  as  low  as  0.5,  but  LDR  values  above  -20dB  can  only  realistically 
occur  for  wet  tumbling  ice  particles.  Such  high  values  are  restricted  to 
the  bright  band.  Raindrops  give  rise  to  a  very  low  cross-polar  return. 


Figure  2.  ZDR  values  as  a  function  of  axial  ratio  for  various  precipitation 
particles,  assuming  the  major  axis  is  horizontal. 


Figure  3.  LDR  values  for  randomly  tumbling  particles  as  a  function  of 
axial  ratio. 


3.  COMPARISON  OF  ECHOES  WITH  AND  WITHOUT  BRIGHT-BANDS 

An  RHI  through  widespread  stratiform  precipitation  is  displayed  in  Figure 

4,  where  the  reflectivity  does  not  exceed  30dBZ  in  the  rain  but  reaches 
40dBZ  at  2km  altitude  in  the  bright  band.  At  this  height  the  oblate 
melting  snowflakes  give  a  clearly  visible  bright  band  in  ZDR  with  values 
reaching  2dB.  However  automatic  recognition  of  the  ZDR  bright  band  can  be 
difficult.  ZDR  values  in  the  rain  at  10-20km  range  reach  0.5dB,  and  in 
heavier  rain  can  be  much  higher.  We  also  note  in  Figure  4  the  positive 
ZDR  values  above  3km  altitude;  this  low  Z  region  presumably  containing 
aligned  high  density  ice  crystals. 

It  is  much  easier  to  identify  the  Z  bright  band  from  the  LDR  data.  The 
maximum  values  of  Z  coincide  with  the  peak  LDR  of  about  -15dB,  which  is 
consistent  with  wet  tumbling  snowflakes  having  an  axial  ratio  of  about  0.5 
(Figure  3).  In  contrast  LDR  values  in  the  rain  are  near  the  antenna  limit 
of  -32dB,  and  reach  about  -27dB  in  the  low  Z  ice  region  above  3km  where  the 
ZDR  indicated  high  density  crystals.  Figure  4  also  shows  that  ground 
clutter  results  in  LDR  values  above  -lOdB  near  to  the  ground.  Because  LDR 
involves  measuring  the  low-power  cross-polar  return  it  is  much  more 
susceptible  to  ground  clutter  than  is  Z  or  ZDR. 

Figure  5  illustrates  an  example  of  a  heavy  shower  with  a  horizontal  extent 
of  about  25km.  The  Z  values  in  the  rain  are  higher  than  in  Figure  4,  but 
the  enhancement  of  Z  in  the  bright  band  is  over  lOdB.  Again  the  presence 
of  the  bright  band  can  be  most  easily  identified  by  the  values  of  LDR  which 
exceed  -16dB. 

A  vigorous  shower  with  no  bright  band  is  depicted  in  Figure  6  and  it  is 
clear  that  the  polarization  parameters  have  a  quite  different  character: 
the  LDR  values  are  much  lower  than  for  the  bright  band  case  in  Figures  4 
and  5.  We  believe  that  the  data  in  Figure  6  indicate  the  presence  of 
graupel.  The  dry  tumbling  graupel  gives  negligible  LDR,  but  melting  occurs 
at  about  2km  altitude  and  LDR  rises  to  about  -25dB.  This  weak  "LDR  graupel 
bright  band"  is  consistent  (Figure  3)  with  tumbling  wet  ice  with  an  axial 
ratio  of  0.8.  In  the  heavy  rain  the  LDR  is  just  detectable  and  is 
explicable  in  terms  of  a  canting  angle  of  about  5  degrees.  Values  of  ZDR 
are  low  for  the  tumbling  dry  ice,  but  rise  monotonically  as  the  graupel 
melts  and  assumes  the  equilibrium  shape  of  the  large  raindrops.  This 
vertical  profile  in  ZDR  should  be  contrasted  with  the  bright  band  case  in 
Figures  4  and  5,  where  a  maximum  in  ZDR  is  caused  by  the  low  density  oblate 
wet  snowflakes,  which  subsequently,  on  complete  melting,  collapse  to  more 
spherical  raindrops. 

It  should  be  emphasised  that,  in  the  UK  at  least,  the  presence  of  a  bright 
band  is  not  restricted  to  stratiform  clouds  with  low  Z.  In  some  showers  Z 
values  can  reach  50dBZ  in  the  bright  band,  while  others,  with  no  bright 
band,  have  lower  peak  values  of  Z. 

4.  LDR  STATISTICS 

In  order  to  test  our  hypothesis  that  the  peak  value  of  LDR  in  a  vertical 
profile  is  related  to  the  increased  reflectivity  in  the  bright  band,  the 
results  from  scans  on  11  different  days  in  1988  are  summarised  in  Figure  7. 
The  altitude  of  the  maximum  value  in  LDR  (which  is  found  in  all  types  of 
cloud)  is  used  to  fix  the  melting  level.  The  enhancement  of  Z  (delta  Z)  is 
then  estimated  by  comparing  the  Z  at  the  melting  level  with  the  Z  in  the 
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Figure  4.  A  typical  RHI  scan  through  stratiform  cloud  on  29  November  1988 
with  a  bright  band  for  Z,  ZDR  and  LDR. 
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figure  5.  An  RHI  scan  of  Z,  ZDR  and  LDR  through  a  vigorous  shower  with  a 
bright  band.  29  May  1988. 
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Figure  6.  A  typical  RHI  scan  through  a  convective  cloud  on  13  July  1988 
with  no  bright  band  for  Z,  ZDR  and  LDR. 
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rain  500m  below.  Histograms  of  the  enhancement  of  the  reflectivity  in  the 
bright  band  are  plotted  for  each  4dB  increment  in  LDR.  For  most  vertical 
profiles  the  peak  LDR  values  are  in  the  range  -14  to  -18dB  and  in  these 
cases  the  enhancement  of  Z  is,  on  average,  about  lOdB.  Less  common,  in 
this  UK  sample,  are  the  peak  values  of  LDR  in  the  range  -18  to  -22dB  and 
-22  to  -26dB  where  the  Z  enhancement  is  essentially  zero  and  no  bright  band 
is  present. 

It  should  be  stressed  that  in  1988  there  were  no  observed  cases  of  deep 
vigorous  convection.  Measurements  of  LDR  at  3cm  (8)  suggest  that  hail  in 
wet  growth  can  give  high  LDR  values,  although  the  measurements  at  this 
wavelength  are  affected  by  propagation  problems.  For  the  observations 
discussed  in  this  paper  the  depth  of  the  bright  band  is  greater  than  the 
beamwidth  of  the  Chilbolton  radar.  In  future  we  will  analyze  the  effect 
on  the  LDR  measurements  for  a  one  degree  beamwidth  radar  which  is  only 
partially  filled  by  the  bright  band. 

5.  DISCUSSION  AND  CONCLUSION 

Several  factors  need  to  be  considered  if  the  LDR  technique  is  to  be  implem¬ 
ented  on  conventional  C-band  radars.  The  LDR  measurement  has  the  advantage 
that  no  fast  switching  of  the  transmitted  signal  is  needed.  However,  the 
cross-polar  signal  power  is  very  low,  and  a  signal  to  noise  ratio  of  more 
than  20dB  is  required  if  an  LDR  of  down  to  -20dB  is  to  be  detected.  We 
also  note  that  the  low  power  cross-polar  signal  is  much  more  sensitive  than 
Z  and  ZDR  to  ground  clutter  contamination.  Finally,  and  most  importantly, 
the  LDR  signal  is  affected  by  propagation  and  attenuation  problems  at 
shorter  wavelengths.  Correction  algorithms  may  be  possible  at  C-band. 

These  observations  suggest  that  high  values  of  the  linear  depolarization 
ratio  are  associated  with  the  presence  of  melting  snowflakes  and  can  be 
used  to  identify  the  bright  band.  This  parameter  may  also  be  of  use  in 
identifying  anomalous  propagation  and  ground  clutter,  both  of  which  should 
have  high  values  of  LDR. 
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Summary 

A  principle  source  of  error  in  rainfall  rates  derived  from  the 
radar  reflectivity  (Z)  is  caused  by  the  enhanced  radar  return  due  to 
melting  snowflakes  in  the  bright  band.  We  present  observations  of 
two  S-band  polarization  radar  parameters  which  enable  the  bright  band 
to  be  simply  identified.  The  technique  involves  transmission  of 
horizontally  and  vertically  polarized  radiation  and  reception  of  the 
co-polar  and  cross-polar  return  signals.  The  linear  depolarization 
ratio  ( LDR)  is  defined  as  the  ratio  of  the  cross-polar  to  the  co- 
polar  return  and  p(H,V)  as  the  correlation  of  the  time  series  of  the 
horizontal  and  vertical  co-polar  return.  The  high  values  of  Z  in  the 
bright  band  are  accompanied  by  values  of  LDR  above  -18dB  and  p(H,V) 
below  0.8.  In  echoes  where  no  bright  band  is  present  the  values  of 
LDR  are  everywhere  below  -20dB  and  the  correlation  is  always  close  to 
unity.  We  discuss  potential  problems  in  implementing  these 
techniques  for  C-band  radars  with  one  degree  beamwidths,  and  also 
consider  how  they  could  be  used  to  identify  spurious  echoes  from 
ground  clutter  and  anomalous  propagation. 


1.  INTRODUCTION 

A  ma]or  source  of  error  in  deriving  the  rainfall  rate  (R)  from  the  radar 
reflectivity  (Z)  arises  from  the  enhanced  radar  return  occurring  in  the 
melting  layer  or  ’bright  band'.  The  value  of  Z  typically  increases  by 
lOdB  when  low  density  snowflakes  become  wet  and  scatter  microwaves  as  if 
they  were  giant  raindrops.  This  error  in  Z  in  an  empirical  Z-R 
relationship  would  lead  to  a  fivefold  overestimate  of  the  rainfall. 

Most  radar  networks  scan  in  PPI  mode  to  obtain  a  complete  spatial 
coverage,  and  in  truly  stratiform  rain  the  bright  band  should  be 
recognisable  as  a  concentric  ring  of  enhanced  reflectivity  centred  on  the 
radar.  Smith  (1)  has  suggested  an  automatic  means  of  identifying  the 
bright  band  by  comparing  the  range  of  the  maximum  values  of  Z  at  a 
particular  azimuth  for  two  different  beam  elevations.  However,  in  practice 
the  height  of  the  bright  band  can  change,  the  precipitation  is  never 
stratiform  and,  in  the  UK  at  least,  quite  vigorous  showers  often  have 
bright  bands.  Ground  based  rain  gauges  can  provide  localised  real-time 
information  for  correcting  bright  band  errors  (2).  In  this  paper  we 
demonstrate  a  means  of  uniquely  identifying  the  bright  band  using  two  new 
polarization  parameters.  We  shall  also  consider  how  such  techniques  could 
be  implemented  on  C-band  radars  with  smaller  antennas. 


2.  THE  CHILBOLTON  POLARIZATION  RADAR 

The  Chilbolton  radar  operates  at  S-band  and,  with  a  25m  dish,  is  the 
largest  steerable  meteorological  radar  in  the  world,  having  a  beamwidth  of 
only  a  quarter  of  a  degree.  Earlier  reports  (3,4,5)  have  considered  the 
implementation  and  interpretation  of  differential  reflectivity  (ZDR),  but 
here  we  discuss  observations  made  in  1988  of  two  new  parameters  (LDR  and 
p ( H ,- V ) )  which  are  described  below. 

2.1  DIFFERENTIAL  REFLECTIVITY 

The  differential  reflectivity  (ZDR)  provides  an  estimate  of  mean 
hydrometeor  shape.  It  is  defined  as 

ZDR  =  10  log( ZH/ZV)  (1) 

where  ZH  and  ZV  are  the  radar  reflectivities  measured  at  horizontal  and 
vertical  polarizations  respectively.  For  small  raindrops  or  tumbling  ice 
particles,  ZH  and  ZV  are  equal  and  ZDR  is  zero.  Positive  values  of  ZDR 
occur  for  oblate  particles  of  high  dielectric  constant  when  ZH  exceeds  ZV. 
In  heavier  rain  ZDR  is  positive  and  reflects  the  mean  shape  (and  hence  the 
size)  of  the  raindrops.  The  ZDR  of  ice  is  more  complex  (5).  Because  of 
the  low  dielectric  constant,  dry  snowflakes  have  a  ZDR  close  to  zero,  but 
wet  snowflakes  can  have  high  positive  values.  Graupel  tends  to  tumble  and 
so  be  associated  with  a  zero  value. 

2.2  LINEAR  DEPOLARIZATION  RATIO 

The  linear  depolarization  ratio,  LDR,  is  a  measure  of  the  hydrometecr 
fall  mode  and  appears  to  be  an  excellent  indicator  of  wet  ice.  LDR  is 
defined  as: 

LDR  =  10  log( ZVH/ZH)  (2) 

where  ZVH  is  the  (horizontal)  cross-polar  return  from  a  vertically 
polarised  transmitted  pulse,  and  ZH  (as  in  Equation  1)  is  the  co-polar 
(horizontal)  return  for  horizontally  polarised  transmission. 

A  cross  polar  return  occurs  only  when  oblate  hydromete.ors  fall  with  their 
major  or  minor  axis  at  an  angle  to  the  vertical.  Computations  of  LDR  for 
tumbling  oblate  spheroids  are  plotted  in  Figure  1  and  are  found  to  be 
consistent  with  the  Chilbolton  observations  (6).  Snowflakes  have  such  a 
low  dielectric  constant  that  even  if  they  are  very  oblate  their  LDR  is 
below  the  antenna  limit  of  -3ZdB,  oblate  dry  hail  or  graupel  could  have  a 
value  up  to  -2QdB  if  the  axial  ratio  were  as  low  as  0.5,  but  LDR  values 
above  -20dB  can  only  realistically  occur  for  wet  tumbling  ice  particles. 
Such  high  values  are  restricted  to  the  bright  band.  Raindrops  give  rise  to 
a  very  low  cross-polar  return. 

2.3  CO-POLAR  CROSS  CORRELATION 

The  estimates  of  ZH  and  ZV  in  equation  (1)  are  made  from  the  true  linear 
average  (over  210msec)  of  the  return  at  one  75m  gate  from  64  successive 
pulse  pairs  transmitted  with  alternate  horizontal  and  vertical 
polarization.  The  standard  error  in  ZDR  is  reduced  to  O.ldB  by  spatial 
averaging  over  four  adjacent  75m  gates.  The  co-polar  cross  correlation 
( p( H, V) )  is  the  correlation  of  these  two  time  series  in  ZH  and  ZV. 

Observations  (7)  show  that  the  correlation  is  generally  close  to  unity 
in  rain  and  dry  ice,  with  low  values  being  confined  to  the  bright  band. 
Low  values  of  correlation  are  thought  to  indicate  that  a  variety  of 
hydrometeor  shapes  is  present.  We  believe  that  the  low  values  in  the  bright 


band  are  caused  by  the  coexistence  of  oblate  half-melted  snowflakes  with 
nearly  spherical  raindrops. 

3.  COMPARISON  OF  ECHOES  WITH  AND  WITHOUT  BRIGHT-BANDS 

An  RHI  through  stratiform  precipitation  is  displayed  in  Figure  2;  where  the 
reflectivity  does  not  exceed  30dBZ  in  the  rain  but  reaches  40dBZ  at  2km 
altitude  in  the  bright  band.  At  this  height  the  oblate  melting  snowflakes 
give  a  clearly  visible  bright  band  in  ZDR  with  values  reaching  2dB.  How¬ 
ever  automatic  recognition  of  the  ZDR  bright  band  can  be  difficult.  ZDR 
values  in  the  rain  at  10-20km  range  reach  0.5dB,  and  in  heavier  rain  can  be 
much  higher.  We  also  note  in  Figure  2  the  positive  ZDR  values  above  3km 
altitude;  this  low  Z  region  presumably  containing  aligned  high  density  ice 
crystals. 

It  is  much  easier  to  identify  the  Z  bright  band  from  the  LDR  data.  The 
maximum  values  of  Z  coincide  with  the  peak  LDR  of  about  -15dB,  which  is 
consistent  with  wet  tumbling  snowflakes  having  an  axial  ratio  of  about  0.5 
(Figure  1).  In  contrast  LDR  values  in  the  rain  are  near  the  antenna  limit 
of  -32dB,  and  reach  about  -27dB  in  the  low  Z  ice  region  above  3km  where  the 
ZDR  indicated  high  density  crystals.  Figure  2  also  shows  that  ground 
clutter  results  in  LDR  values  above  -lOdB  near  to  the  ground.  Because  LDR 
involves  measuring  the  low-power  cross-polar  return  it  is  much  more 
susceptible  to  ground  clutter  than  is  Z  or  ZDR.  The  bright  band  can  also 
be  identified  via  the  correlation  parameter.  Although  data  is  limited  to  a 
5km  range  window  values  of  correlation  below  90%  are  restricted  to  the 
bright  band.- 


Figure  1.  LDR  values  for 
axial  ratio. 


randomly  tumbling  particles  as  a  function  of 


A  vigorous  shower  with  no  bright  band  is  depicted  in  Figure  3  and  it  is 
clear  that  the  polarization  parameters  have  a  quite  different  character: 
the  correlation  is  high  everywhere  while  the  LDR  values  are  much  lower  than 
for  the  bright  band  case  in  Figure  2.  The  different  vertical  profiles 
through  these  two  clouds  are  plotted  in  Figures  4  and  5.  We  believe  that 
the  vigorous  shower  with  no  bright  band  in  Figures  3  and  5  contains 
graupel.  The  dry  tumbling  graupel  gives  negligible  LDR,  but  melting  occurs 
at  about  2km  altitude  and  LDR  rises  to  about  -25dB.  This  weak  "LDR  graupel 
bright  band"  is  consistent  (Figure  1)  with  tumbling  wet  ice  with  an  axial 
ratio  of  0.8.  In  the  heavy  rain  the  LDR  is  just  detectable  and  is 
explicable  in  terms  of  a  canting  angle  of  about  5  degrees.  Values  of  ZDR 
are  low  for  the  tumbling  dry  ice,  but  rise  monotonically  as  the  graupel 
melts  and  assumes  the  equilibrium  shape  of  the  large  raindrops.  These 
profiles  should  be  contrasted  with  the  bright  band  case  in  Figure  4,  where 
a  maximum  in  ZDR  is  caused  by  the  low  density  wet  snowflakes.  In  the 
graupel  case  the  correlation  is  everywhere  above  901  as  there  is  no  great 
variety  of  shapes  present. 

Figures  6  and  7  display  how  the  value  of  LDR  is  related  to  the  enhanced 
value  of  Z  in  the  bright  band.  The  altitude  of  the  maximum  value  in  LDR 
(which  is  present  in  all  types  of  cloud)  is  used  to  fix  the  melting  level. 
The  enhancement  of  Z  (delta  Z)  is  then  estimated  by  comparing  the  Z  at  the 
melting  level  with  the  Z  in  the  rain  500m  below.  In  Figure  6,  for  the 
bright  band  case  in  Figures  2  and  4,  the  delta  Z  enhancement  in  the  bright 
band  for  ranges  out  to  60km  is  about  lOdB,  and  is  accompanied  by  an  LDR 
value  of  about  -15dB.  In  Figure  7  these  parameters  are  plotted  for  both 
the  graupel  cloud  in  Figures  3  and  5  at  a  range  of  10-20km,  and  for  a 
second  shower  beyond  40km  which  does  have  a  bright  band.  In  the  graupel 
clouds  the  mean  delta  Z  enhancement  is  close  to  zero  and  LDR  values  are  in 
the  range  -20  to  -25dB,  while  for  the  more  distant  cloud  the  bright  band 
increase  in  Z  of  about  lOdB  is  associated  with  higher  LDR  values  of  -15dB. 

It  should  be  emphasised  that,  in  the  UK  at  least,  the  presence  of  a  bright 
band  is  not  restricted  to  stratiform  clouds  with  low  Z.  In  some  showers  Z 
values  can  reach  50dBZ  in  the  bright  band,  while  others,  with  no  bright 
band,  have  lower  peak  values  of  Z. 

4.  LDR  STATISTICS 

In  order  to  test  our  hypothesis  that  the  value  of  LDR  is  related  to  the 
increased  reflectivity  in  the  bright  band,  the  results  from  scans  on  11 
different  days  in  1988  are  summarised  in  Figure  8.  Histograms  of  the 
enhancement  of  the  reflectivity  in  the  bright  band  are  plotted  for  each  4dB 
increment  in  LDR.  For  most  vertical  profiles  the  peak  LDR  values  are  in 
the  range  -14  to  -18dB  and  in  these  cases  the  enhancement  of  Z  is,  on 
average,  about  lOdB.  Less  common,  in  this  UK  sample,  are  the  peak  values 
of  LDR  in  the  range  -18  to  -22dB  and  -22  to  -26dB  where  the  Z  enhancement 
is  essentially  zero  and  no  bright  band  is  present. 

It  should  be  stressed  that  in  1988  there  were  no  cases  of  very  deep 
vigorous  convection.  Measurements  of  LDR  at  3cm  (8)  suggest  that  hail  in 
wet  growth  can  give  high  LDR  values,  although  the  measurements  at  this 
wavelength  are  affected  by  propagation  problems.  We  hope  to  examine  such 
cases  in  the  future.  For  the  observations  discussed  in  this  paper  the 
depth  of  the  bright  band  is  greater  than  the  beamwidth  of  the  Chilbolton 
radar.  In  future  we  will  analyze  the  effect  on  LDR  and  correlation 
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3.  A  typical  RHI  scan  through  a  convective  cloud 
on  13  July  1988  with  no  bright  band  for  Z, 
ZDR,  LDR  and  p(H, V) . 
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Figure  6.  The  enhancement  of  (delta  Z)  of  Z  in  the  Figure  7.  The  enchancement  of  (delta  Z)  of  Z  in  the 
bright  band  ( + )  compared  to  the  rain  below  (Q  )  as  a  bright  band  (♦)  compared  to  the  rain  below  (  Q  |  as  a 
function  of  LDR  in  the  bright  band  (  O  )  for  the  scan  in  function  of  LDR  in  the  bright  band  ( O )  for  the  scan  in 
Figure  2.  Figure  3. 
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Figure  8.  Histograms  of  the  enhancement  of  the  reflectivity  in  the  bright 
band  (delta  Z)  as  a  function  of  LDR  in  the  bright  band. 


measurements  for  a  one  degree  beamwidth  radar  which  is  only  partially 
filled  by  the  bright  band. 

5.  TECHNICAL  ASPECTS  OF  IMPLEMENTATION 

The  Chilbolton  radar  is  a  narrow  beam  research  radar  which  can  make 
measurements  of  unrivalled  polarization  purity.  Several  factors  need  to  be 
considered  if  the  techniques  are  to  be  implemented  on  conventional  C-band 
radars. 

The  LDR  measurement  has  the  advantage  that  no  fast  switching  of  the 
transmitted  signal  is  needed.  However,  the  cross-polar  signal  power  is 
very  low,  and  a  signal  to  noise  ratio  of  more  than  20dB  is  required  if  an 
LDR  of  -20dB  is  to  be  detected.  We  also  note  that  the  low  power 
cross-polar  signal  is  much  more  sensitive  than  Z  and  ZDR  to  ground  clutter 
contamination.  Finally,  and  most  importantly,  the  LDR  signal  is  affected 
by  propagation  and  attenuation  problems  at  shorter  wavelengths.  Correction 
algorithms  may  be  possible  at  C-band. 

If  the  co-polar  cross  correlation  is  to  be  measured  then  rapid 
pule-to-pulse  switching  of  the  polarization  of  the  transmitted  signal  is 
required.  However,  the  correlation  is  unaffected  by  attenuation  and 
propagation,  provided  there  is  sufficient  signal  level.  A  signal  to  noise 
ratio  of  20dB  is  needed  if  the  correlation  is  to  be  measured  to  1%.  At 
shorter  wavelengths  the  correlation  measurement  will  be  degraded  as  the 
time  between  transmitted  pulses  approaches  the  decorrelation  time  of  the 
return  signal.  ZDR  measurements  generally  need  long  dwell  times,  but  it  is 
possible  to  estimate  the  correlation  from  shorter  time  series  and  so  scan 
more  rapidly. 

6.  CONCLUSION 

These  observations  suggest  that  the  linear  depolarization  ratio  and/or  the 
co-polar  cross  correlation  can  be  used  to  identify  the  bright  band.  These 
parameters  may  also  be  of  use  in  identifying  anomalous  propagation  and 
ground  clutter,  both  of  which  should  have  high  values  of  LDR  and  low 

correlations. 
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ABSTRACT 

The  passage  of  aircraft  through  clouds  can  trigger  lightning  which  would  not  occur 
naturally.  We  describe  such  a  triggered  event  which  occurred  while  the  cloud  was 
simultaneously  being  scanned  by  the  Chllbolton  polarization  radar.  This  radar  measures  new 
radar  parameters:  the  differential  reflectivity  (ZDR)  which  provides  an  estimate  of  the  mean 
shape  of  the  precipitation  particles,  and  the  linear  depolarization  ratio  (LDR)  which 
reflects  their  fall  mode.  Using  this  information  it  is  possible  to  distinguish  between 
clouds  which  Contain  snowflakes  and  those  where  graupel  or  small  hail  pellets  are  present. 
Measurements  of  particle  type  made  with  the  aircraft  confirm  the  radar  inferences  that  the 
triggered  lightning  occurred  when  the  aircraft  was  traversing  a  region  of  hail  pellets. 
This  conforms  with  our  knowledge  of  charge  generation  mechanisms  within  convective  clouds, 
which  predict  that  appreciable  electric  fields  should  be  restricted  to  those  regions  of 
clouds  containing  graupel.  It  appears  that  these  radar  parameters  provide  a  means  of 
remotely  locating  clouds  which  pose  a  threat  of  triggered  lightning  to  aircraft. 


1.  INTRODUCTION 

There  is  evidence  that  many  lightning 
strikes  to  aircraft  are  triggered  by  the 
presence  of  the  aircraft  itself.  The 
existing  field  in  the  cloud  is  insufficient 
to  Initiate  lightning  naturally  but  the 
presence  of  a  large  conductor  raises  the 
field  above  the  critical  level  (1).  Once 
lightning  occurs,  many  techniques  are 
available  to  locate  its  precise  position 
and  to  study  its  evolution  in  time  and 
space.  However,  detection  of  electric 
fields  within  clouds  which  are  high,  but 
insufficient  to  initiate  natural  lightning, 
is  much  more  problematic.  The  electro¬ 
static  field  of  the  electric  dipole 
structure  of  the  cloud  falls  off  with  the 
inverse  cube  of  the  distance  from  the 
cloud,  so  that,  apart  from  very  restricted 
areas.  it  is  impractical  to  install  a 
network  of  ground  based  field  sensors. 

In  Section  2  of  this  paper  we  describe  a 
new  polarization  radar  technique  which  is 
able  to  identify  the  different  types  of 
precipitation  particle  present  within 
clouds.  Examples  of  radar  observations  of 
different  clouds.  including  one  which 
triggered  a  lightning  flash  to  the 
aircraft,  are  presented  in  Section  3.  The 
aircraft  measurements  of  particle  type 
during  the  penetration  when  the  triggering 
occurred  are  compared  with  aircraft 
inferen  es  in  Section  4,  and  the 
implications  of  these  observations  are 
discussed  in  Section  S. 

2.  THE  CHILBOLTON  POLARIZATION  RADAR 

The  Chllbolton  radar  (2),  situated  in 
Hampshire,  operates  at  S-band  (10cm)  and, 
with  a  25m  dish,  is  the  largest  steerable 
radar  in  the  world,  having  a  beamwidth  of 
only  a  quarter  of  a  degree.  It  can  make 
measurements  of  the  conventional  radar 
reflectivity  factor,  Z,  with  unrivalled 
resolution,  but  is  also  able  to  make 
additional  polarization  measurements  which 
help  identification  of  the  precipitation 
particles. 


The  conventional  radar  reflectivity,  Z,  is 
proportional  to  ND,  where  N  is  the  con¬ 
centration  of  particles  of  diameter  D, 
summed  over  all  sizes.  Z  is  usually  ex¬ 
pressed  in  units  of  dB  relative  to  the 
signal  from  a  1mm  raindrop  per  cubic  meter. 
The  reflectivity  of  a  raindrop  is  7dB 
higher  than  the  equivalent  mass  of  ice,  but 
from  Z  alone  it  is  not  possible  to 
distinguish  rain  from  ice.  Neither  can  Z 
be  used  to  differentiate  between  the 
various  forms  of  frozen  hydrometeors  (snow, 
hail,  hailstones  etc),  or  to  measure  the 
sizes  and  concentrations  of  raindrops. 

2.1  Differential  Reflectivity 

The  differential  reflectivity  (ZDR) 
provides  an  estimate  of  mean  hydrometeor 
shape.  It  is  defined  as: 

ZDR  =  10  log(ZH/ZV)  (1) 

where  ZH  and  ZV  are  the  radar 
reflectivities  measured  at  horizontal  and 
vertical  polarizations  respectively.  For 
small  raindrops  or  tumbling  ice  particles, 
ZH  and  ZV  are  equal  and  ZDR  is  zero.  The 
theoretical  values  of  ZDR  for  oblate 
particles  witb  their  minor  axes  aligned  in 
the  vertical  are  plotted  in  Figure  1;  ZDR 
increases  with  greater  oblateness  and 
higher  dielectric  constant.  In  heavier 
rain  ZDR  is  positive  and  reflects  the  mean 
shape  (and  hence  the  size)  of  the 
raindrops.  The  ZDR  of  ice  is  more  complex 
(3).  Because  of  the  low  dielectric 
constant,  dry  snowflakes  have  a  ZDR  close 
to  zero,  but  wet  snowflakes  can  have  high 
positive  values.  Graupel  tends  to  tumble 
and  so  be  associated  with  a  zero  ZDR  value. 

2.2  Linear  Depolarisation  Ratio 

The  linear  depolarization  ratio,  LDR,  is  a 
measure  of  the  hydrometeor  fall  mode  and 
appears  to  be  an  excel'lent  indicator  of  wet 
ice.  LDR  is  defined  as: 

LDR  *  10  log(ZVH/ZH)  (2) 

where  ZVH  is  the  (horizontal)  cross-polar 
return  from  a  vertically  polarized 
transmitted  pulse,  and  ZH  (as  in  equation 


1)  is  Che  co-polar  (horizontal)  return  tor 
horizontally  polarized  transmission. 

A  cross  polar  return  occurs  only  when 
oblate  hydrometeors  tall  with  their  major 
or  minor  axis  at  an  angle  to  the  vertical. 
It  particles  tall  with  their  axes  aligned 
in  the  vertical  or  horizontal  then  there  is 
no  cross-polar  return  and  LOR  is  minus 
infinity.  Computations  for  randomly 
tumbling  oblate  spheroids  are  plotted  in 
Figure  2  and  are  found  to  be  consistent 
with  Chilbolton  observations  (4).  The 
values  of  LOR  rise  as  the  particles  become 
either  more  oblate  or  of  higher  dielectric 
constant.  Snowflakes  have  such  a  low 
dielectric  constant  that  even  if  they'  are 
very  oblate  their  LOR  is  below  the  antenna 
limit  of  -32dB.  Oblate  dry  hail  or  graupel 
could  have  a  value  of  up  to  -20dB  if  the 
axial  ratio  were  as  low  as  0.5,  but  LOR 
values  above  -20dB  can  only  realistically 
occur  for  wet  tumbling  ice  particles.  Such 
high  values  are  restricted  to  the  melting 
layer.  Raindrops  give  rise  to  a  very  low 
cross-polar  return. 

The  basis  of  our  method  for  distinguishing 
between  different  forms  of  ice  is  the  value 
of  LOR  when  the  particles  start  to  melt  and 
become  coated  with  water.  In  stratiform 
clouds  and  some  showers  this  melting  layer 
is  associated  with  peak  LDR  values  of 
about  -15dB,  consistent  (Figure  2)  with 
snowflakes  having  a  rocking  or  spinning 
motion  and  an  axial  ratio  of  about  0.5.  In 
contrast,  in  some  vigorous  showers  which  we 
believe  contain  graupel,  the  LDR  in  the 
melting  layer  reaches  a  maximum  value  of 
only  -25dB,  which  we  Interpret  as  tumbling 
water  coated  graupel  pellets  with  an  axial 
ratio  of  about  0.8. 

2.3  Possible  Effect  of  Static  Blectric 
Fields  on  Radar  Parameters. 

Computations  suggest  that  the  fields 
required  to  trigger  lightning  discharges 
(300kV/m)  are  also  sufficient  to  change  the 
shape  of  large  raindrops  (5)  and  to  alter 
the  orientation  of  ice  crystals  (6).  This 
implies  that  the  electric  field  might  be 
sensed  by  its  effect  on  the  ZDR  of 
raindrops  and  the  ZDR  or  LDR  of  ice 
crystals.  However,  although  it  might  be 
possible  to  detect  sudden  changes  in  these 
parameters  coincident  with  a  lightning 
flash,  it  appears  impossible  to  separate 
the  component  of  ZDR  due  to  the  steady 
electric  field.  For  example,  if  a  high 
ZDR  is  associated  with  raindrops,  then, 
unless  we  have  an  Independent  measurement 
of  drop  size,  we  cannot  quantify  the 
distorting  effect  of  the  electric  field. 
Similarly  with  ice  crystals,  in  the  upper 
regions  of  clouds  where  Z  is  low,  the 
values  of  ZDR  are  often  positive  and  highly 
spatially  variable,  presumably  due  to  the 
habit  of  the  crystals;  but  again,  unless  we 
have  a  means  of  knowing  the  ZDR  oi  these 
crystals  in  the  absence  of  a  field,  we 
cannot  determine  the  component  due  to  the 
electric  field. 

Large  horizontal  components  of  the  electric 
field  could  cause  ice  crystals  to  tilt,  and 
in  theory  this  should  lead  to  an  anomalous 
value  of  LDR.  Unfortunately,  the  radar 
return  from  ice  crystals  is  usually 
restricted  to  low  Z  regions,  otherwise  it 
is  masked  by  larger  ice  particles.  For  low 
Z  regions,  the  cross-polar  return  from  dry 
ice  particles  will  generally  fall  below  the 
level  of  radar  sensitivity. 


3.  RADAR  OBSERVATIONS  OF  CLOUDS  AND 
TRIGGERED  LIQfTNING 

A  vertical  section  (RHI)  of  the  three  radar 
parameters  through  a  shower  believed  to 
contain  only  snowflakes  and  aggregates  is 
shown  in  Figure  3.  In  this  case  a 
pronounced  'bright  band'  of  enhanced 
reflectivity  in  Z  at  an  altitude  of  1.5km 
is  caused  by  snowflakes  becoming  wet  and 
reflecting  microwaves  as  if  they  were  giant 
raindrops.  This  bright  band  is  clearly 
visible  ' in  the  RHI  with  this  narrow 
beamwidth  research  radar,  but  most  radars 
scan  in  azimuth  (PPI)  to  obtain  greater 
areal  coverage,  and  in  such  a  case  the 
bright  band  in  Z  is  much  more  difficult  to 
identify.  In  Figure  3  the  oblate  melting 
snowflakes  are  associated  with  a  bright 
band  in  ZDR  with  values  reaching  2dB,  but 
these  values  of  ZDR  cannot  be  taken  as 
automatic  indication  of  the  presence  of 
melting  snowflakes.  In  this  example  the 
value  of  ZDR  in  the  rain  is  over  O.SdB, 
and  in  heavier  rain  values  up  to  5  or  6dB 
are  possible.  We  also  note  in  Figure  3  the 
slightly  positive  values  of  ZDR  above  3km 
altitude;  this  low  Z  region  presumably 
containing  aligned  high  density  ice 
crystals . 

Tl.a  presence  of  the  melting  snowflakes  is 
much  easier  to  identify  from  the  LDR  data 
in  Figure  3,  where  peak  values  of  LDR  reach 
about  -15dB.  In  contrast,  LDR  values  in 
the  rain  are  near  the  antenna  limit  of 
-32dB,  and  reach  about  -27dB  in  the  low  Z 
ice  region  above  3km  where  the  ZDR 
indicated  the  presence  of  high  density 
crystals.  The  high  values  of  LDR  very 
close  to  the  ground  in  Figure  3  are  caused 
by  ground  clutter  which  can  easily  corrupt 
the  low  power  cross-polar  return. 

A  vigorous  shower  is  depicted  in  Figure  4 
and  it  is  clear  that  the  polarization 
parameters  have  a  quite  different 
character;  the  LDR  values  in  the  bright 
band  are  much  lower  than  in  Figure  3.  we 
believe  that  this  cloud  contains  graupel. 
The  dry  tumbling  graupel  high  in  the  cloud 
gives  negligible  LDR,  but  melting  occurs  at 
about  2km  and  LDR  rises  to  about  -25dB; 
this  weak  *LDR  graupel  bright  band*  is 
consistent  with  tumbling  wet  ice  with  an 
axial  ratio  of  0.8.  In  heavy  rain  the  LDR 
is  Just  detectable  and  is  explicable  in 
terms  of  a  raindrop  canting  angle  of  about 
5  degrees.  Values  of  ZDR  are  low  for  the 
tumbling  dry  ice,  but  rise  monotonically  as 
the  graupel  melts  and  assumes  the 
equilibrium  shape  of  the  large  raindrops. 
The  Z  in  Figure  4  is  greater  than  50dBZ, 
but  the  absolute  magnitude  of  Z  cannot  be 
taken  as  an  indication  of  the  presence  of 
graupel;  values  of  Z  in  the  bright  band  in 
Figure  3  exceed  40dBZ  but  some  graupel 
showers  have  Z  values  lower  than  40dBZ. 

The  cloud  which  caused  the  triggered 
lightning  event  is  displayed  in  Figure  5. 
The  values  of  Z  are  lower  than  in  Figure  4 
and  both  types  of  LDR  bright  band  are 
present.  Between  65  and  73km  range  the 
peak  values  of  LDR  are  about  -25dB  and 
suggest  the  presence  of  melting  graupel. 
but  outside  these  ranges  the  higher  values 
of  LDR  in  the  melting  layer  are  indicative 
of  snowflakes  or  aggregates.  At  these 
ranges,  if  LDR  is  low,  the  cross-polar 
power  is  below  the  level  of  detectability 
by  the  radar. 


The  Meteorological  Office  C-130  aircraft 
was  penetrating  the  cloud  in  Figure  5  at  a 
height  where  the  temperature  was  just  below 
zero,  and  was  flying  on  a  radial  towards 
the  radar  on  the  saae  azimuth  as  the  radar 
scan.  As  the  aircraft  penetrated  the 
region  where  we  infer  graupel  was  present, 
it  was  3truck  by  lightning,  even  though  the 
value  of  Z  was  only  about  40dBZ.  A  value 
of  50dBZ  is  often  assumed  to  be  required 
for  lightning.  Data  from  the  UK 
Electricity  Supply  Industry  Lightning  Flash 
Locator  (Personal  Communication.  L  J 
Scott)  shows  that  the  first  cloud  to  ground 
lightning  over  Southern  England  was 
detected  from  this  cloud  six  minutes  after 
the  strike  to  the  aircraft.  We  conclude 
that  the  aircraft  triggered  the  lightning 
flash  and  that  it  did  not  discharge  to 
ground.  In  the  next  Section  we  analyze  the 
aircraft  measnrements  of  particle  type. 

4.  AIRCRAFT  OBSERVATIONS  OF  ICS  PARTICLE 
TYPE 

The  Meteorological  Office  C-130  aircraft 
was  equipped  with  a  PMS  2D  precipitation 
probe;  an  optoelectronic  device  consisting 
of  a  laser  light  source,  focusing  optics,  a 
linear  32  element  photodiode  array  and  high 
speed  electronics  to  record  the  images  of 
the  precipitation  particles.  The 
photodiode  separation  is  0.2mm,  and  the 
resolution  in  the  direction  of  the  aircraft 
flight  is  dependent  on  aircraft  speed  but 
was  also  about  0.2mm.  Manual  recognition 
of  the  many  images  is  tedious,  error  prone 
and  subjective,  so  an  automatic  system  of 
classifying  the  shapes  following  Ouroure 
(7)  was  implemented. 

There  are  several  stages  in  This  analysis 
system.  First.  the  centre  of  mass  of  each 
images  was  found,  and  then  the  radius  of 
the  image  a  function  of  angle  was 
calculated.  The  mean  radius,  the  axial 
ratio  and  orientation  could  then  be 
deduced.  Using  a  fast  Fourier  transform 
the  power  spectra  of  the  radial b  was 
calculated  and  normalised  with  respect  to 
the  mean  radius  of  the  image;  this  spectra 
13  a  size  Independent  measure  of  the 
periodicities  in  the  edge  of  the  particle 
Image.  Classification  involves  a 
comparison  of  the  normalised  power  in  the 
second.  sixth,  and  all  the  higher 
harmonics.  Aggregates  have  more  power  in 
the  higher  harmonics,  but  graupel  particles 
are  recognised  by  their  smoother  profile. 

As  a  check  on  the  reliability  of  the  Image 
analysis  programs.  the  value  of  Z  was 
calculated  from  the  aircraft  images, 
assuming  the  particles  to  be  dry  ice,  and 
compared  with  the  Z  observed  by  the  radar. 
The  comparison  is  plotted  in  Figure  6, 
where  the  solid  squares  are  the  values 
computed  from  the  images  and  the  solid  line 
is  a  polynomial  fit  through  these  points. 
Because  the  values  of  Z  were  changing 
rapidly  with  height,  and  the  scan  was 
separated  by  two  minutes  from  the  aircraft 
lightning  strike,  there  is  some  uncertainty 
in  the  precise  values  of  Z  at  the  aircraft 
position.  Accordingly,  the  vertical  lines 
in  the  Figure  represent  the  range  of  values 
of  Z  for  the  data  point  at  the  aircraft 
flight  level  and  the  radar  observation  300m 
above  and  below  this  height.  Good 
agreement  la  obtained.  Confirmation  of  the 
aircraft  position  was  obtained  by  a  direct 
'hit'  on  the  previous  REI  radar  scan. 


Figures  7  and  3  show  images  during  this 
penetration  which  the  analysis  program  has 
classified  as  graupel  and  aggregates.  The 
shape  of  the  images  has  been  adjusted  in 
the  computations  to  account  for  any  changes 
in  aircraft  speed.  These  Figures  confirm 
that  the  classif ication  algorithm  is 
operating  satisfactorily.  Figure  9 
compares  the  predicted  values  of  LDR  when 
these  ice  particles  become  wet  with  the 
actual  values  of  LDR  observed  by  the  rada: 
at  the  melting  layer.  The  axial  ratio  of 
each  particle  was  determined,  and  then 
assuming  that  they  were  randomly  tumbling 
and  wet,  a  reflectivity  weighted  value  of 
LDR  was  computed.  These  values  are  plotted 
as  the  open  squares;  the  bold  line  is  a  fit 
through  these  squares .  The  observed  values 
of  LDR  in  the  melting  layer  are  represented 
by  the  solid  squares  in  the  Pigure;  we  note 
that  the  overall  observed  variation  of  LDR 
with  distance  agrees  well  with  the  values 
predicted  from  the  images.  The  oblate 
aggregates  predominate  where  LDR  is  high, 
whereas  the  more  spherical  graupel 
particles  are  more  common  where  LDR  is 
lower.  The  actual  observe,  values  of  LDR 
are  consistently  lower  than  the 
computations,  but  this  is  consistent 
(Figure  2)  with  the  ice  particles  becoming 
slightly  less  oblate  as  they  acquire  a 
water  coating. 

5.  DISCUSSION  AND  CONCLUSION 

Aircraft  and  radar  comparisons  confirm  that 
the  LDR  value  in  the  melting  layer  can  be 
used  to  differentiate  regions  containing 
graupel  from  those  where  snowflakes  are 
found.  The  triggered  lightning  occurred  in 
a  region  where  graupel  was  predominant.  It 
is  widely  believed-  (8)  that  charge 
separation  in  clouds  occurs  when  small  ice 
crystals  collide  with  and  separate  from 
graupel  particles;  where  snowflakes  are 
present  the  ice  crystals  will  be  collected. 
Accordingly,  we  would  expect  appreciable 
electric  fields  to  be  limited  to  regions  of 
graupel . 

There  are  two  possible  reasons  why  the 
presence  of  the  aircraft  could  raise  the 
electric  field  to  a  level  needed  to  trigger 
lightning.  Firstly,  the  presence  of  the 
conducting  aircraft  itself,  and  secondly, 
the  charge  acquired  by  the  aircraft  as  a 
result  of  collisions  with  ice  crystals.  The 
precise  mechanism  is  unclear.  Collisions  of 
the  aircraft  with  ice  crystals  generally 
charge  up  the  aircraft  negatively  (9)  which 
would  appear  to  make  the  initiation  of 
streamers  and  subsequent  triggering  of 
lightning  more  difficult  (10). 

The  use  of  the  linear  depolarization  ratio 
(LDR)  for  detecting  clouds  where  there  is  a 
danger  of  lightning  appears  to  have  great 
potential.  The  results  reported  here  have 
been  made  at  S-band  with  a  very  narrow 
beamwldth  radar;  further  research  is  needed 
before  Implementing  the  technique  for 
C-band  radars  with  smaller  antennas. 
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Figure  1.  ZDR  values  as  a  function  of 
axial  ratio  for  various  precipitation 
particles.  assuming  the  major  axis  is 
horizontal . 


Figure  2.  LOR  values  for  randomly  tumbling 
particles  as  a  function  of  axial  ratio. 
Raindrops  do  not  tumble,  but  the  'water* 
curve  is  applicable  to  wet  ice  particles. 
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Figure  3.  An  RHI  scan  of  the  radar  reflectivity  (Z),  the  differential  reflectivity  ( ZDR ) 
and  Linear  Depolarization  Ratio  (LDR)  through  a  shower.  The  high  values  of  LDR  at  l . 5km 
altitude  are  interpreted  as  due  to  melting  snowflakes. 
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RHI  SCAN  ON  13/07/88  AT  155423  UT 
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Figure  4.  An  RHI  scan  of  Z,  ZDR  and  LDR  through  a  shower  believed  to  contain  graupel  or 
soft  hail  pellets.  Values  of  LDR  peak  at  about  -25dB  at  l . 5km  altitude  as  the  graupel 
starts  to  melt. 
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RHI  SCAN  ON  13/07/88  AT  130435  UT 

TAPE  5130  RASTER  28  SCAN  1  A Z  269.50deg 
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Figure  5.  An  RHI  scan  through  the  cloud  which  triggered  the  lightning  flash  to  the 
aircraft.  The  lightning  occurred  at  70km  range  where  the  LOR  indicates  the  presence  of 
graupel . 


Figure  7.  Particles  identified  as  graupel  Figure  9.  A  comparison  of  the  U)R  values 
by  the  analysis  package  (shape  uncorrected  observed  by  the  radar  in  the  melting  layer, 
for  airspeed).  with  those  c Deputed  from  the  aircraft  probe 

images. 


